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ABSTRACT 

A  quantitative  study  is  made  of  the  effect  of  successive 
correlated  binary  collisions  on  the  transport  properties  of  gases. 

In  particular,  the  triple  collision  transport  integrals  determining 
the  first  density  coefficients  of  thermal  conductivity  and  viscosity 
are  derived  for  a  gas  of  hard  spherical  molecules  and  estimates  are 
presented  for  these  integrals.  By  applying  the  method  to  a  two 
dimensional  gas  of  hard  disks,  a  logarithmic  density  dependence  of 
the  transport  coefficients  is  demonstrated. 

J  ‘  1  i 

In  addition,  an  analysis  is  made  of  experimental  data  for 
the  transport  properties  as  a  function  of  density.  It  is  shown  that 
the  theoretically  predicted  density  dependence  is  at  least  consistent 
with  the  experimental  information.  Experimental  values  for  the  first 
density  coefficients  of  thermal  conductivity  and  viscosity  are  reported, 
together  with  an  assessment  of  their  precision. 
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CHAPTER  I  , 

PURPOSE  OF  THE  INVESTIGATION  ^ 

1.1  INTRODUCTION 


Transport  properties  of  gases  are  directly  related  to  the  transfer 

of'  energy  and  momentum  between  colliding  molecules.  In  dilute  non-, ionized 

gases  the  transport  properties  are  independent  of  the  density.  Roughly 

speaking  the  number  of  molecules  available  for  the  transport  of  thermal 

energy  and  momentum  is  proportional  to  the  density,  but  the  distance  over 

which  this  transport  is  effected,  the  mean  free  path,  is  in  a  dilute  gas 

inversely  proportional  to  the  density. 

The  theoretical  description  of  the  transport  coefficients  of  a  dilute 

gas  is  based  on  the  Boltzmann  equation,  which  describes  the  rate  of  change 

of  the  molecular  velocity  distribution  as  a  result  of  binary  collisions 

between  the  molecules .  Based  on  the  Boltzmann  equation  the  Chapman-Enskog 

theory,  provides  a  reliable  procedure  to  obtain  the  transport  coefficients  of 
1  2 

dilute  gases.  *  This  theory  is  sufficiently  well  understood  that  in  turn 

experimental  transport  coefficients  of  dilute  gases  are  used  to  obtain  infor- 

3 

mation  on  the  molecular  Interaction  potential. 

Experience  shows  that  the  transport  coefficients  do  vary  with  the 
density,  when  a  gas  is  compressed.  From  a  molecular  point  of  view  this 
density  dependence  has  two  origins.  First,  the  effect  of  triple  and  higher 
order  collisions  on  the  transport  properties  can  no  longer  be  neglected. 


+)  For  a  survey  of  this  report  the  reader  is  referred  to  the  introductory 
Sections  1.1,  1.2,  3.1  and  4.1  and  to  the  Sections  2.8,  3.6  and  4.5  which 
summarize  the  results  (see  table  of  contents). 
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Secondly,  the  difference  in  position  of  molecules  involved  in  a  collision  - 
(spatial  inhomogeneity)  leads  to  an  additional  density  dependent  contribution, 
sometimes  referred  to  as  collisional  transfer. 

Prediction  of  the  density  dependence  of  the  transport  properties  is 
a  task  for  statistical  mechanics.  For  this  purpose  the  Boltzmann  equatibn 
has  to  be  generalized  to  include  the  effects  mentioned  above.  Up  to  the 
present  approximate  theories  have  been  considered  only.  The  oldest  attempt 
is  the  theory  of  Enskog  for  hard  spherical  molecules.  In  this  theory  the 
spatial  correlation  between  the  molecules  is  approximated  by  the  equilibrium 
radial  distribution  function  and  the  velocity  correlation  is  disregarded. 
Recently,  several  approximate  theories  have  been  considered'  to  describe  the 
first  density  correction  to  the  transport  coefficients  of  gases. ^  All 
these  theories  attempt  to  describe  the  transport  coefficients  in  terms  of 
binary  collisions  and  an  explicit  evaluation  of  the  multiple  collision  dynamics 
is  avoided.  As  a  result  it  is  difficult  to  assess  the  reliability  and 
predictive  power  of  such  descriptions. 

It  is  well  known  that  the  equation  of  state  of  a  gas  can  be  represented 
by  a  power  series  in  the  density  p,  the  virial  expansion: 

pV  =  RT  [1  +  Bp  +  Cp2  +  .  .  .]  .  (1.1-1) 

The  various  virial  coefficients  in  this  expansion  can  be  related  to  molecular 
cluster  integrals  which  are  integrals  over  the  configurations  of  n  simultaneously 
interacting  molecules.  Thus  the  characteristic  size  determining  these  cluster 
integrals  is  the  range  of  the  molecular  interaction. 
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.  To  describe  the  density  dependence  of  the  transport  coefficients  a 

8 

similar  power  series  expansion  was  proposed  by  Bogoliubov.  To  discuss  the  effect 
of  multiple  collisions  we  can  distinguish  between  two  types  of  multiple  collisions: 

1.  Genuine  multiple  collisions  in  which  all  n  molecules  are  interacting 
simultaneously.  It  is  evident  that  such  collisions  contain  again  the  size 
of  the  molecular  interaction  range  as  a  characteristic  length. 

2.  Multiple  collisions  that  consist  of  a  succession  of  correlated  binary 
collisions  between  n  particles.  After  a  binary  collision  the  velocities  of 
two  molecules  will  remain  correlated  until  at  least  the  next  collision. 

Thus  one  can  expect  that  the  characteristic  length  associated  with  successive 
correlated  binary  collisions  is  no  longer  the  interaction  range,  but  the 
mean  free  path.  However,  in  contrast  to  the  interaction  range,  the  mean 
free  path  is  itself  dependent  on  the  density.  As  a  result  it  is  no  longer 
evident  that  the  effect  of  such  multiple  collisions  can  be  represented  by  a 
power  series  in  the  density.  Actually,  as  discussed  in  Chapter  III,  it  turns 
out  that  these  multiple  collisions  lead  to  a  representation  of  the  transport 
coefficients  which  contains  terms  logarithmic  in  the  density  p.  In  particular, 
for  the  thermal  conductivity  A  and  the  viscosity  n  we  obtain 

2  2 

X  =  +  A^p  +  AjP  &np  +  Agp  +  .  .  . 

2  2  (1.1-2) 

n  =  n  +  tup  +  n«p  £np  +  tup  +  .  .  . 

O  1  Z  i 

To  predict  the  transport  properties  of  a  moderately  dense  gas  it  becomes  our 
task  to  develop  theories  and  procedures  by  which  the  various  coefficients  A^ 
and  can  be  computed  for  a  given  molecular  interaction  potential. 
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1.2  SCOPE  OF  THIS  RESEARCH 

The  present  report  describes  the  first  results  of  an  effort 
to  obtain  qualitative  estimates  of  the  effect  of  multiple  collisions,  in 
particular  of  triple  collisions,  on  the  transport  properties.  Evaluation  of 
the  triple  collision  effects  is  needed  to  predict  the  first  density  coefficients 
A^  and  ti  in  (1.1-2).  The  analysis  is  highly  desirable  for  a  number  of  reasons: 

(1)  Reliable  predictions  of  the  density  dependence  of  transport  properties 
should  be  based  on  the  rigorous  statistical  mechanical  equations  which  contain 
the  dynamics  of  multiple  collisions. 

(2)  Quantitative  information  on  the  effect  of  multiple  collisions,  at  least 
for  some  molecular  models,  is  needed  to  make  an  assessment  of  the  validity  of 
any  approximate  theory. 

(3)  In  order  to  study  the  effect  of  nonanalytic  terms  in  the  density  expansion 
for  the  transport  coefficients  a  quantitative  estimate  of  the  multiple  collision 
effects  is  essential. 

(4)  As  will  be  discussed  in  Chapter  IV,  the  first  density  coefficient  of  viscosity, 
in  (1.1-2),  becomes  negative  at  high  temperatures.  The  sign  of  the  corresponding 

coefficient  A^  for  the  thermal  conductivity  at  high  temperatures  is  not  definitely 
clear.  These  first  density  coefficients  A^  and  r,^  are  the  sum  of  two  contributions: 
a  collisional  transfer  contribution  (related  to  binary  collisions)  which  is  usually 
positive  and  a  triple  collision  contribution  which  is  negative.  Thus  even  to 
predict  the  sign  of  the  first  density  correction  at  high  temperatures  in  addition 
to  the  collisional  transfer  contribution  a  quantitative  estimate  of  the  triple 
collision  term  is  very  important. 

Since  a  study  of  the  effect  of  triple  collisions  turns  out  to  be  rather 
complicated,  we  limit  ourselves  in  this  report  to  a  study  of  the  effects  for  a 
very  simple  molecular  model,  namely  two  and  three  dimensional  hard  spheres. 
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The  transport  properties  describe  the  linear  response  of  the  gas  to  the 
presence  of  gradients.  The  dilute  gas  values  and  no  in  (1.1-2)  are  obtained 
from  the  solution  of  the  linearized  Boltzmann  equation.  As  a  result  they  can 
be  expressed  in  terms  of  binary  collision  integrals.  (For  hard  spheres  these 
results  are  summarized  in  Section  2.2).  For  repulsive  intermolecular  potentials 
the  first  density  corrections  are  in  principle  to  be  determined  by  solving  the 
Choh-Uhlenbeck  Integral  equation  which  is  a  generalization  of  the  linear  Boltzmann 
equation  to  include  the  effect  of  triple  collisions.  In  Chapter  II  we  develop 
a  procedure  to  solve  this  equation  for  a  gas  of  hard  spheres.  The  first  density 
corrections  and  are  related  to  well  defined  triple  collision  transport 

integrals  and  preliminary  values  for  these  triple  collision  integrals  are 
presented.  The  results  are  discussed  in  Section  2.8  and  compared  with  the 
approximate  predictions  of  the  theory  of  Enskog. 

To  study  the  possible  non-analytic  density  dependence  of  the  transport 
coefficients,  we  consider  in  Chapter  III  a  two  dimensional  model  gas  of  hard 
disks.  The  advantage  of  this  model  is  that  the  logarithmic  density  dependence 
can  be  evaluated  with  the  same  techniques  as  developed  in  Chapter  II.  The 
coefficient  of  the  term  logarithmic  in  the  density  is  evaluated  explicitly  for 
this  model  and  turns  out  to  be  non  vanishing. 

Finally,  in  Chapter  IV  we  analyze  experimental  data  for  the  transport 
coefficients  of  a  number  of  gases.  In  this  chapter  we  investigate  whether  the 
theoretically  proposed  density  function  (1.1-2)  is  consistent  with  the  experimental 
data.  Based  on  the  analysis,  first  density  coefficients  of  thermal  conductivity 
and  viscosity,  together  with  an  assessment  of  their  precision,  are  reported  as 
a  function  of  temperature. 
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CHAPTER  II 

TRIPLE  COLLISION  CONTRIBUTION  TO  THE  TRANSPORT  COEFFICIENTS 
FOR  A  GAS  OF  HARD  SPHERES 

2.1  INTRODUCTION 


In  the  kinetic  theory  of  gases  the  transport  coefficients  are  usually 
related  to  the  velocity  distribution  function  which  specifies  the  state  of 
the  gas.  This  velocity  distribution  function  f(r,v;t)  is  defined  such  that 
f(r,v;t)dvdr  determines  the  number  of  molecules  at  a  time  t  with  position 
coordinates  in  the  range  dr  around  ?  and  with  velocities  in  the  range  dv"  around 
v.  The  number  density  n  and  the  mean  velocity  u  at  r  and  t  are  related  to 
this  distribution  function  by 

n(r;t)  =  /  f(r,v;t)  dv  ,  (2.1-1) 

nO?;t)  u(r;t)  =  f  vf (r,v;t)  dv  .  (2.1-2) 

^  Jb  afe 

Denoting  the  thermal  velocity  of  a  molecule  by  V  -  v  -  u,  the  average  local 
kinetic  energy  is 

~  n(r;t)  kT(r;t)  =  f  V2  f(r,v;t)  dv  ,  (2.1-3) 

1  2-* 

where  m  is  the  mass  of  the  molecules.  Since  each  molecule  contributes  ymV  V 
to  the  flow  of  thermal  energy,  the  total  heat  flux  vector  reads 

Q  =  |  m  /  V2  V  f(r,v;t)  dv  .  (2.1-4) 

The  transport  of  momentum  -  or  pressure  tensor  reads 

^  =  m  /  V  V  ffr’.vjt)  dv  ,  (2.1-5) 

In  the  equilibrium  state  the  distribution  function  reduces  to  the  Maxwell- 
Boltzmann  distribution 

2 

f(v)  - "  (2sr>3/2exp[-  Is1  •  (2-1-6) 
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In  equilibrium  the  heat  flux  vector  (2.1-4)  vanishes  and  the  pressure 
tensor  becomes  diagonal 


(2.1-7) 


where  I  is  the  unit  tensor  and  p  the  hydrostatic  pressure. 

When  the  gas  is  slightly  out  of  equilibrium  the  heat  flux  vector 
and  the  off  diagonal  elements  of  the  pressure  tensor  are  proportional. to 
the  gradients  in  temperature  T  and  fluid  velocity  ”u,  respectively: 


9T 

q  =  -AVT  =  -A  , 
o 

"A  =t  9  " 

P  =  pl-2rig»u  , 


(2.1-8) 

(2.1-9) 


where  |* IS  is  the  symmetric  traceless  tensor  derived  from  the  dyadic  |»-'u.1 
These  linear  relations  define  the  coefficient  of  thermal  conductivity  X  and 
the  coefficient  of  shear  viscosity  .  On  comparing  these  definitions 
with  the  molecular  expressions  (2.1-4)  and  (2.1-5)  for  the  fluxes  the 
transport  coefficients  are  related  to  the  distribution  function  f.  For 
this  purpose  f  is  expanded  in  terms  of  the  gradients  of  the  macroscopic 
variables  and  only  terms  linear  in  the  gradients  are  retained.  In  zeroth 


approximation  f  reduces  to  the  local  Maxwell  distribution 

3/2 

w  ~  uj  i 
2kT  i  » 


m  \  '  _ r  m(v  -  ’ll)2 


Eo  ■  n  (W>  expl- 


(2.1-10) 


where  the  macroscopic  variables  still  depend  on  r  and  t.  in  first  approximation 

f  -  f  (1  +  *)  ,  (2.1-11) 


+)  In  this  report  we  do  not  consider  the  bulk  viscosity,  since  both  its 
dilute  gas  value  and  its  first  density  coefficient  vanish. 
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where  <t  is  linear  in  the  gradients  of  T  and  u.  For  the  time  scale  under 
consideration  f,  and  consequently  $,  does  not  vary  explicitly  with  time, 
but  only  via  the  dependence  on  the  macroscopic  variables  nCr'jt),  “u(r;t) 
and  T(r;t): 


3f  .  3f  9n  3f  3u  3f  3T  f 

Tt  3n  St  +  3u  >t  +  3T  9t  *  (2.1-12) 

From  the  fact  that  $,  like  f  itself,  is  a  scalar  function,  one 


concludes  that  $  must  have  the  form 

31nT 


*(V)  =  -A(V)  -  B(V)  :|»"u  , 


(2.1-13) 


where  A(V)  and  B(V)  are  scalar  functions  of  the  molecular  velocity  V.  To 
ascertain  that  n,  "u  and  T  in  f  keep  satisfying  their  definitions  (2.1-1), 
(2.1-2)  and  (2.1-3)  in  the  presence  of  gradients,  the  function  $  is  subject 
to  the  auxiliary  conditions 


/  f  0  dv  =  0 
o 

/  Vf  $  dv  =  0 
o 

I  V2f  *  dv  =  0 
o 


Specifically,  (2.1-15)  implies  that 


(2.1-14) 

(2.1-15) 

(2.1-16) 


/  f  V  A(V)  d v  =  0 
o 


(2.1-17) 

When  we  substitute  this  form  of  the  distribution  function  into  the 
molecular  expressions  (2.1-4)  and  (2.1-5)  for  the  fluxes  and  compare  the 
result  with  the  definitions  (2.1-8)  and  (2.1-9)  of  the  transport  coefficients, 
we  obtain 


X  =  |  /  dV  fQ(V)  (—•  -  |)  V2A(V)  ,  (2.1-18) 

n  =  /  dV  fQ(V)  V°V:V°V  B(V)  .  (2.1-19) 

Thus  it  becomes  our  task  to  assess  the  effect  of  the  triple  collisions  on 
the  functions  A(V)  and  B(V). 
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2.2  TRANSPORT  COEFFICIENTS'  IN  THE  DILUTE  STATE 

The  first  density  correction  to  the  transport  coefficients  is 

not  only  caused  by  a  triple  collision  contribution  to  the  kinetic  fluxes 

(2.1-4)  and  (2.1-5),  but  also  by  an  additional  collisional  transfer 

contribution  from  binary  collisions.  The  latter  term  results  from  the 

difference  in  position  of  two  molecules  Involved  in  a  binary  collision. 

For  a  gas  of  hard  spheres  the  theory  of  Enskog  presents  an  estimate  for 

4 

both  contributions.  It  turns  out  that  the  collisional  transfer  contribution 

9 

to  the  first  density  contribution  as  evaluated  by  Enskog  is  exact.  There¬ 
fore,  only  the  effect  of  triple  collisions  remains  to  be  evaluated  and  it 
is  sufficient  for  our  present  considerations  to  limit  ourselves  to  the 
so  called  spatially  homogeneous  case.  That  is,  the  distribution  function 
f  is  considered  as  a  function  of  v"  and  t  and  the  gradients  of  the  macroscopic 
variables,  but  its  explicit  dependence  on  the  position  ^  is  disregarded. 

Before  discussing  the  effect  of  triple  collisions  we  first  summarize 
how  the  transport  coefficients  of  a  dilute  gas  are  obtained  when  only  binary 
collisions  are  taken  into  account.  This  will  enable  us  to  elucidate  how  the 
triple  collision  transport  integrals  are  generalizations  of  the  well  known 
binary  collision  transport  integrals. 

In  a  dilute  gas  the  time  variation  of  f  is  described  by  the  Boltzmann 

equation 

3f (v^  ; t) 

- - - J(ff)  ,  (2.2-1) 

where  the  collision  term  J(ff)  represents  the  effect  of  binary  collisions. 

For  a  gas  of  hard  spherical  molecules  this  collision  term  J(ff)  is  given  by^ 
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J(ff)  =  /dVg/dk  a2(v“21.k)  [f  (v.[;t)f  (v2;t)  -  f  (v^jtjf  (v2;t)  ]  .  (2.2-2) 

Here  a  is  the  diameter  of  the  molecules,  v^  the  relative  velocity  v2  -  v^ 

and  k  the  perihelion  vector  specifying  the  geometry  of  a  binary  collision 

.a 

between  molecules  1  and  2.  For  hard  spheres  this  perihelion  vector  k  is 
the  unit  vector  in  the  direction  from  the  center  of  particle  2  to  the 
center  of  particle  1  at  the  time  of  collision. 


Figure  1.  The  perihelion  vector  k  of  a  binary  collision. 

The  distribution  function  f(v^;t)  changes  with  time,  since  during 
a  time  interval  dt  certain  molecules  with  velocity  v^  will  collide  with 
molecules  with  velocity  v2»  The  number  of  such  collisions  with  perihelion 
vector  k  will  be  proportional  to  f (v^;t)f (v2;t)  and  the  differential  cross 
section  cj2(v2- .k)dk.  On  the  other  hand,  molecules  with  velocities  "v^  will 
collide  with  molecules  v'  such  that  after  the  collision  their  velocities  will 
become  and  Thus  v^  and  v2  are  the  initial  velocities  for  a 


binary  collision  with  perihelion  vector  k  and  final  velocities  v^  and  v2< 
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FUTURE 


PRESENT 


PAST 


V/  =  £  (12)  V, 


=  <£(12)7, 


Figure  2.  Schematic  representation  of  a  binary  collision.  The  lines  represent 
the  particle  trajectories  and  the  circle  encloses  a  region 
in  which  a  binary  collision  occurs .+^ 

It -is  convenient  to  relate  the  velocities  "v^  to  the  velocities'^  by  an  operator  ^(12) : 

=  cf(12)v1  ;  7*  =^(12)^  .  (2.2-3) 


+)  '  Using  the  concept  of  inverse  collisions  one  usually  considers  "v^,  v^  as  the 
initial  velocities  for  a  binary  collision  with  a  perihelion  vector  -k  and  final 
velocities  v^,  v^ •  Since  a  generalization  of  the  concept  of  inverse  collisions 

g 

to  triple  collisions  is  not  obvious,  we  use  here  the  Bogoliubov  form  of  the 

i 

Boltzmann  equation,  in  which  the  consideration  of  inverse  encounters  is  avoided. 
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In  general  the  substitution  operator  $ (1.  .  .s)  associated  with  s  particles 
_8 

s 


is  defined  as 


eftl.  •  -s)  =  lim  s  (1.  .  .3)  n  S  (i)  , 
t  -*■  «  i=l 

‘where  S.  (1.  .  .s)  transforms  the  coordinates  and  the  velocities  of  s 


(2.2-4) 


particles  into  those  at  a  time  t  later,  if  the  s  particles  move  under  their 
mutual  interaction.  The  effect  of  the  operator  cf  (12)  on"v^  and”^ 
schematically  represented  by  Fig.  2.  If  the  positions  and  the  velocities  of 
the  particles  are  such  that  they  did  not  collide  in  the  past,  the  operator 
J*  (12)  reduces  to  one.  The  particles  traverse  the  trajectories  indicated 
in  Fig.  2  from  bottom  to  top.  However,  since  the  present  Is  determined 
by  the  past,  we  actually  retrace  the  history  of  the  particles  in  the  past+\ 
that  is,  we  read  the  diagrams  always  from  top  to  bottom.  The  effect  of  the 
substitution  operators  related  to  triple  collisions  will  be  indicated  later 
by  similar  diagrams.  In  terms  of  the  substitution  operator  & (12)  the 
Boltzmann  equation  can  be  written 


of  (v^t) 

at 


o2/dv2/dk  (v21  .k)  T  (12)f  (v^;  t)f  (v“2;  t)  , 


with 


(2.2-5) 


T(12)  =  £ (12)  -  1  .  (2.2-6) 

For  a  derivation  of  the  transport  coefficients  it  is  sufficient  to 
consider  the  linearized  Boltzmann  equation.  That  is,  we  approximate  f  by 
fQ(l  +  $)  as  in  (2.1-11),  eliminate  the  time  derivatives  of  n,  u  and  T  in 


+)  This  can  be  circumvented  by  changing  the  direction  of  time  as  will  be 
done  in  a  later  stage. 
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(2.1-12)  using  the  Euler  equations  and  retain  only  terms  linear  in  the 
gradients.^-  The  resulting  linearized  Boltzmann  equation  which  can  be 
found  in  many  places  in  the  literature,  reads 

V*<V>  ■  -  !>  +  Sr  ?i°V  ^  * 


(2.2-7) 


where  I 2  is  a  linearized  binary  collision  integral  operator 

I2(*(Vi))  =  -a2/dV2/dk(v21.k)  f0(V1)  fo(V2)  T(12)t<5(V1)  +.*(^)]  .  (2.2-8) 

Since  the  total  kinetic  energy  is  invariant  under  a  binary  collision,  the 

operator  T(12)  commutes  with  fQ(V^)fo(V2) . 

Using  the  specific  form  (2.1-13)  for  the  perturbation  function  <5  we 

obtain  linearized  integral  equations  for  the  functions  A(vp  and  B(V1) 

.2 


^  /  \  mV  _  c  ^ 

Io(V-rA'°'  (Vn ))  =  fQ(V1)  (-g±  -  §)  Vx  * 


L2v*r 


2KT 


(2.2-9) 


IjCvbv/0’  (Vj))  -  f,^)  Sy  . 


where  the  superscript  (o)  indicates  that  we  consider  only  the  contribution 
of  zeroth  order  in  the  density  here. 

To  solve  these  integral  equations  one  customarily  approximates  the 


solution  by  a  finite  sum  of  Sonine  polynomials  S 


(m) 


A  <»!>  =  -  32— I  \  00  S  wp  , 


no 


k=l 

1/2  N-l 
m 


B  (°^  /v  )  =  ^  ~~  (  m  )  l  b^°^(N)  S^V(W2) 

CV  16  2',irkTJ  kT  /  Dk  w  &5/2lV  * 
no  k=o 

where  we  have  introduced  a  dimensionless  velocity 

m  i/2  ^ 

y  =  (.  m.  A  V 

1  ^kT'  i  ’ 


(2.2-10) 


(2.2-11) 
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The  Sonine  polynomial  (x)  is  defined  as  the  coefficient  of  tm  in  the 
expansion 


l  S^(x)  tm  =  — - — ; ■  .  exp{-  -- }■. 

(l-t)n+1  1_t 


(2.2-12) 


m=o 


They, become  the  eigenfunctions  of  the  linearized  Boltzmann  operator 
when  the  molecular  interaction  potential  obeys  an  inverse  fourth  power  law 
(Maxwell  molecules).  The  number  N  of  Sonine  polynomials  taken  into  account 
in  the  approximate  solution  (2.2-10)  is  referred  to  as  the  order  of  the 
Enskog  approximation.  If  we  substitute  the  solution  (2.2-10)  into  the 
expressions  (2.1-18)  and  (2.1-19)  for  the  transport  coefficients,  we  obtain 

3m  1/2 


A  =  A  (N)  =  (j£I)  '  a<°)(N)  , 

o  o  - .  I  mr  1  * 

6  4a 

1/2 

n  =  n  (N)  =  — ^  (2^)  b<o)(N)  . 

O  0  ,  ,  IT  o 


(2.2-13) 


16a  ‘ 


(o) 


(o) 


The  values  of  the  coefficients  a^  } (N)  and  b^  J  (N)  are  determined  by 
2 

a  variational  method  ;  they  should  satisfy  a  set  of  linear  equations 


I  a£0)(N)  =  641  (i  =  1,  .  .  ,N) 

k“l 


(2.2-14) 


b'f  -  \0  (2  -  0.  .  .  .  B-l)  , 

k=o 

where  6  is  the  Kronecker  delta  and  where  and  b^4^  are  matrix  elements 

representing  inner  products  associated  with  the  linearized  Boltzmann  operator : 
1 


(°)  — — kT —  /dW  if  S^^(W^)  I  (M  S^^(W^)) 

k£  2/2tT  2  2  "'  ll  !,3/2l  r  2H  3/21  1;;  * 

mn  a 


(2.2-15) 


,<o) 

k£ 


~  -4V;d* i  *°A  s5/>i> :  • 

mn  a 
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The  equations  have  been  normalized  such  that  for  hard  spheres  a|°^ (1)  = 
b<o)(l)  =  a|^  =  b^  =  1.  The  operator  I2  is  a  symmetric  operator, 
so  that 


(o) 
£k  * 


(2.2-16) 


For  hard  spheres  the  matrix  elements  can  be  readily  evaluated. 
Using  the  generating  function  (2.2-12)  for  (W^) ,  the  elements  a^^  and 
b5°^  are  the  coefficients  of  s^t^  in  the  expansions'^ 


(1  -  |(s  +  t)}1/2st(l-st)“3-  i{l  -  -|(s  +  t)}"1/2s2t2(l  -  st)-2, 

(2.2-17) 

{1  -  y(s  +  t)}1/2(l  +  |st  +  is2t2)(l  -  st)-4-  |{1  -  ~(s  +  t)}”1/2st(l  -  st)-3 
-  |g{l  -  |(s  +  t)r3/2s2t2(l  -  st)-2. 

The  coefficients  up  to  the  third  Enskog  approximation  are  listed  in  Table  I 

and  the  corresponding  corrections  a|°^and  b^0^  to  the  transport  coefficients 

in  Table  II.  For  Maxwell  molecules  the  first  Enskog  approximation  would 

suffice.  From  Table  II  it  is  evident  that  the  expansion  still  converges 

fast  for  hard  spheres;  the  second  Enskog  approximation  amounts  to  only  a 

few  percents.  For  more  realistic  molecular  interactions  the  convergence 
2 

is  even  faster. 


00  QO 

r  I  a<°>  sV 

k-1  l-l  k* 


W  CO  , 

E  I  b<°}  s  t* 

k=o  £=o 
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Table  I 


The  matrix  elements  a^°^  and  b£°\  defined  by  (2.2-15)  for  hard 


spheres 


(o) 

1 

Co) 

1 

a12 

“  “  4 

a13 

32 

(o) 

45 

(o) 

103 

a22 

16 

a23 

=  “  128 

(o) 

5657 

a33 

1024 

b(0) 

1 

b(o) 

1 

01 

"  ”  4 

02 

32 

V  (o) 
bll 

205 
-  48 

rs 

O  CM 

& 

163 

128 

b(o) 

11889 

22 

1024 
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Table  II 


The  successive  approximations  (N)  and  b^0^ (N)  determining  the 


solution  of  the  linearized  Boltzmann  equation  for  hard  spheres 


First  Enskog  approximation:  N  =  1 


a{o)(l)  =  1 


b(0> 


(1)  =  1 


Second  Enskog  approximation:  N  =  2 


1  + 


3 

202 


6 

101 


Third  Enskog  approximation:  N  =  3 


(37)2 
88  x  7439 


b(o)(3)  =  1  + 
o 


(347)2 

808  x  145043 


(°)  (3)  _  + - 

2  K  '  11  11  x  7439 


6  42681 

101  101  x  145043 


148 

7439 


1388 

145043 
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2.3  THE  THEORY  OF  ENSKOG 


The  transport  coefficients  \  and  n  deduced  from  the  Boltzmann 

o  o 

equation  are  independent  of  the  density.  A  first  attempt  to  generalize 
the  Boltzmann  equation  to  account  for  the  density  dependence  of  the 

4 

transport  coefficients  was  made  by  Enskog.  Strictly  speaking,  the  number  of 

molecules  experiencing  a  binary  collision  with  perihelion  vector  k  will 

be  proportional  to  the  pair  distribution  function  f ^ (r^ ^r^-ck.v^jt) . 

(2) 

In  equilibrium  the  pair  distribution  function  f^  '  is  related  to  the  one 

(1) 

particle  distribution  function  f  =  f  by 

o  o  1 

fo2)  »*i  “  a^»^)  -  g(<0  f0(vi)f0(v2)  »  (2.3-1) 

where  g(a)  is  the  equilibrium  radial  distribution  function  at  a  separation 
a.  This  equilibrium  radial  distribution  function  is  independent  of  the 
velocities  and  depends  only  on  the  relative  distance  r2^(  =  o  at  a  collision). 
Enskog  assumed  the  same  relationship  between  the  pair  distribution  function 

9 

and  f,  when  the  system  is  not  in  equilibrium  : 

f^Cr^Vi.r'i  -  ofe,v2;t)  =  g(a)f (v^;t)f (v2;t)  ,  (2.3-2) 

where  g(a)  is  now  the  local  equilibrium  radial  distribution  function 
depending  on  the  local  density  and  temperature.  Thus,  again  disgarding  the 
spatial  dependence  of  the  velocity  distribution  function  f  in  our  present 
discussion,  Enskog  generalized  the  collision  term  J(ff)  in  (2.2-1)  to 

JE(ff)  =  a2g(a)/dv2/d‘k(v^1.k)  T(12)f(^i;t)f(^2;t)  =  g(a)J(ff)  .  (2.3-3) 

If  we  consider  again  a  perturbation  linear  in  the  gradients 

f  =  fo(l  +  *)  ,  (2.3-4) 
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the  linearized  Boltzmann  equation  (2.2-7)  is  to  be  replaced  with 


g(o)I2($  (Vx))  =  -  fQ(V1) 


•"  mV2 

/  1  5^  ^  91nT  ,  m  3 

( - j)  v  —  +  T-  V  V  :  g=r-  u 

2kT  2  1  3r  kT  1  1  3rl 


(2.3-5) 


E  — *•  —1  -A 

This  equation  is  identical  to  (2.2-7),  if  we  identify  $  (V^)  with  g  (o)$(V^), 
Consequently,  the  theory  of  Enskog  predicts  that  the  (kinetic)  transport 
of  a  dense  gas  of  hard  spheres  are  related  to  the  dilute  gas  values  by 


E  ^o 

A  =  g(a)  '  n  "  g(a) 


E  \ 


(2.3-6) 


The  value  g(o)  of  the  radial  distribution  function  is  directly 
related  to  the  equation  of  state 


p  =  nkT{l  +  bpg(a)}  ,  -  (2.3-  ) 

where 

bp  =  ~  no3  t  (2.3-8) 

is  the  covolume  of  the  spheres.  The  pressure  p,  as  well  as  g(a)  can  be 
expanded  into  a  power  series  in  terms  of  the  density 

g(o)  =  1  +  bp  +  .  .  .  .  (2.3-9) 

The  coefficient  of  the  first  density  dependent  correction  term  in  g(a)  is 
determined  by  the  so-called  excluded  volume  of  a  third  particle;  that  is, 
the  volume  where  a  third  particle  would  overlap  with  the  two  particles  1 
and  2  simultaneously. 
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Figure  3.  The  excluded  volume  determining  the  first 
density  correction  to  g(c). 

The  power  series  (2.3-9)  implies  a  power  series  for  the  transport 
coefficients  (2.3-6): 

A  =  Aq  +  A^bp  +  •  •  •  p 

(2.3-10) 

n  =  n  +  n-,bp  +  .  .  . 

o  1 

In  the  sequel  we  shall  refer  to  A^  and  A^  as  the  first  density  coefficient 
of  thermal  conductivity,  resp.  viscosity.  From  (2.3-6)  and  (2.3-9)  it  is 
evident  that  in  the  Enskog  theory  A^  and  are  approximated  by 

>  (2-3'u) 

independent  of  the  order  of  the  Enskog  approximation  in  Aq  and  nQ- 
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The  Enskog  theory  was  based  on  an  ad  hoc  assumption  (2.3-2).  A 

systematic  procedure  to  take  into  account  successive  higher  order  con- 

8 

tributions  in  the  density  was  proposed  by  Bogoliubov.  For  the  first 

12 

density  correction  the  theory  was  worked  out  by  Choh  and  Uhlenbeck 

13  14 

and  further  elucidated  by  Green  and  Cohen.  Limiting  ourselves  to  the 
first  density  correction,  the  Boltzmann  equation  (2.2-1)  is  to  be  replaced 
with 

Sf (v.  ;t) 

- ^ -  =  J(ff)  +  K(fff)  ,  (2.4-1) 

where  K(fff)  represents  the  effects  of  triple  collisions  and  is  related 
to  the  product  of  three  one  particle  distribution  functions: 

K(fff)  =  /dv2dv3/dr2dr3012T(12;3)f  (v^sOf  (v2;t)f(^3;t)  .  (2.4-2) 

12 

Here  012  is  an  differential  operator  introduced  by  Choh  and  Uhlenbeck 
_  1  **12  3_  ,  1  **12  l 

°12  m  3r^  ‘  3v!^  m  3r^  ’  3v“2  *  (2.4-2) 

where  <j>32  is  the  intermolecular  potential  between  particles  1  and  2.  For 

hard  spheres  <J)^2  is  a  step  function  and  consequently  032  contains  a  6-function. 

13  14 

The  operator  x(12;3)  is  a  combination  of  substitution  operators  ’ 

x (12 ; 3)  =  ef(123)  -  $ (12)  ^(13)  -  J*(12)J*(23)  +«f(12)  ,  (2.4-4) 

where  ^(123),  defined  by  (2.2-4),  is  a  generalization  of  the  two  particle 
substitution  operator  cf(12)  and  includes  the  interaction  between  three 
particles.  Apart  from  neglecting  collisions  between  three  particles  that 
actually  occur,  the  Boltzmann  equation  took  into  account  some  binary 
collisions  that  actually  will  not  occur,  due  to  the  interference  of  a 
third  particle.  Thus  in  x(12;3)  we  subtract  certain  combinations  of  two 
particle  streaming  operators. 
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To  derive  the  transport  coefficients  we  again  linearize  the 


equation  by  substituting  f  =  fQ(l  +4).  As  a  result  the  linearized 


Boltzmann  equation  (2.2-7)  is  replaced  with 
I2($(V^))  +  bpl3(*(^_))  =  -  fo(v1)  |  (■ 


tariff  +!iVi41tf]>(2-4-5) 


where  1^  is  a  linearized  ternary  collision  integral  operator 


(2.4-6) 

I3(<t(Vk1))  =  ’  ip/dv2dv3/dr2di3ei2T(12;3)fo(V1)fo(V2)fo(V3)  [*(^1)+*(^)+*(V3)] 

For  convenience  a  factor  (bp)  ^  has  been  entered  in  the  definition  of  I^, 
so  that  1^  has  the  same  dimension  as  l2>  From  the  right  hand  side  of 
(2.4-5)  it  is  evident  that  $(V^)  will  have  again  the  form  (2.1-13).  To 
solve  the  equation  we  expand  the  functions  A(V^)  and  B(V^)  in  terms  of 
the  density: 


A(Vl)  =  A(o)(V1)  +  bpA(1)(V1)  , 

B(VX)  =  B(°)(V1)  +  bpB(1)(V1)  . 


(2.4-7) 


The  corresponding  expansion  for  the  transport  coefficients  was  introduced 
in  (2.3-10).  From  the  expressions  (2.1-18)  and  (2.1-19)  we  see  that  the 
first  density  coefficients  A^  and  are  given  by 


A 


n 


1 


1 


k  1  5  2  (1^ 

-!/dvi£o<v<-2kT  -!>vA  <v  • 

"  To/dV0(V  :  vJv(1)(v1)  . 


(2.4-8) 
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The  integral  equations  for  (V^)  and  B^(V^)  are  derived  by 
substituting  (2.4-7)  into  (2.1-13)  and  (2.4-5).  Since  (V^)  and 
B^°^(V^)  satisfy  the  linearized  Boltzmann  equation  (2.2-9),  we  obtain 


i2(v^1a(1)(v1))  =  -  i3(v1a(o)(v1))  , 
I3^lV<1)(Vl))  =  -I3^iViB<o)(V1))  ,  • 


(2.4-9) 


neglecting  terms  quadratic  in  the  density.  In  these  equations,  to  which 
we  refer  as  Choh-Uhlenbeck  integral  equation,  the  new  ternary  collision 
operator  I3  operates  on  the  functions  A^°^  (V^)  and  B^(V-)  that  are  known 
from  the  solution  of  the  linearized  Boltzmann  equation.  The  kernel  of  the 
Choh-Uhlenbeck  integral  equation  is  the  same  as  that  of  the  linearized 
Boltzmann  equation.  In  principle,  one  could  solve  this  equation  by  expanding 
A^  ^  and  B^  ^  again  in  terms  of  Sonine  polynomials.^-"*  However,  we  are  not 
interested  in  A^  ^  and  B^^  as  such,  but  only  in  the  moments  (2.4-8)  of  these 
functions,  which  can  be  written  explicitly  in  terms  of  A^  and  B^  +\ 

For  this  purpose  we  substitute  (2.2-9)  into  (2.4-8) 

h  =  i  fd\  v(1)(vi)*I2(V<o)(vi))  » 

(2.4-10) 

rl  =  ll  /d¥?i  B(1)(VX)  :  I2(V°Vk1B(0)(V1))  . 


+)  The  author  is  indebted  to  Dr.  M.  H.  Ernst  for  this  remark. 


23 


AEDC-TR-69-68 


Since  I2  is  a  symmetric  operator,  we  can  interchange  and  A^,  resp. 

and  in  (2.4-10).  Thus  it  follows  from  (2.4-9)  that  X^  and 

are  related  to  A^°^  and  B^°^  by 

(2.4-11) 

By  substituting  the  Sonine  polynomial  expansion  (2.2-10)  for  A^0^ 

and  B^°^  (with  known  coefficients)  into  (2.4-11)  we  can  evaluate  X^  and 

in  successive  Enskog  approximations.  Introducing  again  dimensionless 

1/  2 _ ^ 

velocities  =  (m/2kT)  V^,  we  obtain  in  the  Nth  Enskog  approximation 


Xn(N)  N  N  .  v  ..  /.  v 

— i. _ =  — E  E  CN)  a  °  (N)  a/^  • 

xoa)  Z=1  \  w  aa  w  “w  ’ 


(2.4-12) 


n.  (N)  N-l  N-l  ..  ,  . 

00  b‘  00  V  ■ 

The  new  matrix  elements  a^^  and  b^^  are  inner  products  associated  with 
the  linearized  triple  collision  integral  operator  1^: 


(1) _ 1  kT 


(k)  /tt2> 


-Sr2/dWlWlS3/2(ffl)-I30,l83/iCWl>)  * 
2v2ir  mn  o 


(2.4-13) 


(D  =  1  kT 

ki  _  rz —  2  2 

2/2tt  mn  o 


/dwx  wJifjsS*  (wj)  :i3  (wjw^s^^wj)) 


These  matrix  elements  are  the  Immediate  generalizations  of  the  matrix  elements 
a^°^  and  b^°\  defined  by  (2.2-15),  that  determine  the  dilute  gas  values 
Xq  and  dq;  they  were  first  introduced  in  a  previous  publication.^  For 
convenience,  we  have  divided  the  first  density  coefficients  X^(N)  and  t]^(N) 
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by  the  dilute. gas  values  Aq(1)  and  tiq(1) »  so  that  the  ratios  in  (2.4-12) 

are  Independent  of  the  temperature. 

It  is  convenient  to  symmetrize  the  integrands  in  (2.4-13).  The 

16 


operator  can  be  symmetrized  as  shown  by  Green 


3 


!-(*(?,»  =  -  /dV9dvl0(123)f  (V  )f  (V  )f  (V  )  Z  4 >(V  )  (2.4-14) 

3  1  2bp  2  3  o  1  o  2  o  3  n=1  n  > 

where  0(123)  is  a  configurational  triple  collision  integral  operator 


0(123)  =  /dr2dr3[012T(12;3)  +  613t(13;2)  +  623t(23;1)]  . 


(2.4-15) 


The  integrands  in  (2.4-13)  can  be  further  symmetrized  by  interchanging  the 
velocity  variables.  Thus  we  obtain 

1/2  v 


,(1) 

k£ 


Vs  (|r)i,"/dW1dW2dW3exp  [-(wj  +  W2  +  W32)] 


32tt  a 


•  \  V3/2  » 

m=l  n=l 


(2.4-16) 


,(1)  _  1 

bk£ - 


1/2 


32ir  a 


%-j  (|T)  ZdW1dW2dW3  exp[-(W1  +  W2  +  W3)] 


•  K  t?A)Wj:oa22'>  l%s5/2(wn> » 


m=l 


n=l 


where  we  used  the  property  that  0(123)  commutes  with  the  kinetic  energy 
and  therefore  with  the  product  of  the  Maxwell  distributions.^ 

There  exists  a  simple  relationship  between  the  rigorous  expressions 
(2.4-11)  for  A^  and  and  the  approximate  expressions  (2.3-11)  given  by  the 

theory  of  Enskog.  The  operator  I3>  defined  in  (2.4-6),  contains  an  integration 
over  all  possible  positions  of  particle  3.  If  we  limit  this  Integration 
to  those  values  of  ?3  for  which  both  r33<a  and  r23<a*  that  i8«  the  excluded 
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volume  o£  Fig.  3t  the  operator  reduces  to  as  shown  In  a  previous 
paper. ^  Thus  the  rigorous  expressions  (2.4-11)  reduce  to  the  values 
predicted  by  Enskog,  when  the  configurational  Integral  Is  restricted  to 
those  configurations  that  determine  In  equilibrium  the  first  dentlty 
correction  to  the  radial  distribution,  function  g(a). 

2.5  CLASSIFICATION  OF  THE  TRIPLE  COLLISION  EVENTS 


The  transport  coefficients  and  nQ  In  the  dilute  state  were 
determined  by  the  binary  collision  integrals  (2.2-15).  To  obtain  the 
first  density  corrections  and  we  need  to  evaluate  the  triple 
collision  integrals  (2.4-16).  Apart  from  the  fact  that  the  dynamics  of 
triple  collisions  is  more  complicated  than  that  of  an  individual  binary 
collision,  there  is  another  major  difference  between  the  two  sets  of 
transport  collision  integrals.  In  the  binary  collision  integral  operator 
I 2»  defined  in  (2.2-8),  the  configurational  integration  does  not  extend 
over  all  possible  positions  particle  2,  but  the  integration  is 

restricted  to  a  two  dimensional  surface,  specified  by  the  perihelion 

M  Jk 

vector  k.  For  a  given  k  the  operator  T(12)  in  (2.2-8)  can  be  evaluated 
at  any  point  along  the  trajectories  of  the  two  particles  (provided 
that  the  two  particles  are  receding) ,  since  one  needs  only  the  velocities 
before  and  after  the  binary  collision.  On  the  other  hand,  the  triple 
collision  integral  operator  0(123),  defined  in  (2.4-15)  is  still  a  six 
dimensional  integral  and  contains  the  differential  operators  0^ ;  it  has 
not  yet  been  expressed  in  terms  of  the  velocities  of  the  particles  before 
and  after  certain  collision  events.  It  is  possible  to  convert  the  six 
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dimensional  triple  collision  integral  (2.4-15)  into  a  five  dimensional 
surface  integral  which  is  much  more,  similar  to  the  two  dimensional 
surface  integral  form  of  the  binary  collision  integral  operator  ^ 

A  derivation  of  the  surface  integral  form  of  the  triple  collision  integral 
operator  is  presented,  in  Appendix  A. 

It  turns  out  that  the  triple  collision  integral  can  be  decomposed 
into  a  sum  of  six  surface  integral  operators,  each  of  which  is  associated 
with  a  particular  triple  collision  event.  To  elucidate  this  result  we 
give  first  a  description  of  the  six  triple  collision  events  which  are 
schematically  represented'  in  Fig.  4.  The  positions  of  the  particles  (on 
the  five  dimensional  surface)  over  which  the  integration  has  to  be  carried 
out  are  indicated  by  the  black  dots.  The  direction  of  the  velocities  is 
indicated  by  the  arrows.  The  lines  represent  trajectories  and  are  labeled 
with  the  corresponding  particle  number.  The  circles  represent  (binary) 
collisions  between  the  particles  whose  trajectories  are  enclosed.  The  - 
operators  in  (2.4-15)  transform  the  present  phases  of  the  particles  into 
the  phases  that  the  particles  would  have  had  in  the  previous  past 
according  to  a  particular  dynamical  history.  Therefore,  we  retrace  the 
motion  backwards  in  time  and  we  shall  always  read  the  diagrams  from  top 
to  bottom,  just  as  for  the  binary  collision  diagram  of  Fig.  2.  The 
sequence  of  successive  binary  collisions  will  be  ordered  as  they  are 

encountered  in  this  procedure . 

13 

Following  Green  we  distinguish  between  real  and  hypothetical 
collisions.  Any  binary  collision  encountered  when  we  retrace  the  trajectories 
of  the  particles  taking  into  account  the  interaction  between  all  three 
particles  is  called  a  real-  collision.  Any  binary  collision  that  is  encountered 
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only  when  the  effect  of  some  other  collision  on  the  trajectories  of  the 
particles  is  disregarded  is  called  a  hypothetical  collision.  In  addition, 
we  distinguish  between  interacting  and  non-interacting  collisions.  In  a 
non-interacting  collision  the  particles  are  aimed  to  collide,  but  we  con¬ 
tinue  the  free  trajectories  of  the  particles  without  taking  into  account 
the  interaction  between  the  particles.  In  fig.  4,  non-interacting  collisions 
are  indicated  by  a  dotted  circle.  It  is  evident  from  the  definition  that 
all  collisions  encountered  in  the  diagrams  of  Fig.  4  after  a  non-interacting 
collision  are  hypothetical  collisions. 

With  this  nomenclature  the  triple  collision  events  of  Fig.  4  can 
be  described  as  follows  (always  reading  the  diagrams  from  top  to  bottom!). 

R1  (Recollision  of  the  first  kind) : 

A  real  collision  between  1  and  2  is  followed  by  a  real  collision 
between  1  and  3,  followed  by  a  real  collision  between  1  and  2  again. 

R2  (Recollision  of  the  second  kind) : 

A  real  non-interacting  collision  between  1  and  2  is  followed  by  a 
hypothetical  interacting  collision  between  1  and  3,  followed  by  a 
hypothetical  collision  between  1  and  2; 

Cl  (Cyclic  collision  of  the  first  kind): 

A  real  collision  between  1  and  2  is  followed  by  a  real  collision 
between  1  and  3,  followed  by  a  real  collision  between  2  and  3. 

C2  (Cyclic  collision  of  the  second  kind): 

A  real  non-interacting  collision  between  1  and  2  is  followed  by  a 
hypothetical  interacting  collision  between  1  and  3,  followed  by  a 
hypothetical  collision  between  2  and  3. 

HI  (Hypothetical  cyclic  collision  of  the  first  kind): 

A  real  collision  between  1  and  2  is  followed  by  a  real  non-interacting 
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collision  between  1  and  3,  followed  by  a  hypothetical  collision 
between  2  and  3. 

H2  (Hypothetical  cyclic  collision  of  the  second  kind) : 

A  real  non-interacting  collision  between  1  and  2  is  followed  by 
a  hypothetical  non-interacting  collision  between  1  and  3,  followed 
by  a  hypothetical  collision  between  2  and  3. 

As  shown  in  Appendix  A  the  configurational  triple  collision  operator 
0(123),  introduced  in  (2.4-14),  can  be  replaced  with 


0(123)  =  6  •  0(12;13)  , 
with  0(12; 13)  given  by  (A-30) 


(2.5-1) 


0(12;  13)  =  aVdljdk^  dr  £  1/12;  13)  . 


(2.5-2) 


The  points  on  the  five  dimensional  surface  are  specified  by  the  variables 

•A  4 

-k^,  k£  and  x;  k^  is  the  perihelion  vector  of  the  first  collision  (between 
1  and  2) ,  k ^  the  perihelion  vector  of  the  second  collision  (between  1  and 
3)  and  t  the  time  between  the  first  and  the  second  collision.  V^(l)  is  the 
relative  velocity  of  3  with  respect  to  1  between  the  first  and  the  second 
collision.  The  summation  over  y  in  (2.5-2)  has  to  be  taken  over  the  diagrams 
y  =  Rl,  R2,  Cl,  C2,  HI,  H2.  The  operator  T^(12;13)  is  only  different  from 
zero  for  phase  points  whose  dynamical  history  can  be  represented  by  the 


diagram  y.  For  such  phase  points  the  definition  of  T^(12;13)  is  given  in 
Table  III. 
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Table  III 

The  triple  collision  operators  T  (12; 13) 


u 

T  (12;13) 

R1 

«f(123)  -  S(  12)  #(13) 

R2 

-j(13)cf(12)  +  #(13) 

Cl 

#(123)  -  #(12)  #13) 

C2 

~#13)  # (23)  +#(13) 

HI 

-/(12)  # (23)  +  #(  12) 

H2 

# [23)  -  1 

These  expressions  for  the  operators  T^(12;13)  in  terms  of  the  streaming 
operators  J* are  presented  in  Fig.  4  below  each  diagram;  they  are  direct 
generalizations  of  the  operator  T(12),  defined  in  (2.2-6),  for  the  binary 
collision  term.  Operating  on  any  function  of  the  velocities,  T^(12;13) 
transforms  that  function  into  the  difference  between  its  value  for  the 
final  velocities  at  the  bottom  of  the  diagram  and  its  value  for  the  velocities 
between  the  second  and  the  third  collision.  The  effect  of  the  operator  has 

the  opposite  sign,  when  the  number  of  non-interacting  collisions  is  odd. 

17  18 

Note  that  in  R1  and  Cl  a  fourth  collision  may  occur;  ’  in  that  case  the 
trajectories  have  to  be  retraced  until  this  fourth  collision  has  also  been 
completed.  We  consider  only  those  diagrams  for  which  the  first  collision 
refers  to  particles  1  and  2  and  the  second  collision  to  particles  1  and  3. 

Any  other  diagram  obtained  by  permutation  of  the  particle  numbers  is 
accounted  for  by  the  factor  6  in  (2.5-1). 

In  the  previous  Section  we  pointed  out  that  the  expression  for  and 
ri^  reduce  to  the  predictions  of  Enskog,  if  the  configurational  integral  is 
restricted  to  the  excluded  volume  of  particle  3.  To  ascertain  which  triple 
collision  events  are  taken  into  account  in  the  Enskog  theory,  we  consider 
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diagram  H2  in  some  detail.  In  this  diagram  the  first  collision  (between 
1  and  2)  is  a  non-interacting  collision.  Such  a  non-interacting  collision 
implies  that  during  a  certain  time  interval  particles  1  and  2  move  inside 
each  other,  when  their  free  trajectories  are  followed  from  top  to  bottom. 

We  call  such  configurations  overlapping  configurations.  When  the  time 
between  the  first  and  the  second  collision  is  sufficiently  small,  particle 
3  will  be  aimed  to  collide  with  1,  while  1  and  2  are  still  overlapping. 

Since  the  collision  between  1  and  3  is  also  non-interacting,  1  and  3  overlap 
during  a  short  time  interval  after  the  second  collision.  When  in  addition 
the  time  between  the  second  and  the  third  collision  is  small,  it  is  possible 
that  3  will  collide  with  2,  while  particles  1  and  2  as  well  as  particles  1 
and  3  are  still  overlapping.  Thus,  such  a  collision  process  leads  to  the 
excluded  volume  configuration  (Fig.  3  if  we  interchange  1  and  3)  considered 
in  the  Enskog  theory.  Therefore,  as  mentioned  in  a  previous  paper, ^  it  can 
be  concluded  that  the  Enskog  theory  neglects  all  triple  collision  events 
shown  in  Fig.  4,  except  for  some  collisions  of  the  type  H2.  Among  all 
collisions  H2,  the  theory  of  Enskog  includes  only  those  for  which  the  time 
between  the  first  and  the  last  collision  is  so  small  that,  when  3  starts  to 
collide  with  2,  1  and  2  as  well  as  1  and  3  are  still  overlapping. 

The  surface  integral  form  (2.5-1)  for  0(123)  should  be  substituted 

into  the  expressions  (2.4-16)  for  and  b^J^ .  In  doing  so  we  take  the 

diameter  a  of  the  spheres  as  the  unit  of  length.  Since  the  velocity  was 

1/2 

expressed  in  units  (2kT/m)  ,  the  corresponding  dimensionless  time  x*  is 


x* 


1  ,2kT*/2 

—  ( - )  x 

a  m 


(2.5-3) 
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Thus  the  integrals  (2.4-16)  for  and  become  in  dimensionless 


form 


(1)  3/2"  «*  *  * 

>x/  —  rjk  T.T  Jf.T  Jl.7 


k£ 


16n 


-6/dW1dW2dW3exp[-(W‘  +  W*  +  W3>] 


.  E  W  sJ2(W*)..  0*(12;13)  EWS^o(W^)  , 
.  m  3/2  m  ,  n  J/z  n  1 


m=l 


m  3  vT  «*  ■* 


3 

E 

n=l 


(2.5-4) 


k£ 


16tt 


-6/dW1dW2dW3exp[-(W‘  +  +  W3)] 


.  I«°w  s*>(»  : 

.mm  5/2  m 
m=l 


)  :  0*(12;13) 


I«°WS  <«(w2) 
,  n  n  5/2  n 
n=l 


with 

0*(12;13)  =  E/d^d^/dT*  W31(l).k2  Ty(12;13)  .  (2.5-5) 

y 

For  the  evaluation  of  these  triple  collision  integrals  it  is  more  con¬ 
venient  to  describe  the  collisions  between  the  molecules  in  terms  of  forward 
motion,  rather  than  retracing  the  motion  backwards  in  time.  Such  a  description 
is  obtained  if  we  change  the  sign  of  all  velocities.  The  integrands  are 
symmetric  in  the  velocities  and  therefore  invariant  under  this  transformation, 
provided  that  we  replace  t  with  -t  in  the  definition  (2.2-4)  of  the 
operators.  For  the  diagrams  in  Fig.  4  it  means  that  we  change  the  direction 
of  the  velocities  as  shown  in  Fig.  5,  but  still  read  the  diagrams  from  top 
to  bottom.  In  the  sequel  of  this  chapter  we  shall  refer  to  Fig.  5  in  dis¬ 
cussing  the  triple  collision  events. 
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2.6  ANALYSIS  OF  T1IE  TRIPLE  COLLISION  TRANSPORT  INTEGRALS 


In  the  preceding  Section  the.  triple  collision  integral  was 
decomposed  into  a  sum  of  integrals  associated  with  the  six  triple  collision 
events  of  Fig.  5.  However,  there  exists  a  close  relationship  between  the 
diagrams  of  the  first  kind  (Rl,  Cl  and  HI)  and  those  of  the  second  kind 
(R2,  C2,  H2).  If  in  Rl,  Cl  or  HI  the  velocities  after  the  first  collision 
are  taken  as  the  initial  velocities,  we  obtain  the  diagrams  of  the  second 
kind.  Minor  exceptions  to  this  correspondence  will  be  discussed  later. 

Using  this  relationship  we  cat}  combine  the  contributions  from  the  diagrams 
of  the  first  kind  with  those  of  the  second  kind.^ 

As  an  example  let  us  consider  the  recollision  events  Rl  and  R2 
in  Fig.  5.  To  specify  the  effect  of  the  operators,  we  indicate  the 

velocities  after  the  i-th  collision  in  Rl  by  ^(i)  (m  a  1,2,3).  Thus  the 

a 

initial  velocities  are  W  (0),  the  velocities  after  the  first  collision  W  (1), 

m  m 

the  velocities  after  the  second  collision  If  (2),  the  velocities  after  the 

m 

third  collision  W  (3„)  and  after  a  fourth  collision,  if  that  occurs,  W  (4_). 

m  R  ra  R 

The  index  R  indicates  that  the  velocities  are  related  to  the  integration 

variables  by  a  recollision  event.  With  this  notation  we  can  write  for  any 

pair  of  functions  j  and  ) 

v  m  n 


*(k)(V'TRl(12i13)  “ 

■  *(k)<SL(o)).[*(l><w,(JB))  -  ta)(w  (2))]  . 

m  n  k  n 


(2.6-1) 
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If  (as  usually)  only  three  real  collisions  occur  j  =3:  if  four  collisions 

occur  j  =  4.  The  velocities  in  R2  are  indicated  by  the  same  notation, 

except  that  the  initial  velocities  are  no  longer  W  (0)  but  are  identified 

m 

with  W  (1) .  Thus 
m 

4>(k)(Wm).TR2(12;13)  $(£)(Wn)  = 


(2.6-2) 

=  -  ^(k)(Wm(l)).[$(£)(Wn(3R))  -  $(£)(Wn(2))]  . 

In  order  to  combine  these  two  terms  we  refer  to  the  coordinate 
system  of  Fig.  6  for  a  description  of  the  collision  process  Rl.  Particle  1 
is  taken  at  rest  at  the  origin  after  the  first  collision.  The  Z-axis  is 
taken  in  the  direction  of  W21(l)  =  W2(l)  -  W^l).  The  X-axis  is  chosen 

such  that  the  perihelion  vector  k^  of  the  first  collision  is  in  the  XZ 

a 

plane  as  indicated  in  the  figure;  0^  is  the  angle  between  -k^  and  the 
Z-axis.  The  perihelion  vector  k2  of  the  second  collision  lies  in  general 
outside  the  XZ  plane.  In  this  coordinate  system  particle  2  came  originally 
from  A.  At  the  time  of  the  first  collision  particle  2  is  at  B  and  particle 
3  at  C  with  relative  position  vectors^ 

=  -k^  ;  r“31(l)  =  -k2  -  T*i?31(l); 


r32(l) 


t*w31(i)  . 


(2.6-3) 


In  order  for  a  recollision  to  occur,  particle  3  should  collide  with  1  at  D 
after  a  time  x*  such  that  particle  1  will  collide  again  with  2. 


+) 


In  this  report  W. 


ij 


and  r 


ij 
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In  the  collision  event  R2  the  first  collision  is  a  non-interacting 
collision.  Let  the  perihelion  vector  of  this  non-interacting  collision  be 
k^,  and  the  time  between  the  first  and  the  second  collision  x*'.  This 
implies  that  in  Fig.  6  particle  2  came  originally  from  A'.  The  relative 
position  vectors  at  the  time  that  the  first  collision  begins  are  now 
indicated  by  r^  (0)  .  Thus 

^(0)  =  -k^  ;  ?31(0)  =*  -k2  -  T^W^Cl)  ; 


CO)  =  £[  - 


(2.6-4) 


We  require  that  particle  3  collides  with  1  at  D  such  that  1  and  2  will 
collide  again.  The  same  event  as  that  described  earlier  is  obtained  if 
k^  and  x*'  are  chosen  such  that+^ 

k|  =  k x  -  2  W^dJ.kj  5  t*'  =  x*  +  XQ*  ,  (2.6-5) 


with 


(2.6-6) 


This  implies  for  the  initial  positions  (0)  in  (2.6-4) 

A 

r21(0)  =  -kx  +  2  W21  (!).€,_  ;  r^CO)  =  -k2  -(x*  +  t*)W31(1)  5 


(2.6-7) 


r32  (  ®  =  (0)  -r01  (0) 


-31 


21' 


We  take  k^  and  x*,  defined  by  (2.6-5)  as  the  new  integration  variables  for 

the  integrals  associated  with  the  diagrams  of  the  second  kind.  Then  the 

velocities  W  (i)  in  (2.6-2)  are  the  same  functions  of  the  integration  variables 
m 


1  -A 

+)  In  this  report  a  superscript  indicates  a  unit  vector:  a  =  a/a. 
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as  those  in  (2.6-1)  and  one  does  not  need  to  retrace  the  diagrams  R1 
and  R2  separately. 

The  integrands  associated  with  the  other  diagrams  of  the  first  and 

second  kind  can  be  combined  by  the  same  procedure.  Using  the  velocities 

W  (1)  after  the  first  collision  as  integration  variables  we  obtain 
m 


ak£5  =  "  ~6fd\  (1)d^2  d)d”3 a)  eXp  [-(Wl(1)  +  W2(1)  +  W3(1)] 
16it 

^  ^  ^  3  3 

.  /dk1dk2;dT*|W21(l).k1|w31(l).k2  E  E 


(2.6-8) 


l"SV2<<(1»-  tVV8J/*‘,£‘3R»  -  V»S3/»»  + 


+  «m(0>S3/2»m(0))-  '  "n(2)S3/X<2))}  +  <C1) 

-  "„(1)S3/2(W»(1))-  {V3C>S3/X(V>  ‘  V^X*2”  >  +  (C2> 

-  \<°>S3/X<°»-  («„<VS352<Wn(V>  '  VX/X'1”*  +  (H1) 

+  V1)S3/X(1»'  fVVSvX(V>  -  V1)S3/2(Wn<1))}1  .  (H2) 

and  a  similar  expression  for  b^^ . 

To  evaluate  the  various  terms  in  the  integrand  we  only  have  to  retrace 
three  triple  collision  events  to  which  we  refer  as  recollisions,  cyclic 
collisions  and  hypothetical  cyclic  collisions.  Nevertheless,  we  specify  in 
(2.6-8)  whether  a  term  is  associated  with  a  diagram  yl  of  the  first  kind  or 
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y 2  of  the  second  kind.  The  reason  is  that  there  are  some  differences  in 
the  region  of  integration  over  the  variables.  For  Instance,  in  C2  and 
H2  it  is  possible  that  3  starts  to  collide  with  1  while  1  and  2  are  still 
overlapping.  This  phenomenon  occurs  when  0<x*'<x*o.  For  such  a  phase 
point  there  exists  no  diagram  of  the  first  kind  as  a  counter  part.  Thus 

&  A 

it  follows  from  (2.6-5)  that  negative  values  of  x*  such  that  -tq<t  <0 
are  to  be  included  in  the  contribution  from  C2  and  H2,  but  not  in  any  of 
the  other  terms.  Apart  from  the  phase  points  leading  to  overlapping  con¬ 
figurations,  the  integration  region  associated  with  a  diagram  of  the  second 
kind  is  very  similar  to  that  associated  with  the  corresponding  diagram  of 
the  first  kind.  The  remaining  differences  arise  from  the  effect  of  various 
possible  fourth  collisions  between  the  three  particles,  the  effect  of  which 
is  small. 

The  integration  region  will  be  specified  by  formulating  the  conditions 
that  the  integration  variables  should  satisfy.  Some  of  these  conditions 
should  be  imposed  on  all  integrals,  namely 


W21(1)*kl<0  * 

w3iU).k2>0  . 

Condition  (2.6-9)  guarantees  that  the  first  collision  and  (2.6-10) 
second  collision  will  take  place.  Furthermore  as  noted  above 

T*  >  _x*  (y  =  C2,  H2) 


(2.6-9) 
(2.6-10) 
that  the 


(2.6-11) 
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and 

t*>0  (y  -  Rl,  R2,  Cl,  HI)  .  (2.6-12) 

(A  negative  value  of  t*  Is  not  allowed  for  R2,  since  It  implies  that  the 
third  collision  would  precede  the  second  collision) . 

The  requirement  that  a  third  collision  should  take  place  leads  to  a 
second  class  of  conditions.  They  differ,  of  course,  whether  we  consider  a 
recollision,  a  cyclic  collision  or  a  hypothetical  cyclic  collision  and  are 
most  conveniently  expressed  in  terms  of  the  positions  and  velocities  at 
(just  after)  the  time  of  the  second  collision.  The  initial  positions  for  a 
diagram  of  the  first  kind  were  given  by  (2.6-3)  and  those  for  a  diagram  of 
the  second  kind  by  (2.6-7).  The  positions  at  the  time  of  the  second  collision 
are  the  same  for  all  diagrams: 

*r21(2)  =  -t±  +  t*W21(1)  ;  ?31(2)  =  -£,  ;  r^tf)  =  -  VT*"21(1)  .  (2.6-13) 


The  basic  equations  describing  a  collision  between  two  hard  spheres  are 
presented  in  Appendix  B.  Applying  the  inequalities  (B-3)  and  B-4)  we  see 
that  for  recollisions  (y  =  Rl,  R2) 


{Wn(2)  x  rn(2)}2  -  W21(2)  <  0  t  (2.6-14) 

W2i(2).?2i(2)  <  0  ,  (2.6-15) 


for  cyclic  collisions  (y  =  Cl,  C2) 

(W32(2)  x  r32(2)}2  -  W22(2)  <  0  , 

W32(2).r32(2)  <  0  , 

and  for  hypothetical  cyclic  collisions  (y  =  HI,  H2) 

{W32(l)  x  r32(2)}2  -  W32(l)  <  0  # 

W32a).r32(2)  <  0  , 


(2.6-16) 

(2.6-17) 

(2.6-18) 

(2.6-19) 
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r^CO)  >  1  (y  -  R2,  C2,  H2)  .  (2.6-23) 

For  the  cyclic  and  hypothetical  cyclic  collisions,  particles  2  and  3  should 
also  not  overlap  at  the  time  of  the  second  collision,  since  the  third  collision 
should  not  precede  the  second  one. 

r32<2^  >  1  •  (2.6-24) 

The  triple  collision  events  were  defined  in  Section  2.5  together 
with  a  distinction  between  real  and  hypothetical  collisions.  This  distinction 
leads  to  some  minor  additional  conditions: 

1)  In  order  that  the  first  collision  be  real,  the  initial  phases  in 
Fig.  5  should  be  such  that  the  particles  did  not  experience  a  collision  in 
their  past  history. 

ii)  In  order  that  in  Rl,  Cl  and  HI  the  second  collision  be  real, 
particle  3  should  not  collide  with  2  before  colliding  with  1. 

iii)  In  order  that  in  Rl  and  Cl  the  third  collision  be  real,  no  other 
collision  should  interfere  along  the  path  of  particle  2. 
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Although  these  conditions  together  guarantee  the  proper  triple 
collision  events,  some  of  them  are  redundant.  Therefore,  not  all  con¬ 
ditions  need  to  be  imposed  explicitly.  For  instance,  when  in  R1  the 
phases  of  the  particles  are  such  that  the  third  collision  (between  1  and 
2)  will  occur  it  is  no  longer  possible  for  3  to  collide  with  2  before  the 
second  collision.  On  the  other  hand,  in  Cl  this  possibility  has  to  be 
excluded  explicitly. 

We  made  an  effort  to  reduce  the  number  of  conditions  to  be  imposed. 

The  conditions  retained  are  presented  in  Table  IV.  The  conditions  (u-1) 
through  (y-5)  determine  the  basic  phenomenon.  The  remaining  conditions 
exclude  possible  fourth  collisions  that  would  interfere  with  the  triple 
collision  series  desired. 

We  note  that  the  integration  regions  associated  with  R1  and  R2  are 
identical,  except  for  the  exclusion  of  a  possible  fourth  collision  in  the 

18  20 

previous  past  ((Rl-6)  and  (R2-6)),  the  probability  of  which  is  very  small.  * 
The  related  integrands  in  (2.6-8)  only  differ  in  the  possibility  of  a 
fourth  collision  succeeding  the  third  one.  Except  for  possible  negative 
values  of  t*,  the  integrations  associated  with  Cl  and  C2,  resp.  HI  and  H2, 
are  also  very  similar. 

The  position  vectors  r  (n)  are  related  to  the  integration  variables 
by  (2.6-3),  (2.6-7)  and  (2.6-13).  The  velocities  W^  (0)  and  W^  (2)  are 

A  *4 

related  to  the  integration  variables  k^,  k2  and  W^(l)  by 

Wn(0)  =  W21(l)  -  2W21(l).k1k1  , 


w31(0)  =  w31(l)  -  W^fD.kjkj 


(2.6-25) 


W32(°)  =  W32d)  +  W^CQ.k^!  , 
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W21<2)  =W21(1)  -  W31(l).k2k2  , 

W31(2)  =  W31(l)  -  2W31(l).k2k2  ,  (2.6-26) 

W32(2)  =  W32(l)  -  W31(l).k2k2  . 

The  integrals  (2.6-8)  together  with  the  conditions  of  Table  IV 
specify  the  triple  collision  transport  integrals  to  be  evaluated. 
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Table  IV 

Conditions  specifying  the  integration  region  for  the  triple  collision 
transport  integrals. 

Recollisions  of  the  first  kind  R1 


T*  >  0 

(Rl-1) 

WyCD.k,  <  0 

CR1-2) 

-k  -» 

k^.k-  <  0 

(Rl-3) 

r*  w21C1),kl 

W31(1)>k2  >  '  t^2  ~ 

(Rl-4) 

{W21(2)  x  ?21(2))2  -  W2x(2)  <  0 

(Rl-5) 

^32C0)^32(1)  <  0  °r  {^32(0)  X  ?32(1)}2  "  W32(0)  >  ° 

(Rl-6) 

Recollisions  of  the  second  kind  R2 

i*  >  0 

(R2-1) 

W21(l>.fci  <  0 

(R2-2) 

k1.’k2  <  0 

(R2-3) 

W  d).k2  >  £  j* 

(R2-4) 

1  2 

{W21(2)  x  fn(2)}2  -  W21(2)  < 

0 

(R2-5) 

W32(l)  .?32(0)  <  0  or  {W32(l)  x  r“32(0)}2  - 

w32(i)  >  0 

(R2-6) 
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Table  IV  (continued) 
Cyclic  collisions  of  the  first  kind  Cl 


t*  >  0 

W^D.ki  <  0 
W31(l).k2  >  0 
W32(2).?32(2)  <  0 
{W32(2)  x  ^32 (2) }2  -  W22(2)  <  0 

r32(1)  >  1 

W31(0) .^31(1)  <  0  or  {W31(0)  x  ^31(1)}2-  W31(0)  >  0 

3  should  not  collide  with  2  before  colliding  with  1 
2  should  not  recollide  with  1  before  colliding  with  3 

Cyclic  collisions  of  the  second  kind  C2 


T*  >  -T* 

o 

^21(1)-^1  <  0 


W31(l).k2  >  0 


W32(2)‘r32(2)  <  0 
(W32(2)  x  r32(2)}2  -  W22(2)  <  0 


r32(0)  >  1 


r32(2)  >  1 


(l).'r32(0)  <  0  or  {W32(l)  x  r“32C0)}2  -  W2?(l)  >  0 


32’ 


(Cl-1) 

(Cl-2) 

(Cl-3) 

(Cl-4) 

(Cl-5) 

(Cl-6) 

(Cl-7) 

(Cl-8) 
(Cl— 9) 

(C2-1) 

(C2-2) 

(C2-3) 

(C2-4) 

(C2-5) 

(C2-6) 

(C2-7) 

(C2-8) 


46 


AEDC-TR-69-68 


Table  IV  (continued) 

Hypothetical  cyclic  collisions  of  the  first  kind  HI 


T*  >  0 

(Hl-1) 

•ft  «* 

W21(1),kl  < 

0 

(HI- 2) 

4 

w31(i).k2  > 

0 

(Hl-3) 

W32(l).r32(2) 

<  0 

(Hl-4) 

{W32(l)  x  r32(2)}2  - 

w32(i)  <  0 

(Hl-5) 

r32(1)  ” 

1 

(Hl-6) 

r32(2)  > 

1 

(HI- 7) 

W3i(0).r3i(1) 

<  0  or  (W31(0)  x  r*31(l)}2  -  W31(0)  >  0 

(HI- 8) 

W32(°) *r32(1) 

<  0  or  {W32(0)  x  r*32(l)}2  -  W32(0)  >  0 

(HI- 9) 

Hypothetical  cyclic  collisions  of  the  second  kind  H2 


T*  >  -  X* 

o 

(H2-1) 

•A  4 

(H2-2) 

W21(1)#kl  <  ° 

i?31  (l)  .t2  >  o 

(H2-3) 

W32(l).f32(2)  <  0 

(H2-4) 

{W32(l)  X  ?32(2)}2  -  W32(l)  <  0 

(H2-5) 

4<°>  >  1 

(H2-6) 

4(2>  >  1 

(H2-7) 
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2.7  EVALUATION  OF  THE  TRIPLE  COLLISION  TRANSPORT  INTEGRALS 
IN  THE  FIRST  ENSKOG  APPROXIMATION 


The  triple  collision  integrals  (2.6-8)  which  contain  still 
fourteenintegrations,  can  be  reduced  to  seven  dimensional  integrals. 
First,  the  operators  T^(12;13)  commute  with  the  total  momentum  and  the 
conditions  of  Table  IV# specifying  the  integration  region, are  independent 
of  the  total  momentum.  Thus  we  can  integrate  over  the  total  momentum 


or  velocity  and  express  the  resulting  integrals  in  terms  of  the  relative 


velocities  V^U)  and  W^(l).  Secondly,  the  integrand  is  a  scalar,  invariant 
under  rotations  of  the  coordinate  system.  Thus  we  can  integrate  over  all 


possible  directions  of  1^(1)  and  evaluate  the  resulting  integrand  with 
taken  in  the  Z-direction  as  in  Fig.  6.  Furthermore,  we  can  integrate  over 


all  directions  of  the  component  of  k^  in  the  XY  plane  and  take  k^  in  the  XZ 
plane  as  in  Fig.  6.  Finally,  for  hard  spheres  the  angle  of  deflection, 


caused  by  a  binary  collision,  is  independent  of  the  magnitude  of  the  relative 
velocity.  This  enables  us  to  integrate  over  all  values  f^Cl)  (=|W2^(1)|) 

■A 

and  to  evaluate  the  resulting  integrand  for  collision  events  for  which  l^Cl), 
is  a  unit  vector. 

For  the  matrix  elements  an<*  »  determining  the  first  Enskog 

approximations  X^(l)  and  n^(l)  in  (2.4-12),  this  reduction  was  carried  out 

15  •* 

in  a  previous  paper.  Denoting  the  angle  between  -k^  and  the  Z-axis  by 
9^  (see  Fig.  6)  we  obtain 
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^d) 

yTT 


/• i  \  o  H-i  (1)  noc  ^  7r/2 

■all  “  3  TTTy  +  ^  U4  /dW31(1)  /delslne1/dk2/dxcose1 

O  2  O 


•*  -»  -11/2  -*  -» 
W31(l).k2E  XJ/X  Z  K  .  T  (12;13)  K  , 


(1)  35  3^  ■*  11  ^ 

TTOT  =  oo  =  30/dW31(1)/  define  /dk  /dxcose 

O  TT  2  O 


-*  -4  -11/2  St  5 

W,,(l).k,E  '  I  L  :  T  (12 ; 13)  L  . 

31  2  y  V 


For  convenience,  we  have  introduced  the  abbreviations 


(2.7-1) 


E  ■  4l  +  4  -  V331  • 


(2.7-2) 


•Am*  2  2  -fc  2  2 

K  =  W21(W21  '  +  W3l(W3l  "  3®)  * 

t  =  «21w21  +  w31w31  -  i(ff21w31  +  w31?21) 


(2.7-3) 

(2.7-4) 


The  operators  T  (12;13)  commute  with  E  which  represents  the  kinetic  energy. 

In  (2.7-1)  t  is  the  scaled  dimensionless  time  ^W^d)  (see  )•  Just 

as  in  (2.6-8)  the  effect  of  the  T^(12;13)  operators  can  be  written  explicitly 

as 


Z  K  .  T  (12;13)  K  =  K(0) 

•  (K(jR) 

-  K(2)J  + 

(Rl) 

M  U  ^ 

Jk 

-  K(l) 

.  (K(3r) 

-  K(2) }  + 

(R2) 

+  K(0) 

.  (K(jc) 

-  K(2)}  + 

(Cl) 

-  K(l) 

•  (K(3C) 

-  K(2)}  + 

(C2) 

o 

v_/ 

1 

.  (K(3r) 

-  K(l) }  + 

(HI) 

•ft 

+  K(l) 

.  {K(3r) 

-  K(l)}  , 

(H2) 

(2.7-5) 
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ft  ft 

,ft 

4 

:  T^(12:13)  L  =  L(0) 

:  (L(Jr) 

-  L(2)}  + 

(Rl) 

* 

A 

a 

-  L(l) 

:  (L(3r) 

-  L(2)  }  + 

(R2) 

ft 

ft 

ft 

+  L(0) 

:  (L(jc) 

-  L(2)}  + 

(Cl) 

ft 

4 

-  L(l) 

:  (L(3C) 

-  L(2)>  + 

(C2) 

ft 

jm 

ft 

-  L(0) 

:  {L(3r) 

-  L(l)}  + 

(HI) 

ft 

ft 

ft 

+  L(l) 

:  (L(3h) 

-  L(l)  }  , 

(H2) 

where  again  each  term  has  to  be  integrated  subject  to  the  conditions  of 

aft  -ft 

Table  IV.  The  argument  n  in  K(n)  and  L(n)  indicates  for  which  velocities 

oft 

WzfCn),  W^^(n)  these  quantities  should  be  evaluated.  The  kinetic  energy 
E  is  invariant  during  the  collision  process  and  thus  independent  of  n. 

aft 

We  specify  the  perihelion  vector  by  its  polar  angle  0 2  with 
the  Z-axis  (see  Fig.  7)  and  by  its  azimuthal  angle  4>  with  the  XZ  plane  as 
initial  plane.  Hence 


c‘2x  =  sin62cosiJ),  k2y  =  sin02sin<{>  ,  k2z  =  cos02  , 


(2.7-7) 


and 


dk2  =  sin62d02dd>  . 


(2.7-8) 


The  velocity  W^(l)  is  specified  by  its  magnitude  w,  its  polar  angle 

i 

6 ^  with  the  Z-axis  and  the  aximuthal  angle  iJj  in  the  XY  plane,  measured  with 
respect  to  the  plane  through  k2  and  the  Z-axis  as  initial  plane.  Thus 

^31  (l^x  =  wsin93  cos^  +  <l»)  > 

W31^y  =  wsine3  +  ^)  t  (2.7-9) 

W,,  (1)  =  wcos0o  , 

31  z  3  * 
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and 


-*  2 

dWsi(l)  =  w  sinQgdwdG^diJ)  . 


(2.7-10) 


As  a  consequence 

«•  ** 

W31^^*^2  =  w(cos92cos®3  +  >  (2.7-11) 

and  condition  (2.6-10)  becomes 

cos@22  =  cos82cos63  +  sin02sin03cos^  >  0  .  (2.7-12) 

In  terms  of  these  variables  the  triple  collision  integrals  (2.7-1)  read 


A1(l) 


2  nl^ 


385  3 


7  7T  /  2  TT  TT  2TT  2lT  »  “ 


FTTY  =  T  TT7TT  +  4j-,  /do  /de  /d9  /d^/dw/dx 

A0(1)  3  n0(1)  7T2  214  O  10  0  30  O  0  -T 


3  -13/2  «* 

w  sin01cos01sin62sin02COS022E  I  K  .  T  (12;  13)  K  , 

y  ■  y 

(2.7-13) 

n^(l)  22  2^  ff/2  ^  ir  2tt  2tt  °°  °° 

=  ~  ~0  ^r^Qg-^/dili/dw/dx 

o  tt2  o  o  o  ooo  — t 


3  —11/2 

w  sin0^cos0^sin02sin02cos022  E  'EL 


Ty (12*13)  L 


with 


2 

E  =  1  +  w  -  wcos0 2  .  (2.7-14) 

Of  course,  in  addition  the  Integration  variables  are  still  subject  to  the 
conditions  of  Table  IV. 

It  does  not  seem  feasible  to  reduce  the  number  of  integrations 
further  analytically  without  making  approximations.  Therefore,  to  carry 
out  the  remaining  seven  integrations  we  take  resort  to  a  numerical  procedure. 
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The  integration  over  the  velocity  w  extends  to  infinity.  In 

addition,  the  time  x  between  collisions  can  also  become  infinite,  when 

two  of  the  molecules  are  moving  in  parallel  directions.  In  order  to 

apply  a  numerical  integration  procedure,  we  transform  the  integral  into 

an  integral  over  a  finite  integration  region.  For  large  values  of  w 

-4 

the  integrand  decreases  as  w  ;  this  suggests  the  transformation 


u  = 


(2.7-15) 


1  +  wJ 

where  0<u<l.  To  reduce  the  integration  over  t  to  a  finite  interval  we 


proceed  as  follows.  For  a  given  k^,  k ^  and  W^(l)  (0^,  8 ^ ,  0^,  <t>,  and  w 
in  (2.7-13))  we  solve  the  inequalities  (p  -  4  )  and  (p  -  5  )  in  Table  IV  for 
the  variable  x  (together  with  the  requirement  x>  -  x  ) : 

TR,Mi„  '  T  ‘  TR,Ma>  -  “•  R2)  > 


TC,Min  *  T  <  TC,Max^  "  01  •  C2J  ’ 


(2.7-16) 


TH,Min  <  T  <  TH,Max(w  ”  H1,  H2)  * 


These  inequalities  determine  the  range  of  x  values  for  which  a  third  collision 
can  occur  when  k^,  k^  and  W^(l)  are  specified.  It  turns  out  that  T^  ^n“0 
independent  of  the  other  variables.  As  a  new  Integration  variable  we  take  x 
such  that 


x 


■  T 


v,Min 


+  ?  AT 


(v  -  R,  C,  H)  , 


(2.7-17) 


where 


AT 

v 


^v,Max  ^v,Min 


(2.7-18) 
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Hence 


^(D 

X^O) 


,  n,  (1)  oor  06  v/2  IT  7T  2t  2tT  1  1 

I  ^  34  /de1/d02/d03/d«j./d^/du/di 

'  it  2;  o  o  o  oooo 


sinO  ^cos0  ^sin0  ,._sln8  ,^wcos6  23 

’  2  13/2 

u  E 


Z  AT^  K  .  T^(12;13)  K  , 


(2.7-19) 


n.  (1) 

^Tj- 


35 

2 

IT 


,5  tt/2  it 


2tt  2tt  1  1 


^jQ  /d0  /d02/d03/d*/dTp/du/d7 
2'l  o  o  o  o  o  o  o 


sin0^cos0^sin02sin03wcos023 

‘  2  11/2 

u  E 


E  AT 


3 

L 


T  (12 ; 13)  L 
P 


where  AT^  =  ATR2  =  ATR,  etc. 

As  mentioned  earlier  the  time  intervals  AT^  become  infinite  for  certain 
points  in  the  integration  region.  For  recollisions  this  happens  when  02-K), 
where 


AT  *  1  *  1  ; 

R  tfn(2)  x  ff21(l)  sin02 

for  hypothetical  cyclic  collisions  when  03-*-O,  where 


AT  *  1  a  1  ; 

H  tf31(l)  x  tf21(l)  sine3 
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and  for  cyclic  collisions  when  ? ,  the  angle  between  W^^)  and  ^^(l) 
approaches  zero 


*TC  * 


*  1 
sin? 


The  behavior  of  the  integrand  near  these  singularities  can  be  characterized 


aR  ae  aH 

sinSr.slnS.  (—7— g  +  .  -  +  -  -  )  . 

2  3  sin02  sin?  sin6^ 

The  coordinate  system  introduced  was  designed  to  remove  any  singular  behavior 
of  the  integrand  for  62  =  0  and  0^  =  0.  However,  although  the  integral  is 
convergent  (see  Section  3.1),  the  integrand  is  singular  near  ?=0  and  there¬ 


fore  not  yet  suitable  for  numerical  integration.  Since 


lim  sin?  =  lim  lim  {  (y  -  07)^cos^0,  +  \I>Z sin^S,}1^  f 
?-K)  e2 -»u/2  4+0 

the  singularity  can  be  removed  by  the  transformation 


1 2  .2.  ,1/2 


pcosx  =  it  -  2&2  t 
psinx  =  i|>  . 


(2.7-20) 


p  and  x  are  the  polar  coordinates  associated  with  the  variables  n  -  26 2  and 


ir  2tt  it  pMax 

f d02  /dip  =  /dx  fp dp  t 

00  00 


(2.7-21) 


with 


p  =  Min  {—7—  ,r~—  }  • 
Max  sinx  (cosxi 


(2.7-22) 


Thus  the  limiting  behavior  of  the  integrand'  around  ?  =  0  or  p  =  0  becomes 
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d0  2dip 


dxdp 


2.  2„  ,1/2  ' 


sine;  ,1  2  2. 

l-^cos  xc°s  9^  +  ^sin  XSin  0 ^ j 

which  is  well  behaved  for  any  value  of  x ■  Introducing  p  such  that 

c  '  *  pMa*  * 

the  final  form  of  the  triple  collision  integrals  to  be  considered  for 
numerical  evaluation  becomes 


(2.7-23) 


X1(l) 

raj 


2  nl(1) 

3  no(l) 


38S  36  ff/2  T  2r  r  1  1  1 

— 2~  ~ /de^/de^/dfjj/dx/dff/du/dx 

tt  2  o  o  ooooo 


If  p„  sin6. cos0- sin0-sin6_  wcos0„o 
Max  1  1  2  3  23 


u 


2  „13/2 


at  k 

y 


T  (12;13)  K 


(2.7-24) 

n.  (1)  -n  -5  tt/2  tt  2ir  b  1  1  1 

■jpTjj  -  •—  ■—  /dei/d63/d«/dX/diy/du/dT 

lv  TT  2  O  0  OOOOO 


p  p„  sin8.cos0.1sin0-sin0,1  wcos8~_ 
Max  112  3 _ 23 

2  -.11/2 

u  E 


3 

I  AT  L 

y  11 


(12; 13)  L  . 


The  seven  dimensional  integration  region  is  specified  by 


0  <  0]L  <  tt/2 
0  <  e^  <  IT 
0  <  d>  <  2tt 

0  <  X  <  it 

0  <  p  <  1 

0  <  u  <  1 

0  <  T  <  1 


(2.7-25) 
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Condition  (y-2)  in  Table  IV  is  taken  care  of  by  0  <9  <tt/2.  Conditions  (y— 4) 

and  (y-5)  are  satisfied  by  0<t<1  (together  with  x>0  for  R  and  t>-tq  for  C 
and  H) .  The  remaining  conditions  will  be  imposed  explicitly.  When  one  of 
these  conditions  is  violated  the  corresponding  term  in  the  integrand  will 
be  taken  to  be  zero.  We  remark  that  the  triple  collision  transport  integrals 
(2.6-8)  together  with  the  conditions  of  Table  IV  represent  a  complete  solution 
of  the  Choh-Uhlenbeck  integral  equation  for  X^  and  in  which  all  dynamical 
effects  are  taken  Into  account.  Since  we  did  not  have  any  reliable  a  priori 
estimate  of  the  effects  of  the  various  triple  collision  events,  we  retained 
the  complete  collision  integrals  (2.7-24)  together  with  all  conditions  of 
Table  IV  in  the  numerical  integration.  From  the  results  obtained  we  can 
suggest  a  posteriori  approximations  that  would  not  change  the  values  significantly 
as  will  be  discussed  later. 

A  computer  program  was  devised  that  yields  the  value  of  the  integrand 
for  any  point  of  the  integration  region  (2.7-25).  The  following  checks  were 
included  In  the  computer  program: 

1.  For  each  triple  collision  event  we  verified  conservation  of  kinetic 
energy . 

2.  For  each  triple  collision  event  we  verified  conservation  of  angular 
momentum. 

3.  For  each  triple  collision  event  subject  to  the  conditions  of 
Table  IV  we  verified  that  all  other  (redundant)  conditions  mentioned  in 
Section  2.6  were  satisfied  as  well. 

To  obtain  an  estimate  for  the  integrals  we  need  to  select  a  sampling 

21 

procedure.  At  first  we  tried  a  method  of  Haselgrove  for  numerical  integration. 
However,  we  did  not  obtain  a  satisfactory  convergence  for  the  estimates.  As  an 
alternative  we  decided  upon  a  modified  Monte  Carlo  quadrature  suggested  by 
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Haber. 


22+) 


Consider  the  integral 


1  1 

I  f(x)  dx  t 

o  o 


(2.7-26) 


over  the  n  dimensional  vector  x.  The  n  dimensional  integration  region  is 
to  be  pictured  as  divided  into  N  =  kn  subcubes.  In  the  i-th  subcube  two 
independent  random  points  are  chosen  x^  and  y^.  In  addition  we  consider  the 
points  x*  and  y^  which  are,  respectively,  the  points  symmetrically  opposite 
to  x^  and  y^  in  the  subcube.  A  first  estimate  of  the  integral  I  is 


N 


J1 OO  =  2^  '  z  {f(x±)  +  f(y±)}  , 


with  an  estimated  standard  deviation 

N 
E 

i-1 

A  second  better  estimate  is  obtained  by 


1/2 


c^i  =  2^  I  E  tf  (Xj)  “  fC^))2]  . 


J  <k)  =  E  { 
■  i=l 


.  N  f(x  )  +  f(x’)  f  (y“. )  +  f(y!) 

1  -  r  1  - s_l -  +  - - —  }  • 


(2.7-27) 


(2.7-28) 


(2.7-29) 


with  an  estimated  standard  deviation 


°2  2N 


4 

i=l* 


f  (f±)  +  f(y;hSi/2 


(2.7-30) 


+)  The  author  is  indebted  to  Dr.  S.  Haber  for  his  advice  and  to 

Mr.  F.  P.  Karriker  for  his  assistance  in  carrying  out  the  computations. 
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To  test  the  reliability  of  this  procedure  for  our  triple  collision 

transport  integrals,  we  first  evaluate  the  seven  dimensional  integrals 

2 

associated  with  diagram  H2  with  the  additional  conditions  r21(3)  <  1  and 
r^(3)  <  1  at  the  time  of  the  third  collision.  As  pointed  out  in  Section  2.5 
under  these  conditions  we  should  regain  the  values  predicted  by  Enskog.  Thus 
for  these  integrals  we  can  compare  the  estimates  with  the  exact  answer.  The 
results  obtained  for  k  =  2,  3,  4,  5  and  6  are  presented  in  Table  V  together 
with  the  corresponding  standard  deviations.  It  turns  out  that  the  estimates 
k  =  2  and  k  =  3  are  not  yet  reliable  which  is  not  surprising.  However,  the 
estimates  for  k  >  4  are  all  within  each  others  standard  deviations  and  do 
agree  with  the  exact  result. 

Subsequently,  we  applied  this  procedure  to  the  complete  triple 
collision  integrals  (2.7-24).  In  practice,  we  subtract  from  H2  the  con¬ 
tributions  contained  in  the  Enskog  theory  and  compute  directly  the  difference 
between  the  complete  triple  collision  integral  and  the  value  predicted  by 
Enskog.  Since  we  expect  possible  cancellations  between  the  contributions 
of  different  diagrams  (confirmed  by  the  results)  we  combine  for  each  phase 
point  all  six  terms  of  the  integrand  to  obtain  a  maximum  precision  for  the 
total  integral.  From  the  test  calculation  we  conclude  that  we  should  dis¬ 
regard  the  estimates  for  k  =  2  and  k  =  3.  In  the  preliminary  computations 
reported  here,  we  obtained  four  successive  estimates  J^(4),  1^(4),  J^(5),  J2(5), 
indicated  respectively  by  i  =  1 ,  2 ,  3 ,  4  in  Tables  VI  and  VII.  The  results 
obtained  for  the  total  triple  collision  integrals  are  presented  in  Table  VI. 

In  addition,  we  made  an  estimate  of  the  contributions  of  the  six  individual 
triple  collision  events;  these  results  are  reported  in  Table  VII.  We  see 
that  the  four  different  estimates  agree  with  each  other,  as  was  the  case  for 
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the  test  calculation  in  Table  V.  Thus  we  conclude  that  our  preliminary 
calculations  seem  to  yield  a  fair  estimate  for  the  triple  collision 
transport  integrals. 
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Table  V 

Test  of  the  numerical  integration  procedure 


xl(1> 

n1(l) 

X  (1)  3  VaJ 

Enskog 

nQ(l) 

Enskog 

Jx(2) 

-0.32  +  0.29 

-0.53  +  0.35 

J2<2) 

-0.13  +  0.15 

-0.49  +  0.21 

^(3) 

-0.16  +  0.02 

-0.54  +  0.05 

1 

«  1 

^  1 
_ L 

-0.17  +  0.04 

-0.55  +  0.06 

Jx(4) 

-0.203  +  0.017 

-0.634  +  0.025 

J2(4) 

-0.211  +  0.012 

-0.634  +  0.019 

J1(5) 

-0.215  +  0.008 

-0.637  +  0.011 

J2(5) 

-0.213  +  0.007 

-0.634  +  0.009 

Jj_(6) 

-0.213  +  0.007 

-0.629  +  0.009 

m  mmmm  »»  m 

exact 

-0.208 

-0.625 
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Table  VI 


Estimates  for  the  triple  collision  transport  integrals 


i 

AjU)  2  n^l) 

n1(l) 

rilT  3  n“UT 
o  o 

1 

-0.191  +  0.010 

-0.612  +  0.012 

2 

-0.199  +  0.010 

-0.614  +  0.010 

3 

-0.190  +  0.004 

-0.598  +  0.005 

4 

-0.192  +  0.004 

-0.601  +  0.004 
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Table  VII 


Estimates 

for  the  contributions 

of  the  separate 

triple  collision  events . 

4 

XjU)  2  *»1(1) 

nx(l) 

u 

1 

Tjrr  3  TTTIT 

o  o 

n  AD 

O 

R1 

1 

-0.048  +  0.007 

-0.042  +  0.006 

R1 

2 

-0.050  +0.005 

-0.039  +  0.004 

R1 

3 

-0.051  +  0.003 

-0.039  +  0.003 

R1 

4 

-0.049  +  0.002 

-0.040  +  0.002 

R2 

1 

+0.044  +  0.006 

+0.035  +  0.007 

R2 

2 

+0.045  +  0.005 

+0.032  +  0.004 

R2 

3 

+0.047  +  0.003 

+0.033  +  0.003 

R2 

4 

+0.046  +  0.002 

+0.032  +  0.002 

Cl 

1 

-0.023  +  0.003 

-0.036  +  0.004 

Cl 

2 

-0.028  +  0.006 

-0.039  +  0.005 

Cl 

3 

-0.035  +  0.006 

-0.037  +  0.002 

Cl 

4 

-0.032  +  0.004 

-0.037  +  0.002 

C2 

1 

+0.057  +  0.004 

+0.147  +  0.006 

C2 

2 

+0.064  +  0.007 

+0.152  +  0.006 

C2 

3 

+0.069  +  0.006 

+0.151  +  0.003 

C2 

4 

+0.069  +  0.004 

+0.152  +  0.003 
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Table  VII  (continued) 


u 

i 

a1(1)  2  nl(l) 

n1(l) 

3  ^717 

rH 

N— ^ 

o 

c: 

HI 

1 

+0.135  +  0.016 

+0.197  +  0.018 

HI 

2 

+0.141  +  0.012 

+0.206  +  0.013 

HI 

3 

+0.150  +  0.008 

+0.219  +  0.009 

HI 

4 

+0.145  +  0.006 

+0.211  +  0.007 

H2 

1 

-0.355  +  0.017 

-0.913  +  0.021 

H2 

2 

-0.371  +  0.014 

-0.926  +  0.016 

H2 

3 

t0.370  +  0.008 

-0.924  +  0.010 

H2 

4 

-0.370  +  0.006 

-0.920  +  0.008 
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In  this  chapter  we  have  shown  how  the  first  density  coefficients 
of  the  transport  coefficients  are  related  to  triple  collision  transport 
integrals.  In  the  preceding  Section  we  obtained  preliminary  numerical 
estimates  for  these  integrals.  Retaining  the  best  estimates  we  conclude^ 


x.d)  MD 

^=-0.593+0.014,  rm- 

o  o 

where,  for  the  error  estimate,  we  have  replaced 

by  a  95  percent  confidence  interval. 

The  triple  collision  integrals  could  be 

associated  with  the  six  triple  collision  events 

estimates  for  the  contributions  of  the  separate 


0.601  +  0.008  9  (2.8-1) 
the  standard  deviation 

decomposed  into  integrals 
of  Fig.  5.  The  best 
events  are 


U  h™ 

Aod) 

R1  -0.076  +  0.007 

R2  +0.067  +  0.007 

Cl  -0.057  +  0.011 

C2  +0.172  +  0.012 

HI  +0.286  +  0.021 

H2  -0.983  +  0.023 


n1(l)  (2.8-2) 

noCl) 

-0.040  +  0.004 
+0.033  +  0.004 
-0.037  +  0.004 
+0.152  +  0.006 
+0.211  +  0.015 
-0.920  +  0.015  . 


The  theory  of  Enskog  neglects  the  effects  of  all  diagrams  except 


for  a  piece  of  diagram  H2.  To  facilitate  a  comparison  with  the  Enskog 
prediction  -5/8  (2.3-11)  we  rewrite  (2.8-1)  as 


X,  (1)  n,  (1) 

\  =  -0.625J1  -  (5  +  2.2)%];  =  -0.625[1  -  (4  +  1.3)%]  .  (2.8-3) 

o  '  o  v 

+)  In  addition,  there  is  a  collisional  transfer . contribution  +6/5  to  X./X  and 

1  o 

+4/5  to  nn/u  . 

1  o 
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Thus,  as  a  first  conclusion,  we  see  that  despite  the  approximations. 

involved,  the  theory  of  Enskog  presents  a  surprisingly  good  estimate  for 

the  first  density  coefficient,  namely  within  5%.  From  (2.8-2)  it  appears 

that  the  complete  contribution  from  diagram  H2  would  have  overestimated 

the  triple  collision  integrals  by  as  much  as  50%.  Somehow,  this  additional 

contribution  is  cancelled  against  contributions  from  the  remaining  diagrams. 

Secondly,  the  difference  between  the  results  (2.8-3)  obtained  from 

the  Choh-Uhlenbeck  integral  equation  and  those  from  the  Enskog  theory  are 

larger  than  the  error  estimates.  Thus  we  confirm  the  theoretical  expectation 

that  the  theory  of  Enskog  is. not  strictly  correct.  Recently,  Alder  and 

Wainwright  obtained  an  estimate  for  the  first  density  coefficient  of  self 

23 

diffusion  from  molecular  dynamics  calculations.  The  deviations  from 
the  theory  of  Enskog  found  here  for  the  thermal  conductivity  and  viscosity 
have  the  same  sign  as  the  deviation  found  by  Alder  and  Wainwright  for  the 
self  diffusion.  We  have  extended  the  method  of  this  report  to  include  the 
self  diffusion;  preliminary  calculations  carried  out  so  far  are  in  sub¬ 
stantial  agreement  with  the  molecular  dynamics  results. 

As  pointed  out  in  Section  2.6  the  integrals  associated  with  R1  and 
R2  could  not  be  combined  into  one  single  integral  due  to  the  effect  of  a 
possible  fourth  collision  preceding  or  succeeding  the  three  correlated  binary 
collisions.  In  a  quantitative  analysis  Cohen,  Murphy  and  Foch  concluded 
that  in  a  succession  of  three  real  collisions  pertaining  to  our  diagrams  R1 

and  Cl,  the  probability  of  a  fourth  collision  is  only  of  the  order  of  one 

24 

or  two  parts  in  a  thousand.  Our  computer  results  confirm  that  the  pro¬ 
bability  of  a  fourth  collision  is  rather  small  both  for  real  and  hypothetical 
collisions.  Thus  to  a  good  approximation,  certainly  within  the  present 


66 


AEDC-TR-69-68 


precision,  we  can  combine  these  two  integrals  and  forget  the  fourth  collisions. 
Similarly,  we  expect  that  in  first  approximation  we  could  combine  Cl  and  C2, 
resp.  HI  and  H2,  provided  we  exclude  the  overlapping  configurations  which 
should  be  treated  separately.  However,  we  have  not  made  an  explicit  quan¬ 
titative  investigation  of  the  contributions  to  C2  and  H2  when  the  overlapping 
configurations  are  excluded. 

The  results  (2.8-2)  reveal  an  appreciable  difference  between  the  ' 
contributions  Rl,  R2  and  Cl  on  one  hand  and  the  contributions  C2,  HI  and  H2 
on  the  other  hand.  We  note  that  in  C2,  HI  and  H2  a  binary  collision  can 
occur,  while  two  other  particles  are  overlapping  (see  Fig.  5).  In  C2 
1  and  3  are  allowed  to  collide  while  1  and  2  are  still  inside  each  other. 

In  HI  particle  2  and  3  are  allowed  to  collide  while  1  and  3  are  overlapping. 

In  H2  particles  1  and  3  are  allowed  to  collide  while  1  and  2  are  overlapping, 
and  2  and  3  are  allowed  to  collide  while  1  and  2  and/or  1  and  3  are  over¬ 
lapping.  Collisions  with  such  overlapping  configurations  are  not  included 
in  Rl,  R2  and  Cl.  Our  computations  indicate  that  the  probability  for  triple 
collision  events  with  overlapping  configurations  is  rather  large  and  they 
determine  the  major  contributions  to  the  triple  collision  integrals.  On  the 
other  hand,  the  contributions  from  successive  correlated  binary  collisions 
that  are  spatially  separated  are  very  small.  In  (2.8-2)  this  is  illustrated 
by  the  results  for  Rl  and  R2  which  do  not  include  overlapping  configurations 
and  whose  contributions  almost  cancel  each  other. 

These  results  suggest  a  posteriori  that  the  Enskog  values  could  be 
considered  as  a  valid  first  approximation  to  the  triple  collision  integrals; 
this  first  approximation  is  determined  by  those  events  in  H2,  where  in  the 
third  collision  2  and  3  are  colliding,  while  both  2  and  3  are  overlapping 
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with  particle  1.  As  a  second  approximation  those  triple  collision  events 
suggest  themselves  which  contain  collisions  between  two  particles  while  one 
of  them  is  overlapping  with  the  third  particle.  This  second  approximation 
contains  contributions  from  C2,  HI  and  H2.  As  a  third  approximation  one 
could  then  consider  the  effect  of  three  successive  binary  collisions  for 
which  all  pairs  are  spatially  separated.  The  various  effects  of  possible 
fourth  collisions  could  be  treated  as  an  even  higher  order  correction.  A 
theoretical  formulation  of  such  an  expansion  procedure  is  presently  being 
proposed  by  Ernst  and  Van  Beijeren+^ . 

We  should  remark  that  we  have  evaluated  X^  and  only  in  the  first 
Enskog  approximation,  that  is,  we  have  taken  only  the  first  Sonine  polynomial 
contribution.  As  discussed  in  Section  2.2,  the  first  Enskog  approximation  is 
a  good  one  in  the  case  of  the  dilute  gas  values  Aq  and  no »  obtained  from  the 
linearized  Boltzmann  equation.  In  Section  3.5  we  shall  make  some  study  of 
the  effect  of  the  higher  order  Enskog  approximation  on  the  solution  of  the 
Choh-Uhlenbeck  equation  for  the  two  dimensional  model.  The  preliminary 
results  lead  us  to  expect  that  this  expansion  is  suitable  to  obtain  the 
first  density  coefficients  and  Since  the  Enskog  values  -5/8  are 

independent  of  the  Enskog  approximation,  we  only  need  the  expansion  into 
Sonine  polynomials  for  the  term  determining  the  deviation  from  the  Enskog 
theory."  We  expect  that  the  second  Enskog  approximation  would  not  change  the 
values  (2.8-1)  by  more  than  a  few  tenths  of  a  percent.  To  justify  this  suggestion 
we  hope  to  obtain  better  estimates  for  the  triple  collision  integrals,  including 
the  effect  of  the  second  Enskog  approximation,  in  the  future. 


+)  The  author  is  indebted  to  Dr.  M.  H.  Ernst  for  a  valuable  discussion  on 
this  point. 
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CHAPTER  III 

LOGARITHMIC  DENSITY  DEPENDENCE  OF  THE  TRANSPORT  COEFFICIENTS 
FOR  A  GAS  OF  HARD  DISKS 

3.1  INTRODUCTION 


Up  to  this  point  the  treatment  of  the  density  dependence  of  the 
transport  properties  was  based  oh  the  theory  of  Bogoliubov.  In  this  procedure 
the  transport  properties  are  expanded  in  a  power  series  in  the  density  p 


A  = 


X  + 
o 


V 


+  ~2P2  + 


(3.1-1) 


Kf  2  +  .  .  .  . 

n  =  n  +  n,p  +  n0p 

Ol  i. 

For  hard  spheres  we  prefer  to  use  bp  as  a  dimensionless  density.  In  this 

expansion  the  coefficients  A^  and  are,  at  least  in  part,  determined  by 

transport  integrals  related  to  successive  correlated  binary  collisions  between 

the  molecules.  The  coefficients  A^  and  were  related  to  collision  events 

between  three  molecules.  As  mentioned  in  Section  2.7  the  time  between  two 

successive  collisions  among  three  particles  can  become  very  large.  If  one 

analyzes  the  triple  collision  events  of  Fig.  5  for  large  values  of  the  time  x 

between  collisions,  it  turns  out  that  the  probability  of  such  an  event  decreases 
_2 

as  t  .  The  reason  is  that,  in  order  that  for  large  values  of  t  a  third  collision 

A 

will  occur,  the  perihelion  vector  k2  of  the  second  collision  must  be  restricted 
to  a  solid  angle  proportional  to  a2/t2.  If  we  then  integrate  over  all  values 
of  t  up  to  a  maximum  value  T,  the  collision  integrals  converge  as  T  \  when  T 
goes  to  infinity: 


'V 


n,  =  constant  +  lim  — 
T-*»  1 


(3.1-2) 


where  the  constant  term  represents  the  finite  value  obtained  for  the  triple 
collision  integrals. 
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In  an  actual  gas  triple  collisions  with  very  large  times  between 
the  separate  binary  collisions  will  not  occur,  since  for  times  larger  than 
the  mean  free  time  7  a  collision  with  another  molecule  would  interfere  and 
destroy  the  triple  collision  correlation.  Such  an  interference  with  another 
molecule  is  not  accounted  for  in  this  expansion,  since  a  power  series  expansion 
in  the  density  amounts  to  an  expansion  in  terms  of  the  properties  of  groups  of 
isolated  molecules  of  2,  3  .  .  .  etc.  particles.  One  would  expect  that  the 
other  molecules  would  exert  a  so  called  collisional  damping  on  the  triple 
collision  correlation,  so  that  the  upper  limit  T  in  (3.1-2)  should  be  identified 
with  the  mean  free  time  t+^ . 

i 

Q 

X^,  =  constant  +  —  •  (3.1-3) 

7 

Since  the  mean  free  time  7  is  in  first  approximation  inversely  proportional 
to  the  density  p,  the  difference  between  (3.1-2)  and  (3.1-3)  does  not  affect 
the  first  density  coefficients  A^  and  n^,  but  is  of  higher  order  in  the  density 
and  would  affect  the  next  coefficients  and  ^  in  (3.1-1). 

The  second  density  coefficients  and  ^  contain  contributions  from 

correlated  collisions  between  four  particles.  It  was  discovered  by  Dorfman 

25  26  27 

and  Cohen  ,  by  Weinstock  and  Goldman  and  Frieman  that  the  probability  of 

-2  -1 

such  events  no  longer  would  decrease  as  t  ,  but  as  t  .A  similar  situation 
is  encountered  in  the  triple  collision  integrals  for  a  two  dimensional  gas, 
since  in  two  dimensions  for  large  values  of  x  the  perihelion  vector  is 
restricted  to  a  planar  angle  proportional  to  a/x.  As  a  consequence,  after 


+)  Strictly  speaking  the  collisional  damping  modifies  the  probability  of  a 

— T  /  T 

multiple  collision  by  a  factor  e  which  leads  to  the  same  result. 
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integration  over  all  values  of  t,  these  integrals  contain  a  contribution  that 
would  diverge  as  Jin  T.  For  the  reasons  mentioned  above  it  is  to  be  expected 
that  in  reality  the  upper  limit  T  should  be  identified  with  the  mean  free 

time  T  which  is  proportional  to  p-1 .  Thus  it  was  suggested  in  a  previous 

paper  that  X£  and  should  be  replaced  by  X^Jlnp  +  ^  and  n^np  +  80 

that 

2  2 

X  =  X  +  Xnp  +  X'p  Jlnp  +  X  p  +  .  .  .  i 

Ol4  Z. 

(3.1-4) 

2  2 

n  =  n  +  tup  +  nip  &np  +  n0p  +  .  .  .  • 

O  ±  Z  L 

for  a  three  dimensional  gas,  And 

X  *»  X  +  X’pJlnp  +  X.p  +  .  .  .  > 

o  l  i 

(3.1-5) 

n  *=  nQ  +  n^pJinp  +n^p  +  .  .  .  t 

for  a  two  dimensional  gas.  This  conjecture  was  confirmed  by  a  partial 

26  29  30 

resummation  of  the  formal  density  series.  ’  ’ 

An  evaluation  of  the  coefficients  X^  and  of  the  logarithmic  con¬ 
tribution  requires  an  analysis  of  four  body  collisions  which  is  rather  cumbersome. 
For  a  two  dimensional  gas,  however,  the  same  physical  effect  can  be  studied 
from  an  analysis  of  the  triple  collision  integrals  developed  in  the  previous 
Chapter.  For  this  reason  we  have  made  an  evaluation  of  the  logarithmic  density 
dependence  of  the  transport  coefficients  for  a  two  dimensional  model  gas 
consisting  of  hard  disks. 
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3.2  SOLUTION  OF  THE 'LINEARIZED  BOLTZMANN  EQUATION  FOR  HARD  DISKS 


To  obtain  the  transport  coefficients  for  this  model  gas  we  consider 
again  a  perturbation  to  the  velocity  distribution  linear  in  the  gradients 


f  -  f  (1  +  *)  % 

o 

where  f  is  the  (two  dimensional)  local  Maxwell  distribution 
o 


f  =n 


2ffkT  exp 


m(v  -  u)  . 
r  2kT  1 


(3.2-1) 


(3.2-2) 


and  $(V)  has  again  the  form  (2.1-13) 


$(V)  =  -A(V)  V-  -B(V)  V°V:g»  u  .  (3.2-3) 

3r 

The  local  density  n(r;t)  and  average  velocity  u(r;t)  are  again  related  to 
the  distribution  function  f  by  (2.1-1)  and  (2.1-2).  Since  in  two  dimensions 
the  molecules  have  only  two  degrees  of  freedom,  the  average  kinetic  energy 
per  molecule  is  kT.  Thus  the  expression  (2.1-3)  for  the  local  temperature 
T(r;t)  has  to  be  replaced  with 


n(r;t)  T(r;t)  =  /  j  mV2f (r;v;t)dv  .  (3.2-4) 

Substitution  of  (3.2-3)  into  the  expressions  (2.1-4)  and  (2.1-5)  for  the 
fluxes  yields  for  the  two  dimensional  transport  coefficients: 

2 

A  -  |  /dV  fo(V)  (2^  -  2)  V2A(V)  ,  (3.2-5) 


n  =  r  /dV  f  (V)  V°V:V°^  B (V)  .  (3.2-6) 

4  O 

Just  as  in  (2.1-18),  we  have  added  in  (3.2-5)  a  term  which  is  zero  according 
to  the  auxiliary  condition  (2.1-17). 
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Before  studying  the  logarithmic  density  contribution  (Xj^  and  in 
(3.1-5)),  we  first  should  obtain  the  dilute  gas  values  Xq  and  ri0.  As  in 
three  dimensions  we  limit  our  discussion  to  the  spatially  homogeneous  case. 

For  a  two  dimensional  gas  of  hard  disks  with  diameter  o  the  Boltzmann  equation 
reads 


3f  (v^jt) 
3t 


=  o/dv^/dkCvi^.lc)  T(12)f  (v1;t)f  (v2;t)  . 


(3.2-7) 


Using  the  same  Chapman-Enskog  procedure  as  in  three  dimensions  we  obtain  for 
the  linearized  perturbation  (3.2-3) 


.  ,  v  mV. 

l2(VlA(0)(V;L))  =  fo(Vl)  (^  -  2)Vj 


(3.2-8) 


I  <*V°V  B^°^(V  ))  =  f  (V  )  —  V°V 

V  lV  V V  kT  V1V1  » 


to  be  compared  with  (2.2-9).  Here  I2  is  the  two  dimensional  linearized  binary 
collision  operator 

I2(<f(V1))=.-o/dV2/dk(v21.k)fo(V1)fo(V2)T(12)[$(V1)  +  4>(V2)]  .  (3.2-9) 


To  solve  the  linearized  Boltzmann  equation  (3.2-8)  we  approximate 
again  the  solution  by  a  finite  sum  of  the  appropriate  Sonine  polynomials 

1/2  N 


k=l 

a  (® )  /y  \  _  1_  (  m  .1/f2  m  N„1  (o)f  .  (k)f  2. 

B  (V1>  “  no  W  2kT  .  Z  bk  (N)S2  (V  ' 

k=o 


(3.2-10) 
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so  that  from  (3.2-5)  and  (3.2-6) 


,  _  ,  _  2  ,k\1/2  (o),m 

A„  =  A  (N)  *  -  (— — )  a.  (N)  , 
o  o  o  mir  1 


n  =  n  (N)  (5^)1/2b^o)(N)  . 

o  O  iCJ  TV  o 


(3.2-11) 


The  coefficients  a('°^  (N)  and  (N)  are  again  determined  by  the  equations  (2.2-14) 


2  a£o)(N)  a<°}  -  (t  -  1 . N)  , 

k=l  *  *'* 


N-i  ,  V  ,  V 

E  b^°;(N)  b^'  =  6ffn  (£  =  0,  .  .  .  ,N-l)  , 


(3.2-32) 


k=0 


where  the  matrix  elements  a^°^  and  b^°^  are  now  two  dimensional  binary  collision 


integrals 


■m  =4-  .  vvfM>> , 

n  a 

1/2 

b^°^  =  — i—  (— )  /dW  tj°W  *  I  (W°W  S(£)(W2)) 

bk£  2  W  J  “w1wiwis>2  •  VW2  'V'  * 

n  a 


(3.2-13) 


If  we  represent  the  Sonine  polynomials  by  the  generating  function  (2.2-12)  and 
carry  out  the  Integrations  we  obtain 
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_  _  (o)  k  Z  1,. 

Z  Z  a'/e  t  *  ,(1  - 

k=l  £=1  L 


[3(1  -  st)"5/2-2(l  -  st)_i-(l  -  st)“3/2] 


1  s+t  "1/2 

[a 


st)”3/2-st(l 


St)  ^(l 


St)"1/2]  , 


00  oo 

. E  E  bk£  8  C 
k“o  Z=o 


(3.2-14) 


1(1  _  s+t)  [15(1  _  st)”7/2-8(l  -  st)“3-3(l  -  st)-5/2  +  2(1  -  st)~2  +  2(1  -  st)"3/2] 


_  in  s+t.~1/2 [3(1  -  st)"5/2-2(l  -  st)"2-(l  -  st)"3/2] 

8U  2  } 

^  "3/2  "3/2  1  -1/2 
-  -^(1  -  -^)  [(1  -  St)  -st(l  -  St)  Ml  -  st)  yu]  . 


The  resulting  values  for  the  matrix  elements  are  listed  In  Table  VIII  up  to 
the  third  Enskog  approximation,  and  the  resulting  corrections  a£°^ (N)  and 
bk°^  ^  IX*  It  turns  out  that  the  expansion  procedure  is  again  very 

satisfactory.  The  convergence  is  almost  as  rapid  as  in  three  dimensions,  as 
is  shown  in  Table  X. 
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Table  VIII 


The  matrix  elements  a 


(o) 
k  Z 


and 


defined  by  (3.2-14)  for  hard  disks. 


39 

16 


.  (o)  »  51 
bll  16 


a(o)  =  -  ^ 
“23  128 


(o)  4433 

a33  =  1024 


1 

32 


(o)  123 

12  =  "  128 


,  (o)  7569 

°22  "  1024 
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Table  IX 


The  successive  approximations  (N)  and  b^0^ (N)  determining  the  solution 

of  the  linearized  Boltzmann  equation  for  hard  disks. 


First  Enskog  approximation:  N  =  1 


a<o)(l)  =  1 


b(o)(l) 

o 


=  1 


Second  Enskog  approximation:  N  =  2 


a{0><2>  -  1+38 


b(o)(2) 

o 


1  +  50 


a(°) (2)  =  — 
a2  K  *  19 


.(o).-.  _  2_ 
bl  (2)  "  25 


Third  Enskog  approximation:  N  =>  3 
a(°)m  -  i  i  1  i  6g2 

a!  w  “  38  152  x  9997 


*(°)m  -  1  +i_  + _ 29_ 

bo  (3)  "  1  +  50  +  200  x 


a<o)<3) 


2  23  x  65 

19  19  x  9997 


b{0)(3) 


2_  31  x  29 

25  75  x  2521 


a<o) (3) 


4x65 

9997 


b<o)(3) 


4  x  29 
3  x  2521 


2521 
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Table  X 


Convergence  of  the  Sonlne  polynomial  expansion  In  determining  the  dilute 
gas  values  of  the  transport  coefficients 


hard  spheres 

hard 

disks 

N 

JZJ 

O 

r< 

VN) 

\,«> 

xo(D 

"o(1) 

xo(D 

no(D 

1 

1.000 

1.000 

1.000 

1.000 

2 

1.023 

1.015 

1.026 

1.020 

3 

1.025 

1.016 

1.029 

1.022 
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3.3  ASYMPTOTIC  SOLUTION  OF  THE  CHOH-UHLENBECK  INTEGRAL  EQUATION 


To  obtain  the  first  density  correction  to  the  transport  coefficients 
we  start  with  a  formal  density  expansion 


A  =  X  +  T,bp  , 

o  1  ' 

r\  =  n  +  n.bp  , 

o  1 


(3.3-1) 


where  bp  is  the  co-area 

bp  =  ■1'irno^  ■  (3.3-2) 

This  density  expansion  is  a  formal  one  and  the  coefficients  and  will 
turn  out  to  depend  logarithmically  on  the  density.  Pursuing  first  the  formal 
procedure,  we  consider  the  related  expansion  for  the  functions  A(V)  and  B(V) 
in  (3.2-5)  and  (3.2-6)  . 


A(V)  =  A (o)(V)  +  bp  A^(V)  ,  (3.3-3) 

B(V)  ■  B(o)(V)  +  bp  B(1)(V)  . 

(V)  and  *B ^  (V)  are  again  determined  by  the  Choh-Uhlenbeck  Integral 
equation  (2.4-9) 


i2(V1A(1)(vi))  -  -I3(V1A(o)(V1))  f 


I2(vJviB(1)(vi))  =  -If3(vJv^B(°)(V1))  , 


(3.3-4) 


where  1^  is  the  linearized  ternary  collision  integral  operator  in  two  dimensions 

3 

Vt(7l))  ■  -  ^/dVVd?2df3812T<12:3)£o<Vl)fo<V2)fo<V3)  *  <3-3-5> 

m=l 
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Substituting  (3.2-8)  into  (3.2-5)  and  (3.2-6)  and  using  the  symmetry  of  the 
operator  we  obtain  similarly  to  (2.4-11) 

X1  "  "  T  /d^l  ^A^  (Vx)  .T3(i?1A(:o)  (V1))  , 


(3.3-6) 


n,  =  -  f-fd^  ^'^(vp^^B^^)) . 


Using  the  solution  (3.2-10)  of  the  linearized  Boltzmann  equation  we  obtain  in 
the  Nth  Enskog  approximation 


X.  (N)  N  N  .  .  ,  .  ... 

TJX T  =  -  z  E  ak°)(N)  a£°)(l0  “S’  A 

oUJ  k=l  £=1  *  16 


n  (N)  H-l  N-l  ,  .  ~ 

TTTiT"-  £  *  bk  w  a  w  b£  ' 

o  k**o  l=o 


(3.3-7) 


where  and  are  now  two  dimensional  triple  collision  integrals 

*U)  =  JL  /dwt  w  s^^wh.T-^s^V?))  , 

ak£  2  ''mtr  . . 

n  a 

b^^  =  — i—  (— )  fdW  W°W  S(k)(W2)-7  (W°W  S^(W2)) 

°k£  2  W  1  1  Ti  i  '•1;3''112  * 

n  a 

They  can  be  reduced  again  to  the  surface  Integral  form  (2.5-4) 

=  /dW1dW2dW3exp[-(W2  +  W2  +  W2)] 

TT  * 

r 

^  2  •*  ^  «&  (o"\  2 

.  I  W  s;  ;(w  ).0*(12;13)  E  W  S.'  '(M  )  , 

-ml  nr  ,  n  1  n  » 


(3.3-8) 


= 


=  -  ;dW1dW2dW3exp[-(W2  +  W2  +  W2)] 

3  3 

.  L  W°W  S^k)(W2):0*(12;13)  I  W°W  (W2)  , 
mm2  m  «  «  "  «  r 


n  n  2 


(3.3-9) 
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with  0* (12; 13)  defined  by  (2.5-5) 

0*(12;13)  =  Z  /diVd^df  W2l,kl  V3iM'*2  Ty(12;13)  •  (3.3-10) 

vi 

To  obtain  the  coefficient  of  the  logarithmic  contribution,  we  evaluate 

the  Integrand  in  the  triple  collision  Integrals  for  large  values  of . the  time 

t*.  For  this  purpose  the  integrand  is  expanded  in  terms  of  1/t*  and  only  the 

leading  terms  of  first  order  in  1/t*  are  retained.  The  triple  collision  integrals 

are  again  a  sum  of  integrals  over  the  six  triple  collision  events  specified  by 

the  conditions  of  Table  IV,  just  as  in  (2.6-8).  However,  for  the  leading  terms 

in  1/t*  we  can  neglect  the  difference  in  integration  region  between  the  diagrams 

of  the  first  and  second  kind.  Part  of  this  difference  was  related  to  the  effect 

of  overlapping  configurations  for  small  values  of  t*  which  clearly  do  not 

affect  the  asymptotic  behavior  for  large  values  of  t*.  The  remaining  differences 

were  related  to  possible  fourth  collisions  interfering  with  the  three  successive 

binary  collisions;  these  effects  lead  to  contributions  of  higher  order  In  1/t* 

31 

as  was  shown  In  a  previous  publication.  Thus  for  an  asymptotic  analysis  it  is 
sufficient  to  consider 


^  ZdT*/dW^ (l)dif2 (l)dW3 (l)exp [-wj (1)  +  W2(l)  +  W2(l))] 

IT  * 


*  3  3  nOOftA 


di^dkjiw  a>.k1|w31(i).k2  I  nivi  <V-Vi <V]R+  (3.3-11) 

m=l  n=l 


[W  S(k)(W2),W  S(£)(W2)] 
ml  m’nl  n'JC 


with  a  similar  expression  for  in  terms  of  W°W  ^ (W2) .  For  convenience, 

we  have  Introduced  the  abbreviated  notation 
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[$,  f]R  o  [*(0)  -  *(1)].[»(3R)  -  V(2)]  , 

I+.  *]<.  -  td>C0)  -  ♦(1)].[?(3C)  -  ¥(2)]  ,  (3.3-12) 

[♦.  -  [<K0)  -  ♦(1)].[Y(3h)  -  ¥(1)]  , 

where  ®(i)  =  $(W(i))  and  i/>(i)  =  i|;(w(i))  are  vector  or  tensor  functions  of 

jk 

the  velocities  W(i)  and  the  product  Is  a  scalar  product.  The  three  terms 
in  (3.3-11)  .should  be  integrated  over  those  phase  points  leading,  respectively, 
to  a  recollision  R,  a  cyclic  collision  C  and  a  hypothetical  cyclic  collision  H 

t 

(see  Fig.  8).  These  Integration  regions  are  specified  by  the  conditions  (y-2) 
through  (y-5)  of  Table  IV. 


Figure  8.  The  three  triple  collision  events  for  the  collision  integrals 
determining  the  coefficient  of  the  term  logarithmic  in  the  density 
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The  leading  term  in  1/x*  can  be  evaluated  analytically.  To  illustrate 

#*(1) 

the  procedure  we  consider  the  recollision  contribution  a^  (R)  in  some  detail. 

Just  as  in  three  dimensions  we  first  integrate  over  the  total  velocity 

W  =  i(w,  +  w“„  +  VO  and  over  all  directions  and  the  magnitude  of  W„.(l).  To 
O  3  1  2  j 

2i 

evaluate  the  remaining  integral  we  take  W2^(l).  in  the  direction  of  the  Z-axis 
and  specify  the  perihelion  vectors  k-  and  k2  by  t*ie  angles  and  ©2  as 
indicated  in  Fig.  7.  The  velocity  W31(l)  (strictly  speaking  the  scaled  ve¬ 
locity  w‘31(l)/W21(l)  15)  is  specified  by  its  length  w  and  its  polar  angle  63< 
Thus , 


*g(l)  /R\ 
ak£  (R) 


.  gr  ~  +v/2  +ir/2  +tr/2  “ 

=  -  ^7 f -  /dx/dW  /dW_. (1)W,. (1)  /  d6  /  d0.  I  d0  / 
T*71  0  21  21  -ti/2  ^tv/2  2r/2  3  1 


dw 


COS0. 


(3.3-13) 


3  ^ 


•  exp[-(W2(l)  +  W2(l)  +  W2(l))  ]w2cos0^cos  (02-03)  E  E  [W^S^  ^  (W^)  ,  ^wn^]R 

m=l  n=l 


Just  as  in  (2.7-1)  x  is  the  scaled  dimensionless  time  x*W21(l).  It  is  con- 
venient  to  consider  the  perihelion  vector  of  the  third  collision  as  a  new 
integration  variable.  We  specify  k^  by  an  angle  0^  as  indicated  in  Fig.  9. 

After  the  first  collision  (between  1  and  2)  particle  1  is  at  the  origin  of 
the  coordinate  system  and  particle  2  is  moving  in  the  Z-direction  with  the 
unit  velocity,  as  shown  in  Fig.  7.  After  the  second  collision  (between  1  and  3) 

•A 

particle  1  will  move  in  the  direction  of  k2>  A  singularity  will  occur  when 

because  for  small  values  of  02  particle  1  will  move  approximately  in  the 
Z-direction  and  be  able  to  collide  again  with  2  even  for  large  values  of  x. 

At  the  time  of  the  second  collision  2  will  have  travelled  a  distance  x  in  the 
Z-direction.  After  the  collision  between  1  and  3  the  velocity  of  particle  1 
is  equal  to  wcos03  in  the  limit  considered.  Therefore,  the  time  x'  required 
for  1  to  overhaul  2  is  approximately  x/(wcos03~l) ,  so  that  the  total  time 
x  +  x'  between  the  first  and  the  last  collision  becomes 
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X  +  X* 


WCOS&2 
wcose3“l  * 


(3.3-14) 


The  positions  of  the  particles  at  the  time  of  the  third  collision  are 
Indicated  In  Fig.  9.  Particle  2  is  at  B  and  has  travelled  a  distance  t  +  t' 

In  the  Z-directlon  when  colliding  with  1  again.  Neglecting  terms  of  order 
0 2  and  constant  terms  when  compared  with  x,  we  have  AC462  x  OCcr0„  (t  +  x ' ) . 

The  distance  BC  Is  essentially  k3x  -  sin0  ,t.  Since  AB  *=  sin0, ,  we  obtain 

WCOS0- 

sin0,  =  -  sin9,-  + - .  ■■■*  x0«  »  (3.3-15) 

4  1  wcos0--l  2 


so  that 


1  wcosOyl 
x  wcos03 


cos0 , d0 ,  • 
4  4 


(3.3-16) 


In  (3.3-13)  we  replace  the  integration  over  by  an  integration  over  0^. 

Since  the  Jacobian  (3.3-16)  Is  proportional  to  1/x,  it  Is  sufficient  to 
evaluate  the  remaining  part  of  the  integrand  for  0^  =  0  (terms  of  order  @2 
are  equivalent  to  terms  of  order  1/x) .  We  can  indicate  this  formally  by  intro¬ 
ducing  a  6-function  6(0^)  and  integrating  again  over  02*  Thus 


(R) 


1  sfl  Ar  -a  9  +ir/2  +*/ 2  +v/2  -hr/2 

—tv  “  /dw  /dW,. (1)W,. (1)  /  dS.cosO.  f  d0.cos0.  /  d0,  /  d0,  I  dw 

*  ^  -*ll  1  -v/2  *  — tt/2  — tt/2  J  1 

cosO^ 


?  2  2  cos(02-e3) 

exp[-W^(l)+W‘(l)+W3(l))]w(wcos03-l)  cosQ-; — 


(3.3-17) 


3  3 

•  E  E  [W  S(k)(W2),  W  S(2,)(W2)]  6(0  )  +  /dT0(l/x2)  . 

--in  m  n  i  n  k.  * 

m=l  n=l 
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The  first  term  in  this  expansion  is  proportional  to  1/t.  As  mentioned 
in  Section  3.1  we  take  as  the  upper  limit  for  the  time  integral  a  mean  free 
time  7  inversely  proportional  to  the  density  bp.  Thus  writing 


(3.3-18) 


we  obtain  for  the  recollision  contribution  (R)  to  the  coefficient  of  the 
logarithm: 


“u  <R) 


^  ^  2  +ir/2  +ir/2  +t/2  +tt/2  00 

/dw  /dC,,Q-*M1^de,cos0-  /  de.cose,  /  de„  /  de,  /dw 

°  21  Z1-tt/21  l-n/2  4  Sw/2  ±ir/2  3  1 

COS6g 


.  exp [- (W^ (1)+W2 (1)+W^ (1) ) ]  w(wcos03-l) 


cos(e2-e3) 

cos03 


.  I  l'  [w  sjk)  (Wz) ,  w  s|l)  (WZ)  ]  6  (02)  . 

-  -  ml  m  ni  n  k  z 

m=l  n=l 


(3.3-19) 


The  contributions  from  the  cyclic  and  hypothetical  cyclic  collisions 

can  be  evaluated  in  a  similar  manner.  We  first  introduce  again  the  perihelion 

vecotr  k3  of  the  third  collision  as  a  new  variable.  The  asymptotic  recollisions 

occurred  when  02^O.  The  asymptotic  cyclic  collisions  will  occur  when  -ir/2 

31 

and  0  ■>  +ir/2;  both  limits  give  an  equal  contribution.  The  asymptotic 
hypothetical  cyclic  collisions  occur  when  93-K).  The  velocity  functions  to  be 
evaluated  before  and  after  the  third  collision  are  symmetric  in  all  three  particles. 
Therefore,  this  difference  before  and  after  the  third  collision  is  the  same 
function  of  k3  regardless  of  the  fact  between  which  two  particles  the  third 
collision  actually  occurs.  As  a  result  the  contributions  from  cyclic  and 
hypothetical  cyclic  collisions  can  be  simply  related  to  the  expression  found 


86 


AEDC-TR-  69-68 


for  recollisions.  After  introducing  as  a  new  variable,  in  the  case  of 
recollisions  the  remaining  integrand  had  to  be  evaluated  for  02=0.  For  cyclic 
and  hypothetical  cyclic  collisions  we  have  to  consider  the  same  integrand  for 
9 2  "  “1T/2,  and  0^=0,  respectively.  Thus  the  total  contribution  reduces  to 


x(1) 

k£ 


1  tX  ,1,  ^  2  +7T/2  +7t/2  +n/2  +it/2  » 

~N -  /dW  /dW„.  (1)W-  (1)  /  d0.cos0  /  d0.cos0.  /  d0„  f  d0_  f  dw 
4»  ir  o  LL  LI  _/0  1  l_i«  4  4  2  3  , 


-rr/2  x  in/2 


n/2  in/2  J  1 
COS0. 


2  0  0  |.COs(0_-0  ) 

expt-(W^(l)  +  W2(l)  +  W3(l))]  w (wcos  0 3~1 )  — —  -  ™ - 


(3.3-20) 


3  3 


m=l  n=i 


The  three  delta  functions  pick  out  the  contributions  from  recollisions,  cyclic 
collisions  and  hypothetical  cyclic  collisions,  respectively."*"^  For  the  matrix 
elements  in  (3. 3-^18)  we  obtain  similarly 


/1\  x  AJX  m*.  Jk  2  +7r/2  +TT/2  +n/2  +n/2  “> 

1  =  / dW  /dW91  (1)W,.  (1)  I  d0n cos0.  f  d0.cos0.  f  d0.  /  d0_  /  dw 

kt  H  o  21  21  1  1  4  4  ,,  2  ,,  3  . 

n  -n/2  -n/2  -n/2  -n/2  1 

cos03 

222  |cos(0  -0  )| 

.  exp[-(W'  (1)  +  wj(l)  +  W,(l)>]  w(wcos0  -1) - -  -  (3.3-21) 

X  L  J  J  COS0- 


'  !,  ^  [6<02)  +  S(92  -  ¥  -  5<63>>. 


+)  When  integrating  the  delta  function  6(02  -  ^-)  the  lower  limit  of  should 

7T 

be  taken  as  -(y  +  e) . 
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The  logarithmic  density  dependence  (-3.3-18)  of  the  matrix  elements 
and  implies  =  X^Jlnbp  +  X^  and  =  n^Jlribp  +  so  that  the 

expansion  (3.3-1)  for  the  transport  coefficients  should  be  replaced  with 


X  -  Xq  +  XjbpJlnbp  +  X^bp  +  .  .  .  , 

(3.3-22) 

n  ■  n  +  n.'bpJinbp  +  n,bp  +  .  .  .  . 

o  1  1 

The  coefficients  X^  and  ri|  are  determined  by  the  new  asymptotic  triple 
collision  integrals  and  : 


Xj(N) 

irm 


k=l  £=1  X 


*■[00 

VI) 


N-l  N-l 
E  E 
k=o  l=o 


(N)b^o)(N)e^  . 


(3.3-23) 


The  coefficients  X^  and  can  be  considered  as  the  sum  of"  three 
contributions  associated  with  the  three  triple  collision  events  of  Fig.  8. 

We  note,  however,  that  in  (3.3-20)  and  (3.3-21)  the  matrix  elements  and 

are  reduced  to  triple  collision  integrals  containing  the  dynamics  of 
recollisions  only.  In  these  integrals  the  velocities  (0)  and  (2)  are 
related  to  the  integration  variables  by  (2.6-25)  and  (2.6-26).  The  final 
velocities  (3^)  are  related  to  the  variable  k^  (or  6^)  by 

tf2i (3r)  =  W21(2)  -  2W21(2).k3k3  , 


<V 


W31(2) 


Wn(2).k3k3 


(3.3-24) 
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3.4  LOGARITHMIC  DENSITY  DEPENDENCE  OF  THE  TRANSPORT  COEFFICIENTS 
IN  THE  FIRST  ENSKOG  APPROXIMATION 


The  logarithmic  density  dependence  of  the  transport  coefficients 

in  the  first  Enskog  approximation  is  determined  by  the  matrix  elements 

* 

j 

o  21 


and  0^.  If  we  carry  out  the  integration  over  W  and  W„.  (1)  we  obtain 


11 


oo 


similarly  to  (2.7-1) 


31 


„<»  -  |  +  fifi cose, '  Tde.'Tdo,  /  «. 

11  3  oo  1024  ,2.^2  1  l.l/2  4  4.,^  2_>/2  3  i 


+ir/2 


+tt/2  +ir/2  00 


cos  6, 


(3.4-1) 


w(wcos0,-l)  |cos(0_-0j I  . 

cos§V  ■  tK-  Ky5(82>  +  s<e2- 1 '  S(83>]  ’ 


E 


11/2 


}(D 

oo 


ds+j/2  +tt/2  +tt/2  +tt/2  » 

^77-^  /  depose.  /  d0.cos0.  /  d0o  /  dO*  / 
64  2  /0  1  1  /0  4  4  2  3  . 

tt  — tt/2  -tt/2  -tt/2  -tt/Z  i 


COS0, 


(3.4-2) 


w(wcos0--l).  |cos(0_-6,)|  ±  ±  _ 

~57F - ‘L-  LVS(V  +  s(e2-  ¥  -  6(63>1  . 


where  E,  K  and  L  are  again  defined  by  (2.7-2),  (2.7-3)  and  (2.7-4).  In  the 
coordinate  system  under  consideration  E  is  given  by  (2.7-14): 


E  =  w  +  1  -  wcos0 


3  ’ 


(3.4-3) 


The  components  of  the  various  velocities  (n)  to  be  substituted  in  K(n)  and 
L(n)  are  listed  in  Table  XI. 


89 


AEDC-TR-69-68 


Table  XI 

i 


The  velocities  to  be  substituted  In  the  triple  collision  Integrals  a 


and  0. 


a) 

k£  ' 


(1) 

kl 


w21(o) 

w31(°) 

“21<2> 

“31 <2> 


<V 


■  X-component 
-2sin0^cos0^ 
wsln0 3-sin0^cos0^ 

0 

wsinO^ 

0 

wslne^ 

2(l-wcos03)sin0^cos0^ 

wsin03+(l-wcos03)sin0^cos0^ 


Z-component 

2 

l-2cos  0^ 

2 

WCOS03“COS  0^ 

1 

WCOS03 

l-WCOS©3 

-WCOS03 

(l-wcos03) (l-2cos^0^) 

2 

-wco803~(  1-wcoge^cos  0^ 
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Omitting  terms  odd  in  the  integration  variables  61  and  0^  we  obtain 


[KrK]R  -  {K(0)  -  K(1)}.{K(3r)  -  K(2) } 


(3.4-4) 


4(l-wcos02)^[4w^sin^02+  (l+wcosO^)  (l-2wcos@3)  ]sin^0jcos^0jsin^0^cos  0^  » 


i,  L]r  -  {L(0)  -  L(1)}:{L(3R)  -f(2)> 


2  2  2  2  2 
18(l-wcos03)  sin  0^cos  0^sin  0^cos  0^  . 


(3.4-5) 


Introducing  these  results  into  (3.4-1)  and  carrying  out  the  integrations  over 


0.  and  0.,  we  obtain 
1  4 


n\  9  m  ign  -Ht/2  -  w(wcos0^-l)  |cos(02“03) 

„  (1)  ±,  a(^)  4.  f  dfl  f  d&  f  dw  .  —  -4- -  - — 

Q-l  -I  O  +  9  Oft  •*  atJ9  *  a09  J  QW  H/2  ftftofl 


*11  3  woo 


2  80  ““2  '  .,3  , 

7T  -n/2  -it/2  1 


cos8. 


COS0, 


(3.4-6) 


[4w^sin^03-2w^cos^03  -  (wcosO^-l)]  [5(02')  +  6(e2“  “  6(63)]  p 


-hr/2  +tt/2  “  w(wcos0_-l)' 


ir  — tt/2  — tt/2  1 

cos0_ 


,9/2 


|cos(02-93) 

COS0„ 


(3.4-7) 


.[fi(e2)  +  fi(e2-  |>  -  6(e)]  . 

31 

The  remaining  integrals  can  be  readily  evaluated.  Thus  we  obtain  in  the 
first  Enskog  approximation 
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r£y  -  '  “S*  -  -  1555  (8276  '  1920  Jr '  4185  *n3)  '  l3-4'8> 


/ 1  \  -1  7T 

raT  -  -  soo  -  -  255  *7  <82°  -  224  vf  -  405  *n3)  •  (3'4-9) 

The  contributions  of  the  three  triple  collision  events  individually  are 
presented  in  Table  XII. 


92 


AEDC-TR-69-68 


Table  XII 


Contributions  to  the  coefficient  of  the  logarithm  in  the  first  Enskog 
approximation. 


diagram  (Fig.  8) 


> '(1)/X  (1) 
1  o 


48  (  5  tk 

25  2  K  9  l/?' 


1  (2764  -  2511  £n3) 


1080  2 
TT 


H 


16  f3 
135  2 

TT 


R  +  C  +  H 


1850  %  <8276-1920,1  -  6185  *-3) 

ir  f  3 


diagram 


n,’(l)/n  (1) 
1  o 


48  i3_  n  14  tt 
25  2  U  "  27 


£ 


~  (452  -  405  £n3) 


H 


64 

225  2 

7T 


R  +  C  +  H 


"  2§5  ^2  (82°  -  224  5  -  405  £n3) 

it  /3 
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3.5  LOGARITHMIC  DENSITY  DEPENDENCE  OF  THE  VISCOSITY 
IN  THE  SECOND  ENSKOG  APPROXIMATION^ 


From  the  results  obtained  In  the  preceding  Section  we  conclude  that 
the  transport  properties  do  depend  logarithmically  on  the  density  with  non¬ 
vanishing  coefficients  X*  and  n^.  Nevertheless,  one  may  wonder  whether 
the  first  Enskog  approximation  provides  a  reasonable  estimate  for  these 
coefficients.  From  the  fact  that  the  Sonine  polynomial  expansion  (3.2-10) 
of  the  solutions  A^°^  (V^)  and  B^(V^)  of  the  linearized  Boltzmann  equation 
is  a  convergent  one.  we  may  expect  that  the  resulting  expansion  for  the 
density  correction  (3.3-6)  to  the  transport  coefficients  converges  as  well, 
provided,  of  course,  that  is  replaced  with  a  finite  operator  1^  by  taking 
the  collislonal  damping  into  account.  However,  in  order  to  conclude  that  a 
reliable  estimate  has  been  obtained  for  X|  and  n|,  it  is  not  sufficient  to 
know  that  the  expansion  procedure  converges  ,  but  we  need  to  verify  whether 
the  next  correction  terms  are  small,  as  was  the  case  for  the  dilute  gas 

values  X  and  n  (see  Table  X).  For  this  purpose  we  have  extended  the 
o  o 

analysis  to  determine  the  corrections  in  the  second  Enskog  approximation.  In 
this  report  we  discuss  the  results  obtained  for  the  viscosity. 

The  second  Enskog  approximation  n^(2)  is  given  by 


hi  (2) 


^  -  -  CCb^o) (2)}20^>  +  b^O>(2)b^°)(2){0^)  +  0^}  +  (b{O>(2))20^>]  .  (3.5-1) 


x<1} 

xki 


jd) 

JkA 


It  can  be  easily  demonstrated  that  the  matrix  elements  a,;*'  and  0^„'  are 
symmetric 


|d>  -  ftd> 


ki 


(3.5-2) 


+)  The  research  of  this  Section  was  carried  out  in  cooperation  with 
Dr.  H.  T.  Wood  at  the  National  Bureau  of  Standards,  Washington,  D.  C. 
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The  reason  is  that  in  Fig.  8  the  same  triple  collision  events  are  obtained 
whether  the  diagrams  are  read  from  top  to  bottom  or  from  bottom  to  top 
(except  for  a  permutation  of  the  particle  numbers) . 

For  the  matrix  element  given  by  (3.3-21),  we  need  to  consider 


l  l  fw°tf  s^-W),  «s<°V)]  =  i  i  [w°w  (3  -  wf),  W°W ■  ]R 
«  .  mm2  m  n  n  2  nR  .  ,  m  m  m  nnn 

m=l  n=l  m=l  n=l 


3  3 

-±ns* 


3  3 


E  E  3[«  ,  «  ]_-  «].  . 

_  -  m  m  n  n  R  m  m  m  hhk 

m*5!  n=*l 


Thus 

- 38 


a) 

oo 


r*  .  --ht/z  -Hit 

IT  *  -TT/2  -IT 


+tt/2  -Hit/2  +tt/2  -Hr/ 2  00 


/  dS.cose.  /  d0_  /  d0„  /  dw 
/2  4  — tt/2  ?„/2  3_1^_ 

COS0- 


cos(02~03) 


exp  [-  (W2  (1)+W2  (1)+W3  (1)  ]  w(wcos03-l)  - 


(3.5-3) 


3  3 

.  E  E  [W°W Vj,  W°W  ]  [6(0,)  +  6{02-  j)  -  <5  (03)  ]  . 

-  -  m  m  m  n  n  k  z  n  ^  J 

m=l  n=l 


As  in  the  preceding  Sections  we  express  the  velocity  functions  in  terms  of 
the  total  velocity  Wq  =  (W^  +  W2  +  ftj)  and  the  relative  velocities  *21  and 


*31  * 


3  ^  oft  a  2  ^ 

EWW  =3WW  +  -s-  L  > 
n=1  n  n  o  o  3  • 

EHH 2  =  -^[9W  W  (3W2+  |fi)  +  (6W2  +  2E)  ?  -  W  » 
,  m  m  m  9  o  o  o3  o 

m=l 


(3.5-4) 


+  24  (ft  L)W  -6W  K  +  4W  .M]  f 

0  O  O  n  1 


(3.5-5) 
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where  we  have  introduced  two  new  tensors 


+) 


«  =  Wn  + Wn  + 


2[W21W21W31  +  W21W31W21  +  W31W21W21  +  W21W31W31  +  W31W21W31  +  W31W31W21'1  9 


N  =  W21W31W31  +  W31W21W21  "  W21’W31(W21W31  +  ^31W21*  * 


(3.5- 


Note  that 


a  a  3 

IsM  =  -  y  K 


In  two  dimensions 


and 


N  =  0  . 


3 

Tr  A.M  =  -  A.K  . 


(3.5- 


(3.5- 


In  (3.5-5)  terms  odd  in  vanish  after  integration  over  Wq.  The  terms  that 

contain  only  W  and  E  are  constants  of  the  motion  and  vanish  also.  The 
o 

2  i 

remaining  terms  in  (3.5-3)  contain  [L,  L]  and  as  a  consequence  they  can  be 

K 

related  to  the  expression  (3.4-2)  for  As  a  result  we  obtain  after 

oo 

integration  over  Wq  and  T^-j^l) 


(1)  _  (1)  _  1  (1) 

eio  "  eoi  -  -  6 


/  B 
4  oo 


(3.5- 


For  the  matrix  element  we  consider 


1  *  (Wm}  »  “Xs!15  (Wn}  ] 

,  .mm2  m  n  n  l  n 

m=l  n=l 


3  3 

-  r^O-  *0**  i  2  ^O , 

1  E  9 [W  W  ,W  W  ]„  -  3  [W  W  W  ,W  W  ]„ 
.  .  m  m  n  n  R  in  m  m  n  n  R 

m=l  n=l 


-3[thr  ,w°w vh  + 

mm  n  n  n  K  m  m  m  nnn  * 


+)  M  is  a  two  dimensional  tensor  of  the  third  rank.  The  scalar  product 


defined  as  T.  Z  Z  M'^^.  Symmetrized  dyadics  are  indicated  by  AB  =  -j(AB 


7) 


8) 

9) 


10) 


'  is 
’+  BA). 
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so  that 


(1)  _  (1)  (1)  .  (1) 
811  ■  36oi  +  3610  9800  + 


rr  ^  ^  7  +ir/2  +7J/2  +tt  f'L  +tt/2  ® 

+  ^4  V  f-/’dWo‘rdW2i(1)W2p)  f  d0!0080  f  de^cose^  I  d02  f  d03  /  dw 

TT  — tt/2  -tt/2  -7t/2  — tt/2  1 


+71 /2 


+tt/2  +tt/2  ® 


COS0. 


2  2  2  |  cos  (0„“9-i  I 

.expf-^d)  +  W2(l)+W3(l))]  w(wcos03-l)  - - 

3  3 

.  Z  Z  [W°W  W l,  W°W  W^]_[5(0J  +  6(0  -  j)  -  6(0  )]  . 

,  ,  mmm  n  n  n  K  z  z  z  j 

m=l  n=l 


(3.5-11) 


To  evaluate  this  integral  we  substitute  (3.5-5)  into  (3.5-11).  We  see  that 

i  3 

upon  integration  over  W  several  terms  can  be  related  again  to  [L,  L]  . 

0  -  3 

Moreover,  as  a  result  of  relations  (3.5-8)  most  terms  involving  K  and  M  can 


be  related  to  [K,  K]  and,  therefore,  to  the  expression  (3.4-1)  for 

K 

-  -i-  6^.  Thus  we  obtain  after  integrating  over  IT  and  W01  (1) 
11  3  oo  91 


m  73  m  rn  315  «/3+7r/2  +1T/Z  "" ,c 

+  a!r  +  4^  ^  f  d6,cos0,  f  d9,cos0,  I  d0„  S  d0,  f  dw 


+tt/2 


21 


+tt/2  +tt/2  “ 


11  48  oo 


' 128  TrV/r1™ v/2  2;/2  3 1 


COS0, 


w(wcos03~l) 

,11/2 


|cos(0  — 03)|  $  dt  „ 

■  Jr-  IM,  tlE[«(92)  +  «<e2-  f  -  «<e3>)  , 


(3.5-12) 


E“'“  3 

where  the  energy  E  is  again  given  by  (3.4-3).  To  evaluate  the  last  term  we 
substitute  the  velocities  of  Table  XI  into 


[M,  M]R  =  {S(0)  -  M(l)}  :  (M(3r)  -%(2)} 

Upon  comparing  with  (3.4-4)  we  find 


+77/2  +r/2  2  5  27+tt/2  +tt/2 

f  dO^os©..  /  d0.cos0.[M,  M]„  =  —  /  d01cos01  f  d0.cos9,  [K,  K]d  (3.5-13) 
/o  i  1  4  4  K  4  <0  1  ±  _  /«  4  4  a. 

—ti/2  -tt/2  — tt/2  -tt/2 
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so  that  this  term  can  also  be  related  to  afP  -  \  8^.  We  obtain  finally 

11  3  oo 


(3.5-14) 


Substitution  of  the  results  for  and  3^  into  (3.5-1)  yields 


for  the  first  density  coefficient 


n*(2) 

Viy 


125ie(1)  +  24a.(P 
_ oo _ 11 

1250 


(3.5-15) 


Thus  the  second  Enskog  approximation  n^(2)  can  be  expressed  in  terms  of  the 
known  results  (3.4-8)  and  (3.4-9)  for  A^(l)/A^l)  and  n^(l)/rio(l) : 


n{(2) 


1251  ni(1)  +  12  Xl(1) 


1250  n  (1)  625  A  (1) 

o  o 


(3.5-16) 


3.6  DISCUSSION  OF  THE  RESULTS 

As  indicated  in  Section  3.1,  we  expect  on  general  grounds  that  the 
successive  correlated  binary  collisions  will  cause  a  logarithmic  density 
dependence  of  the  transport  properties.  For  a  two  dimensional  gas  of  hard 
disks  the  density  dependence  of  the  transport  properties  can  be  represented  by 

A  =  A  +  A'bp£nbp  +  A, bp  +  .  .  .  * 

o  i  i 

(3.6-1) 

n  =  n  +  n,'bp£,nbp  +  n,bp  +  .  .  .  . 

o  i  x 

In  this  chapter  we  developed  a  method  to  determine  the  coefficients  Xj^  and 

Although  the  procedure  was  carried  out  for  hard  disks  in  particular,  the 
method  can  be  used  to  predict  the  logarithmic  contribution  to  the  transport 
properties  for  any  two  dimensional  gas  of  molecules  with  short  range,  spherically 
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symmetric,  forces.  The  coefficients  and  were  expressed  in  terms  of 
(asymptotic)  collision  integrals  related  to  the  recollisions  (R) ,  the  cyclic 
collisions  (C)  and  the  hypothetical  cyclic  collisions  (H)  of  Fig.  8. 

The  numerical  values  obtained  for  and  of  hard  disks  are 
presented  in  Table  XIII.  The  Table  includes  also  the  contributions  from  the 
three  diagrams  R,  C  and  H  individually.  The  values  X £ ( 1)  and  n^(l)  obtained 
in  the  first  Enskog  approximation  are  listed  in  the  second  and  third  columns. 

The  values  obtained  for  the  viscosity  coefficient  n^(2)  in  the  second  Enskog 
approximation  are  presented  in  the  fourth  column.  It  turns  out  that  the 
second  Enskog  approximation  modifies  the  values  found  in  the  first  approxi¬ 
mation  by  about  1%  only.  Thus  we  expect  that  we  have  found  the  correct  order 
of  magnitude  for  the  logarithmic  density  dependence  of  thermal  conductivity 
and  viscosity  for  the  model  considered. 

From  the  last  column  of  Table  XIII  we  see  that  the  second  Enskog 
approximation  affects  the  three  contributions  to  by  the  same  relative 
amount  of ft  1%.  Since  the  total  value  for  n|(l)  is  of  the  same  order  as 
the  contribution  from  the  most  important  diagrams  (R  and  H  in  this  case) , 
approximately  the  same  correction  is  obtained  for  the  total  contribution. 

However,  this  is  not  necessarily  true  in  general.  For  instance,  if  we  evaluate 

the  logarithmic  contribution  to  the  self  diffusion  in  the  first  Enskog  approxi- 

32 

mation,  the  total  coefficient  turns  out  to  be  an  order  smaller  than  the  contri¬ 
butions  of  the  individual  diagrams,  due  to  some  cancellations  between  these  con¬ 
tributions.  Although  in  this  case  the  first  Enskog  approximation  leads  again  to 
a  reasonable  estimate  for  these  individual  contributions,  it  is  no  longer  evident 
that  it  will  lead  to  a  good  estimate  of  the  total  effect.  This  will  be  dis¬ 
cussed  in  a  future  report. 
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From  the  results  of  Table  XIII  we  conclude  that  the  coefficients 

of  the  term  logarithmic  in  the  density  is  non  vanishing.  However,  the 

fact  that  these  coefficients  are  different  from  zero,  can  be  expected  on 

much  more  general  grounds.  Inspection  of  the  collision  integrals  (3.3-20) 

and  (3.3-21)  shows  that  in  any  Enskog  approximation,  the  three  collision 

events  R,  C  and  H  will  lead  to  a  contribution  of  the  form  tt  /3  (a^  -  bRy^) , 

tt/3  (a_  -  b„2.n3)  and  tt  a  ,  respectively.  Thus  the  three  collision 
C  C  ri 

events  lead  to  three  different  irrational  numbers  and  a  cancellation  is  therefore 
unlikely  in  any  higher  order  Enskog  approximation. 

The  results  of  Table  XIII  also  reveal  that  the  sign  of  the  logarithmic 
density  contribution  depends  on  the  particular  form  of  the  collision  integrals 
and  is  not  the  same  for  the  thermal  conductivity  and  viscosity.  Thus  it 
could  be  that  both  the  sign  and  the  magnitude  of  the  logarithmic  density 
dependence  would  be  sensitive  to  the  intermole cular  potential  as  well,  and 
it  does  not  seem  possible  to  predict  the  sign  from  a  priori  considerations. 

As  discussed  in  Section  3.1,  the  same  physical  effects  encountered 
in  the  triple  collisions  of  a  two  dimensional  gas,  will  occur  in  the  quadruple 
collisions  of  a  three  dimensional  gas.  Thus  we  submit  that  our  results  for 
the  two  dimensional  model  gas  provide  significant  evidence  that  also  in  an 
actual  three  dimensional  gas  a  logarithmic  density  dependence  must  be  expected, 
as  expressed  in  equation  (3.1-4). 
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Table  XIII 


Coefficient  of  the  logarithmic  density  dependence  for  the 
transport  coefficients  of  hard  disks. 


diagram 

AjU) 

n[(l) 

n[(2) 

n'(2)  -  nj(l) 

Aod) 

no(D 

R 

+0.00258 

-0.02005 

-0.02002 

-0.2% 

C 

-0.00088 

-0.00551 

-0.00553 

+0.4% 

H 

-0.02080 

+0.04992 

+0.04956 

-0.7% 

R+C+H 

-0.01909 

+0.02436 

+0.02401 

-1.4% 
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CHAPTER  IV  + 

DENSITY  DEPENDENCE  OF  EXPERIMENTAL  TRANSPORT  COEFFICIENTS  ) 

4.1  INTRODUCTION 


In  this  Chapter  we  want  to  investigate  whether  the  theoretically 
predicted  density  dependence  of  the  transport  properties  can  represent 
the  experimental  data.  Up  to  the  present  experimental  transport  properties 
have  usually  been  represented  by  a  power  series  in  the  density  (3.1-1). 


A  = 


A  + 
o 


A1P 


A/  2 

+  a2p  + 


(4.1-1) 


n  =  n  +  n.p  +  n9P  +  .  .  . 

Our  analysis  of  the  two  dimensional  model  gas  demonstrated  that  such  a 
power  series  cannot  be  justified  theoretically  and  we  provided  substantial 
evidence  that  instead  one  should  expect  a  density  dependence  that  contains 
terms  logarithmic  in  the  density 


2  2 

A=A  +  A.p  +  A'p  £np  +  A„p  +  .  .  .  , 
O  1  z  z 

2  2 

n  =  n  +  n^p  +  £no  +  +  .  .  . 


(4.1-2) 


It  is  the  purpose  of  this  Chapter  to  establish  some  criteria  to  test 
whether  meaningful  density  coefficients  A^  and  can  be  deduced  from  the 
experimental  data.  In  particular  we  shall  try  to  determine  the  first 
density  coefficients  A^  and  and  make  an  assessment  of  the  precision 
with  which  they  can  be  obtained  from  the  presently  available  experimental  data. 

In  order  to  determine  the  first  density  coefficients,  one  initially 
looks  at  the  experimental  data  in  the  linear  region.  But  to  decide  which 
data  are  actually  in  the  linear  range,  an  estimate  of  the  higher  order  terms 


+)  This  Chapter  was  prepared  in  collaboration  with  Dr.  H.  J.  M.  Hanley  and 

33 

Dr.  R.  D.  McCarty  of  the  National  Bureau  of  Standards,  Boulder,  Colorado. 
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in  the  density  expansion  is  desirable.  ‘Moreover,  it  is  not  evident 
a  priori  that  the  best  value  for  this  coefficient  is  obtained  from  the 
data  in  the  linear  region  alone.  For  example  to  determine  the  second 
virial  coefficient  B  in  the  equation  of  state  (1.1-1)  it  is  often 

34 

advantageous  to  include  experimental  pVT  data  in  the  quadratic  range. 

For  these  reasons  we  try  to  assess  the  effect  of  higher  order  terms  on  the 
determination  of  the  first  density  coefficient.  We  use  equation  (4.1-2) 
for  this  purpose,  since  it  represents  the  density  expansion  suggested  by 
the  theory.  However,  for  comparison  we  shall  make  a  few  remarks  concerning 
the  results  obtained  when  the  data  are  fitted  to  a  simple  power  series 
(4.1-1). 

In  Section  4.2  we  outline  a  general  procedure  to  analyze  experimental 
transport  properties  as  a  function  of  density.  In  Section  4.3  we  apply 
this  method  to  the  thermal  conductivity  of  neon.  In  Section  4.4  we  present 
first  density  coefficients  and  deduced  from  the  experimental  vis¬ 
cosities  and  thermal  conductivities  of  a  number  of  gases.  In  Section  4.5 
we  summarize  and  discuss  these  results. 
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4.2  CRITERIA  FOR  A  CONSISTENT  REPRESENTATION  OF  THE 
EXPERIMENTAL.  DATA  AS  A  FUNCTION  OF  DENSITY 


We  assume  that  the  theory  predicts  a  density  representation  (not 
necessarily  a  power  series)  with  finite  coefficients 


X(T,p)  =  E  X  (T)f  (p)  , 
k=o 


n(T,p)  =  I  nk(T)fk(P)  , 
k=o 


(4.2-1) 


where  the  density  coefficients  Xk  and  are  a  function  of  the  temperature 
alone.  The  functions  f^(p)  are  ordered  such  that 

lira  fk+1(o)/fk(p)  =  0  .  (4.2-2) 

p-H) 

Theoretically,  (4.2-1)  contains  an  infinite  number  of  terms  representing 
the  density  variation  around  p=0.  In  practice,  one  uses  a  finite  number 
of  terms  to  represent  the  transport  properties  in  a  finite  density  interval 


N-l 

X  (T , p )  =  E  Xk(T)fk(p)  , 
k=o 

N-l 

n(T,p)  =  E  nk(T)fk(p)  . 
k=o 


(4.2-3) 


As  pointed  out  in  Chapter  II,  the  present  theory  suggests  fQ(p)  =  1  (excluding 

2  2 

Knudsen  effects),  f^(p)  =  p,  f£(p)  =  P  £np  and  f^(o)  =  p  .  However,  we  shall 

also  include  in  the  discussion  a  power  series  for  which  ^k(p)  =  P  • 

Recently, investigators  have  discussed  some  of  the  factors  affecting 

the  precision  with  which  the  virial  coefficients  of  the  equation  of  state 

34  35 

can  be  deduced  from  the  experimental  pVT  data.  *  In  trying  to  determine 
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the  experimental  values  for  the  coefficients  Xfc  and  nk  in  (4.2-3)  several 

additional  complications  present  themselves.  The  precision  of  the  second 

virial  coefficient  B  in(l*l- 1 )  can  be  enhanced  by  an  independent  knowledge 

of  the  temperature,  so  that  the  first  virial  coefficient  RT  no  longer 

represents  a  degree  of  freedom  in  analyzing  the  pVT  data.  However,  the 

Chapman-Enskog  values  X ^  and  ri  cannot  be  predicted  with  an  accuracy 

significantly  better  than  the  experimental  accuracy,  presumably  due  to 

insufficient  independent-  information  on  the  intermolecular  potential. 

This  is  despite  the  fact  that  a  detailed  theory  is  available  for  their 

description.  In  practice,  therefore,  the  low  density  values  X  and  n 

'  o  o 

have  to  be  extracted  from  the  experimental  data  themselves. 

Furthermore,  in  order  to  deduce  physically  meaningful  coefficients 
from  experimental  data  with  a  given  accuracy,  one  ideally  needs  to  know 
a  priori  whether  the  chosen  function  is  a  correct  representation  of  the 
data.  This  requires  not  only  a  theoretical  knowledge  of  the  form  of  the 
functions  fk(p),  but  also  an  estimate  of  the  magnitude  of  the  coefficients, 
so  that  one  can  decide  which  terms  can  be  neglected  in  a  given  density  in¬ 
terval.  Again,  this  condition  is  satisfied  for  the  equation  of  state  of 
34 

simple  gases.  For  the  transport  properties,  however,  we  do  not  know 
the  functions  ffc(p)  for  k>3,  nor  do  we  have  reliable  a  priori  estimates  of 
the  magnitude  of  the  coefficients  for  k=l,2,3.  As  a  result  we  shall  be 
forced  to  conclude  that  at  the  present  we  cannot  determine  the  coefficients 
for  k^2  with  enough  confidence  to  make  them  amendable  to  theoretical  inter¬ 
pretation. 

Nevertheless,  in  order  to  proceed,  one  can  lay  down  some  minimum 
criteria  which  have  to  be  satisfied  if  (4.2-3)  is  to  be  compatible  with  the 


105 


AEDC-TR-69-68 


data.  The  first  criterion  is 

(i)  The  standard  deviation  of  the  fitted  data  points.  When  the  standard 
deviation  of  the  points  fitted  with  (4.2-3)  exceeds  the  experimental  error 
of  the  data,  and/or  when  a  curve  representing  the  difference  between 
experimental  and  calculated  points  shows  systematic  deviations,  the  pro¬ 
posed  function  cannot  be  a  correct  representation  of  the  given  experimental 
data.  Although  this  criterion  rules  out  certain  possibilities,  it  is 
rather  weak  in  that  it  still  allows  many  interpolation  formulas  to  represent 
the  data.  In  particular,  we  shall  see  that  in  a  given  density  interval 

the  transport  coefficients  can  always  be  fitted  to  a  polynomial  with  or 
without  a  term  logarithmic  in  the  density. 

To  reduce  the  class  of  possible  correct  functions  we  consider  two 
additional  requirements: 

(ii)  The  coefficients  of  the  fitted  equation  should  be  independent  of  the 
size  of  the  density  interval.  If  (4.2-3)  is  the  correct  equation  to 
represent  the  data  in  a  given  density  interval,  one  may  fit  it  to  the  data 
and  return  values  for  the  coefficients.  When  the  equation  is  fitted  to 
data  in  a  subinterval  the  values  of  the  coefficients  then  obtained  must 

be  equal  to  (within  their  standard  deviation)  to  the  corresponding 
coefficients  found  from  the  fit  for  the  complete  density  interval.  This 
should  be  true  for  any  subinterval  as  long  as  all  coefficients  are  significant. 

(iii)  The  coefficients  should  be  independent  of  the  number  of  terms  in 
the  equation.  Assume  that  equation  (4.2-3)  is  the  correct  equation  to 
represent  the  data  in  a  given  density  interval  and  that  one  determines 
values  for  all  coefficients  (0<_k<N-l).  When  the  density  interval  is  suf¬ 
ficiently  reduced  so  that  the  last  term  (k  =  N-l)  is  no  longer  significant, 

N-2  N-2 

a  refit  in  the  reduced  interval  with  the  equation  £  X,  f^(p)  or  I  ^k^k^^ 

k=o  k=o 
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must  return  the  same  values  (within  their  standard  deviations)  for  the 
first  N-l  coefficients  (0£k<_N-2)  as  are  returned  from  the  fit  with  the 
full  equation  (4.2-3)  in  the  complete  density  interval. 

We  shall  call  any  function  (4.2-3)  satisfying  the  criteria  i,  ii 
and  iii  consistent  with  the  experimental  data.  In  addition  to  these 
criteria,  we  shall  also  require  that 

(iv)  All  coefficients  are  to  be  significant  for  a  proper  fit.  Violation 
of  this  condition  does  not  exclude  the  possibility  that  the  function  under 
consideration  is  theoretically  correct.  However,  the  presence  of  insigni¬ 
ficant  terms  affects  the  values  returned  for  the  significant  coefficients. 

Our  analysis  of  the  transport  data  is  thus  to  investigate  whether 

i 

(4.2-3)  is  consistent  with  the  experimental  data  so  as  to  return  reliable 
estimates  of  the  coefficients  A^  and  n^.  If  a  certain  function  of  the 
density,  linear  in  the  coefficients,  is  a  correct  representation  of  the 
data,  the  best  estimates  of  the  values  of  the  coefficients  can  be  obtained 
by  the  method  of  least  squares.  Of  course,  when  the  data  are  correlated 

3 

and  of  unequal  variance,  these  factors  should  be  properly  taken  into  account. 

In  practice,  we  proceed  as  follows.  Suppose  we  have  several  isotherms 
for  A  and  n.  To  reduce  the  effects  from  internal  inconsistencies,  we  con¬ 
sider  a  set  of  data  from  only  one  experimental  source  at  a  time.  We  have  n 
data  points  A(pi)  or  nfpj)  (i  =  1,  .  .  .n)  at  a  given  temperature  T,  ordered 
such  that  p^+^>p^.  To  ensure  equal  variance  of  the  data  points  we  decide 
whether  or  not  weighting  of  the  data  is  desirable,  basing  our  decision  on  the 
experimental  estimates  of  the  possible  errors  in  A  or  n  and  the  density  p. 

It  turns  out  that  the  effect  of  errors  in  the  independent  variable  p  is 
usually  negligible  compared  to  that  of  the  errors  in  the  experimental 
transport  properties. 
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and  determine  Aq,  A^  or  nQ,  from  a  least  squares  fit  for  the  first  few 
data  points  (is  3  or  4).  We  then  add  the  next  datum  point  and  redetermine 
the  .coefficients  and  note  the  criteria  i-iv.  This  procedure  is  repeated, 
the  criteria  noted  for  each  fit,  until  at  a  certain  number  of  points, 
corresponding  to  densities  up  to  a  density  p^,  equation  (4.2-4)  no  longer 
represents  the  experimental  data  correctly.  Such  a  conclusion  can  be 
drawn  when  the  criteria  are  not  satisfied,  i.e.,  when  the  standard  deviation 
begins  to  increase  and/or  when  the  deviation  curve  begins  to  indicate  a 
systematic  deviation.  Furthermore,  the  values  returned  for  the  coefficients 
begin  to  change  systematically  with  the  size  of  the  density  interval. 

At  this  point,  that  is  when  (4.2-4)  ceases  to  be  valid,  the  next 
density  dependent  term  A2f2(p)  or  n2f2(p)  is  added  to  the  fitting  equation 
and  the  procedure  is  repeated.  Once  again,  the  criteria  are  noted  after 
each  fit  and,  in  addition,  the  significance  test  on  the  new  coefficients 
A 2  or  n2  is  used  to  recheck  that  the  original  estimate  for  the  limiting 
density  for  a  linear  fit  is  valid  (A2  or  n2  should  not  be  significant 
for  PlP^)  ■ 

In  principle,  when  enough  experimental  data  are  available,  this 
procedure  can  be  continued  until  all  terms  in  equation  (4.2-3)  are  significant. 

From  the  results  obtained  by  this  procedure  we  can  obtain  some  insight, 
whether  a  proposed  density  expansion  (4.2-3)  is  consistent  with  the  experimental 


+)  Strictly,  we  should,  and  did  on  occasion,  fit  the  data  first  to  the  equation 
A=A^  or  n=nQ .  However,  usually  there  are  not  enough  data  as  a  function  of  density 
in  the  Chapman-Enskog  regime  to  warrant  this  extra  step. 
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data.  Specifically,  here  we  shall  fit  the  data  to  (a)  the  linear  equation 
(4.2-4),  (b)  to  a  quadratic  polynomial,  (c)  to  a  cubic  polynomial  and 
(d)  to  the  equation 


X  =  X  +  X.p  +  X’p2£np  +  X_p2  » 
o  1  Z  Z 

2  2 

n  =  n  +  n.p  +  nip  u.np  +  n„p  , 

o  l  z  z  ' 


(4.2-5) 


which  contains  the  first  four  terms  of  the  theoretically  predicted  density 
expansion. 


4.3  ANALYSIS  OF  THE  THERMAL  CONDUCTIVITY  OF  NEON 

Ve  discuss  in  some  detail  the  results  of  this  procedure  for  the 

set  of  thermal  conductivity  data  of  neon  obtained  by  Bolk,  Stigter  and 
37 

the  author.  In  the  experiment  care  was  taken  to  measure  the  data  points 
at  regularly  spaced  density  intervals  in  an  appreciable  density  range. 

The  variation  of  the  thermal  conductivity  of  neon  as  a  function  of  density 
is  shown  in  Fig.  10.  The  actual  experimental  values  are  listed  in  Table 


+)  Tables  XIV  through  XVIII  can  be  found  at  the  end  of  this  Section.  In 

this  report  the  thermal  conductivity  is  expressed  in  cal/cm.s.K,  since 

most  authors  have  presented  their  experimental  data  in  this  unit.  For  the 

same  reason  the  viscosity  is  expressed  in  g/cm.s  or  poise.  The  presently 

2 

recommended  SI  values  for  X  and  n  are  W/m.K  and  N.s/m  ,  resp.  The  conversion 

9 

factors  are:  1  cal/cm.s.K  =  418.68  W/m.K;  1  poise  =  0.1  N.s/m  . 
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Figure  10.  Experimental  thermal  conductivities  of  neon  as  a  function  of  density. 
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Because  equation  (4.1-2)  contains  a  term  logarithmic  in  the 
density(  we  prefer  to  express  the  density  in  dimensionless  units.  Therefore, 
we  nave  retained  the  amagat  unit  of  density  used  in  the  original  experi¬ 
mental  tables. ^  The  density  in  amagat  Represents  the  ratio- of  the  actual 
density  to  the  density  p  _  of  the  gas  at  standard  temperature  and  pressure 

bit 

3 

(0°C  and  1  atm.).  For  neon  one  amagat  corresponds  to  0.00089979  g/cm  . 

Three  thermal  conductivity  isotherms  are  available,  namely  at  25.1°C 
S0.1°C  and  75.0°C.  The  experimental  method  suggests  that  any  random  error  is  due 
to  errors  in  the  measurement  of  the  temperature  difference  across  the  gas 
layer  or  of  the  heat  current  through  the  gas  layer,  and  is  therefore  most 

likely  a  relative  error,  that  is,  an  error  proportional  to  X.  This  suggests 

_  36  , 

a  weight  factor  X(P±J  for  the  least  squares  analysis.  Actually,  we 

chose  lX(Pl)/100]  ^  as  the  weighting  factors,  so  that  the  method  of  least 

squares  returns  the  estimated  standard  deviation  of  the  data  points  as 

a  percentage. 

Linear  fit.  The  results  obtained  by  fitting  the  data  to  the 
linear  equation  (4.2-4)  are  presented  in  Table  XV  for  the  three  isotherms. 

The  second  column  gives  the  density  interval  0<p<p  and  the  third 

lUoX 

column  the  number  n  of  data  points  from  Table  XIV  used  In  fitting  the 

equation.  In  columns  four  and  five  the  values  returned  for  the  coefficients 

X  and  X.  are  presented  together  with  their  standard  deviation  cr^  and 
o  1  A0 

a  •  The  last  column  shows  the  estimated  standard  deviation  of  the  data 
1 

points.  It  turns  out  that  up  to  a  density  of  120  amagat  the  coefficients 
within  their  standard  deviation  do  not  depend  on  the  size  of  the  density 
Interval.  The  standard  deviation  a  does  not  increase  with  density.  The 
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DENSITY,  amagat 

Figure  11.  Piot  of  deviations  (Agxp  -  A^)/!^  in  %  for  the 
thermal  conductivity  of  neon  in  the  linear  region. 

deviation  curves  (shown  in  Fig.  11  for  the  largest  density  interval)  confirm 
that  the  linear  equation  fits  the  data  to  within  their  experimental  precision. 
Thus  the  linear  equation  satisfies  the  criteria  i,  ii  and  iv  of  Section  4.2. 
However,  when  data  points  at  densities  beyond  120  amagat  are  included, 

(not  shown  here) ,  the  standard  deviation  begins  to  increase  and  the  coef¬ 
ficients  Xq  and  begin  to  change  systematically.  Thus  we  conclude  that  the 
linear  equation  (4.2-4)  is  consistent  with  the  experimental  data  at  densities 
up  to  120  amagat.  The  most  probable  value  for  the  first  density  coefficient 
*1»  virtually  independent  of  the  temperature,  deduced  by  fitting  the  data  to 
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a  linear  equation  is  0.77x10  7  cal/cm.s .K.am  with  a  standard  deviation 
of  0.02x10  7  in  the  same  units. 

Quadratic  fit.  We  then  first  tried  to  represent  the  data  over  a 
wider  density  range  by  a  polynomial,  as  has  been  the  general  practice  in  the 
literature  up  to  now.  For  this  purpose  we  consider  the  quadratic  equation 


X  =  Xo  +  AlP  +  T2p2  .  (4.3-1) 

The  results  obtained  when  the  data  are  fitted  to  this  equation  are  presented 
in  Table  XVI.  The  quadratic  term  does  not  become  significant"*^  until  p  =  200 
amagat,  in  agreement  with  our  previous  conclusion  that  at  p  =  120  amagat  a 
linear  equation  is  still  sufficient.  We  note  that  in  the  units  chosen 

-3 

is  about  10  times  X^.  Up  to  a  density  p  =  440  amagat  the  quadratic 
equation  satisfies  conditions  i,  ii  and  iv;  that  is,  the  coefficients  are 
independent  of  the  density  interval  and  the  standard  deviation  a  does  not 

A 

increase  with  density.  (At  densities  beyond  440  amagat  we  begin  to  detect 
systematic  deviations.)  However,  the  estimate  for  the  first  density 
coefficient  X^  obtained  by  fitting  the  quadratic  equation,  varies  from 
0.64x10  7  cal/cm.s. K.am  to  0.71x10  7  cal/cm.s. K.am,  in  distinct  disagreement 
with  the  value  (0.77  +  0.02)xl0  7  cal/cm.s  .K.am  deduced  from  the  linear  fit.. 
Thus  the  quadratic  equation  is  the  density  interval  considered  violates 
condition  ill. 

Cubic  fit.  Next  we  consider  the  cubic  equation 


,  *'2  ***  3 

X  =  Xq  +  X^P  +  +  X^P  f 


(4.3-2) 


results  for  which  are  presented  in  Table  XVII.  The  cubic  term  does  not 
become  significant  until  p  ■  560  amagat.  Up  to  p  =  800  amagat  the  coefficients 
do  not  vary  with  the  density  interval  and  the  standard  deviation  o,  is  constant . 


+)  As  a  test  of  significance  for  a  linear  coefficient  Xfc  we  require  that  the 
2  2 

quantity  T  =  )  is  larger  than  the  value  of  the  F  0,c  distribution  for  the 

corresponding  degrees  of  freedom. 


.975 


36 
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(The  fact  that  a.  is  the  same  for  all  equations  reinforces  our  original 

assumption  that  the  random  error  is  approximately  a  relative  error 

independent  of  the  density).  However,  the  cubic  equation  again  does  not 

reproduce  for  the  first  density  coefficient  the  value  of  0.77x10  ^  cal/cm.s .K.am 

suggested  by  the  linear  fit,  nor  does  it  reproduce  for  X^  the  value 
-10  2 

&  0.8x10  cal/cm.s. K.am  estimated  from  the  quadratic  equation  in  Table  XVI. 

From  these  results  we  conclude  that  neither  a  quadratic  equation  in 
a  density  interval  up  to  p  =  400  amagat,  nor  a  cubic  equation  in  a  density 
interval  up  to  p  =  800  amagat,  is  consistent  with  the  experimental  data. 

One  could  relax  conditions  ii  and  iii  by  allowing  for  variations  of 
the  coefficients  within  a  95%  confidence  interval  i  which  is  related  to 
the  standard  deviation  by 

4  =  ot<025  ,  (4.3-3) 

where  t  ^5  *-8  t^ie  2.5  percent  point  of  the  "student  t"  distribution  for 

36+) 

the  corresponding  degrees  of  freedom.  But,  even  so,  the  values 

(0.77  +  0.06)xl0  ^  cal/cm.s .K.am  and  (0.67+0.06)xl0  ^  cal/cm.s. K.am  for 

X.  +  A..  ,  returned  from  the  linear  and  quadratic  fit,  do  not  agree  very 
1  X1 

well,  while  the  values  returned  for  X 2  from  the  quadratic  and  cubic  fit 
are  well  outside  each  other's  95%  confidence  interval. 

A  serious  problem  is,  of  course,  the  possible  effect  of  systematic 
errors  in  the  data.  A  statistical  analysis  can  only  account  for  the  effects  of 
random  errors  and  rat  of  systematic  errors.  The  effect  of  systematic  errors 
could  be  estimated,  e.g.,  by  comparing  sets  of  data  from  a  number  of  dif¬ 
ferent  experimental  sources  for  the  same  gas  at  the  same  temperature.  A 


+)  When  the  degrees  of  freedom  are  large,  t  q2^5?  2. 
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more  decisive  check  could  be  made  if  the  results  were  compared  with  those 

from  experimental  data  with  an  appreciably  higher  precision,  so  that  the 

various  terms  would  become  significant  at  other  (lower)  densities.  The 

effect  of  two  simple  systematic  errors  can  be  easily  assessed.  A 

systematic  error  independent  of  the  density  (e.g.  an  error  in  the  radiation 

constant  )  would  only  modify  the  value  returned  for  Xq  equally  for  all 

equations.  If  there  is  a  systematic  error  proportional  to  X  (e.g.  an 

38 

error  in  the  distance  between  the  hot  and  the  cold  plate  ),  all  coefficients 

will  have  this  same  percentage  error.  Thus,  systematic  errors  of  these  kinds 

would  not  affect  the  conclusions  of  this  analysis. 

Logarithmic  fit.  After  these  negative  results  for  the  polynomial 

equations,  we  fitted  the  data  with  equation  (4.2-5)  that  includes  a  term 

logarithmic  in  the  density  as  predicted  by  the  theory.  The  presence  of 

this  logarithmic  term  leads  to  some  complications  not  encountered  when 

2 

fitting  polynomials.  It  does  not  seem  desirable  to  consider  the  p  £np 
2 

and  the  p  terms  separately.  For  instance,  if  we  change  the  density  units 
such  that  p"  =  p/pq»  the  coefficients  X^  and  X£  of 

X  =  X  +  *X.p  +  Xlp^np  +  X_F2  >  (4.3-4) 

O  X  Z  Z 

are  related  to  X^  and  in  (4.2-5)  by 

x;  =  Xlp  ^  ;  X0  =  Xlp  ^£np  +  X0p^  .  (4.3-5) 

ZZo  ZZooZo 

Therefore,  the  relative  contribution  of  the  two  individual  terms  depends 

on  the  particular  unit  chosen  for  the  density.  In  fact,  if  we  had  chosen 

p  such  that  inp  =  -  X„/X  the  quadratic  term  in  (4.3-4)  would  have 

o  o  Z  Z 
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disappeared  completely.  This  problem  could  be  avoided  by  using  the  equation 

X  =  X  +  X„p  +  X 'p2£n  *  (A. 3-6) 

o  1  2  p 

o 

where  (in  addition  to  Xo,  X^,  X£)  PQ  is  a  constant  to  be  determined. 

However,  equation  (4.3-6)  is  nonlinear  in  the  coefficient  pQ  and  therefore, 

less  suitable  to  least  squares  analysis.  Thus,  in  practice  we  prefer  equation 

2  2 

(4.2-5)  and  test  the  sum  Xp  £np  +  \^P  for  significance. 

The  results  when  the  thermal  conductivity  data  are  fitted  to  this 

2  2 

equation  are  presented  in  Table  XVIII.  The  sum  Xp  £np  +  X^p  becomes 
significant  at  p  =  240  amagat  similar  to  the  quadratic  results  given  in 
Table  XVI.  However,  since  (4.2-5)  still  contains  four  degrees  of  freedom, 
the  individual  terms  become  significant  at  much  higher  densities,  comparable 
to  the  densities  where  all  terms  in  the  cubic  polynomial  became  significant. 
The  fit  is  within  the  experimental  error  up  to  p  =  720  amagat;  at  densities 
beyond  720  amagat  we  detect  systematic  deviations. 


0  100  200  300  400  500  600  700 

DENSITY,  amagat 

Figure  12.  Plot  of  deviations  ^exp-Xcaic)/xcalc  ln  *  for  the  thermal  conductitity 

of  neon  using  equation  (4.2-5) 
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We  note  that  in  tne  density  interval  considered  is  again 
approximately  constant  and  the  deviations  of  the  data  points  are  random, 
as  shown  in  Fig.  12.  The  value  returned  for  the  coefficient  A^  does 
not  vary  with  the  density  interval;  this  is  also  true  for  the  coefficients 
A^  and  A2  when  they  become  significant  individually^.  Thus  the  equation 
satisfies  criteria  i,  ii  and  iv.  The  results  for  the  75°C  isotherm  are 
slightly  out  of  line;  here  the  coefficients  become  significant  at  a  later 
stage.  But  even  here,  when  finally  at  p  =  720  amagat  all  coefficients  are 
significant,  the  coefficients  are  similar  to  those  at  the  other  temperatures. 

Furthermore,  from  the  fit  up  to  p  =  720  amagat  we  note  that,  on 
the  average,  A^  =  0.79  x  10  7  cal/cm.s .K.am  with  a  standard  deviation  of 
0.03x10  7,  in  agreement  with  the  value  (0.77+0.02)xl0-7  cal/cm.s. K.am 
deduced  for  A^  from  the  linear  fit  in  Table  XV.  Thus,  condition  ill  is 
satisfied  as  well.  We  therefore,  conclude  that  equation  (4.2-5)  is  con¬ 
sistent  with  the  experimental  data.  Of  course,  any  physical  conclusions 
from  this  result  should  be  considered  with  the  same  caution  as  mentioned 
earlier  when  discussing  the  results  for  the  polynomial  fits. 

To  obtain  a  realistic  error  estimate  for  A^,  we  prefer  to  replace 

the  standard  deviation  a.  with  a  95%  confidence  interval  A^  .  From  the 

1  A1 

agreement  between  Tables  XV  and  XVIII  we  conclude  that  these  tables  yield 


+)  Upon  a  density  transformation  p=ppQ,  the  standard  deviation  of  the 

coefficient  in  (4.3-5)  transforms  as 

2  4  2  2  2  2 

aT  =  p^ai'£n  p„  +  2ai'i  enp„  +  » 

X2  0  X2  0  X2^2  °  X2 

2 

where  cr  ,  is  the  related  off-diagonal  element  of  the  variance-covariance 
matrix  of  the  coefficients. 
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a  meaningful  value  for  the  first  density  coefficient  A^.  On  the  other 

hand,  it  seems  that  a  quadratic  polynomial  would  have  underestimated  this 

coefficient  A^  by  as  much  as  10  percent. 

We  cannot  draw  a  similar  conclusion  about  the  significance  of 

the  values  returned  for  A^  and  A2*  To  check  the  reliability  of  these 

coefficients  by  our  procedure,  we  would  have  to  add  the  next  term  (or 

2 

combination  of  terms)  beyond  the  p  term  to  the  equation  and  compare  the 
results  for  the  new  equation  with  those  of  Table  XVIII.  The  character  of 
these  terms  have  not  been  investigated  theoretically  in  sufficient  detail 
as  yet,  nor  are  there  enough  experimental  data  to  investigate  their 
effect.  However,  we  think  it  reasonable  to  assume  that  such  higher  order 
terms  will  no  longer  affect  the  value  returned  for  the  coefficient  A^ 
of  the  linear  term. 
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Table  XIV 


37 

Experimental  Thermal  Conductivity  data  for  neon 
25.1°C  50. 1°C  75,0°C 


P 

AxlO4 

P 

AxlO4 

P 

AxlO4 

n 

amagat 

cal/cm.s ,K 

amagat 

cal/cm.s .K 

amagat 

cal/cm.s.K 

1 

0.92 

1.166 

0.86 

1.236 

0.83 

1/296 

2 

10.01 

1.174 

10.06 

1.247 

10.06 

1.302 

3 

10.03 

1.172 

10.62 

1.243 

10.06 

1.307 

4 

10.31 

1.179 

24.97 

1.255 

25.02 

1.316 

5 

25.00 

1.183 

25.01 

1.256 

25.76 

1.320 

6 

25.02 

1.190 

25.06 

1.251 

40.06 

1.325 

7 

40.01 

1.201 

39.97 

1.264 

60.06 

1.340 

8 

59.96 

1.212 

39.99 

1.266 

80.06 

1.357 

9 

79.99 

1.232 

59.95 

1.280 

100.06 

1.375 

10 

100.01 

1.245 

60.02 

1.284 

120.06 

1.390 

11 

120.31 

1.259 

79.96 

1.301 

159.99 

1.426 

12 

120.34 

1.259 

80.21 

1.298 

199.99 

1.463 

13 

160.41 

1.295 

100.00 

1.312 

199.99 

1.465 

14 

200.23 

1.337 

119.96 

1.330 

239.99 

1.511 

15 

200.33 

1.332 

159.91 

1.373 

279.94 

1.552 

16 

240.53 

1.372 

199.85 

1.406 

319.99 

1.599 

17 

279.93 

1.419 

239.89 

1.450 

360.02 

1.653 

18 

319.94 

1.463 

278.63 

1.492 

399.72 

1.708 

19 

359.94 

1.509 

319.89 

1.540 

439.81 

1.768 

20 

359.98 

1.512 

319.89 

1.540 

479.73 

1.824 

21 

399.74 

1.565 

359.90 

1.590 

519.68 

1.893 

22 

439.88 

1.621 

400.02 

1.641 

559.7 

1.962 

23 

479.91 

1.674 

440.05 

1.700 

599.7 

2.042 

24 

519.92 

1.742 

480.03 

1.769 

639.7 

2.110 

25 

559.84 

1.814 

520.03 

1.832 

676.5 

2.194 

26 

599.69 

1.881 

560.04 

1.896 

719.7 

2.298 

27 

640.0 

1.963 

600.1 

1.973 

759.7 

2.398 

28 

680.0 

2.044 

640.2 

2.041 

799.6 

2.493 

29 

720.0 

2.130 

680.2 

2.130 

839.6 

2.602 

30 

760.0 

2.230 

720.2 

2.221 

879.7 

2.725 

31 

799.6 

2.327 

760.1 

2.307 

919.6 

2.851 

32 

839.5 

2.422 

800.2 

2.408 

959.6 

2.981 

33 

840.0 

2.422 

840.1 

2.524 

34 

879.9 

2.547 

880.1 

2.654 

35 

920.0 

2.661 

920.0 

2.777 

36 

960.1 

2.792 

960.1 

2.903 

37 

1000.0 

2.928 

984.6 

2.996 
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Table  XV 


Coefficients  of  X 

=  X  +  X.p  for  neon 
o  1 

(X  in  cal/cm.s 

•K;  p  in  amagat) 

p 

max 

n 

(X  +  a.  )xl0^ 

O  —  A 

O 

(X  +  o  )xl07 

1  “  A1 

ax 

°c 

am 

% 

25 

40 

7 

1.166  +  0.002 

0.86  +  0.10 

0.27 

25 

60 

8 

1.167  +  0.002 

0.79  +  0.06 

0.27 

25 

80 

9 

1.166  +  0.002 

0.81  +  0.04 

0.25 

25 

100 

10 

1.167  +  0.001 

0.79  +  0.03 

0.24 

25 

120 

12 

1.167  +  0.001 

0.77  +  0.02 

0.23 

50 

• 

40 

8 

1.237  +  0.002 

0.71  +  0.06 

0.18 

50 

60 

10 

1.236  +  0.001 

0.75  +  0.04 

0.18 

50 

80 

12 

1.235  +  0.001 

0.79  +  0.03 

0.19 

50 

100 

13 

1.235  +  0.001 

0.78  +  0.02 

0.18 

50 

120 

14 

1.235  +  0.001 

0.78  +  0.02 

0.18 

75 

40 

6 

1.297  +  0.002 

0.79  +  0.09 

0.22 

75 

60 

7 

1.297  +  0.002 

0.73  +  0.05 

0.20 

75 

80 

8 

1.297  +  0.001 

0.74  +  0.03 

0.19 

75 

100 

9 

1.297  +  0.001 

0.77  +  0.03 

0.19 

75 

120 

10 

1.296  +  0.001 

0.77  +  0.02 

0.18 
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Table  XVI 


Coefficients  of  X  =  \ 

2 

+  X.p  +  X.p  for 

O  1  i 

neon 

(X  in  cal/cm. 

s.K;  p  in  amagat) 

t 

p 

max 

n 

(X  +  a,  )xl04 

0  —  A 

0 

(X  +  a  )xl07 

1  -  xx 

(X2  +  ^)xl010 

CTx 

°c 

am 

% 

25 

200 

15 

1.168  +  0.002 

0.71  +  0.04 

0.57  + 

0.22 

0.26 

25 

240 

16 

1.168  +  0.002 

0.71  +  0.04 

0.57  + 

0.16 

0.25 

25 

280 

17 

1.168  +  0.001 

0.68  +  0.03 

0.72  + 

0.12 

0.25 

25 

320 

18 

1.169  +  0.001 

0.68  +  0.03 

0.75  + 

0.09 

0.'24 

25 

360 

20 

1.169  +  0.001 

0.68  +  0.02 

0.76  + 

0.06 

0.23 

25 

400 

21 

1.169  +  0.001 

0.67  +  0.02 

0.79  + 

0.05  . 

0.23 

25 

440 

22 

1.169  +  0.001 

0.66  +  0.02 

0.82  + 

0.04 

0.23 

50 

200 

16 

1.236  +  0.001 

0.71  +  0.04 

0.73  + 

0.20 

0.20 

50 

240 

17 

1.236  +  0.001 

0.71  +  0.03 

0.75  + 

0.14 

0.19 

50 

280 

18 

1.236  +  0.001 

0.71  +  0.03 

0.75  + 

0.10 

0.18 

50 

320 

20 

1.236  +  0.001 

0.71  +  0.02 

0.75  + 

0.06 

0.17 

50 

360 

21 

1.236  +  0.001 

0.71+0.02 

0.76  + 

0.05 

0.17 

50 

400 

22 

1.236  +  0.001 

0.71  +  0.02 

0.76  + 

0.04 

0.16 

50 

440 

23 

1.237  +  0.001 

0.70  +  0.01 

0.79  + 

0.04 

0.16 

75 

200 

13 

1.297  +  0.001 

0.69  +  0.03 

0.69  + 

0.16 

0.16 

75 

240 

14 

1.298  +  0.001 

0.66  +  0.03 

0.90  + 

0.14 

0.18 

75 

280 

15 

1.298  +  0.001 

0.66  +  0.03 

0.87  + 

0.10 

0.17 

75 

320 

16 

1.298  +  0.001 

0.67  +  0.02 

0.86  + 

0.07 

0.16 

75 

360 

17 

1.298  +  0.001 

0.66  +  0.02 

0.90  + 

0.05 

0.16 

75 

400 

18 

1.298  +  0.001 

0.65  +  0.02 

0.92  + 

0.04 

0.16 

75 

440 

19 

1.299  +  0.001 

0.64  +  0.01 

0.95  + 

0.03 

0.16 
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Table  XVII 

r2  3 

_p  +  X_p  for  neon 
o±'/  J 


(X 

In  cal/cm.s.K; 

p  In  amagat) 

t 

p 

max 

n 

(a  +  C 

o  — 

rx 

0 

)xl04 

(X.  +  a.  )xl07 
1  x1 

m  in 

(X,  +  ar  )xl0AU 
i  x2 

+  ay  )xl013 

J  x3 

ax 

°c 

am 

% 

25 

560 

25 

1.168 

+ 

0.001 

0.72  +  0.03 

0.44  +  0.13 

0.59  +  0.17 

0.22 

25 

600 

26 

1.168 

+ 

0.001 

0.71  +  0.02 

0.47  +  0.11 

0.54  +  0.13 

0.22 

25 

640 

27 

1.168 

+ 

0.001 

0.72  +  0.02 

0.45  +  0.10 

0.58  +  0.11 

0.21 

25 

680 

28 

1.168 

+ 

0.001 

0.72  +  0.02 

0.45  +  0.08 

0.58  +  0.09 

0.21 

25 

720 

29 

1.168 

+ 

0.001 

0.71  +  0.02 

0.46  +  0.07 

0.56  +  0.08 

0.20 

25 

760 

30 

1.168 

+ 

0.001 

0.72  +  0.02 

0.43  +  0.07 

0.61  +  0.06 

0.20 

25 

800 

31 

1.168 

+ 

0.001 

0.72  +  0.02 

0.42  +  0.06 

0.61  +  0.05 

0.20 

50 

560 

26 

1.236 

+ 

0.001 

0.74  +  0.03 

0.51  +  0.12 

0.53  +  0.15 

0.18 

50 

600 

27 

1.236 

+ 

0.001 

0.73  +  0.02 

0.51  +  0.10 

0.52  +  0.10 

0.18 

50 

640 

28 

1.236 

+ 

0.001 

0.72  +  0.02 

0.61  +  0.10 

0.38  +  0.11 

0.19 

50 

680 

29 

1.236 

+ 

0.001 

0.72  +  0.02 

0.59  +  0.09 

0.41  +  0.09 

0.19 

50 

720 

30 

1.236 

+ 

0.001 

0.73  +  0.02 

0.54  +  0.08 

0.47  +  0.08 

0.19 

50 

760 

31 

1.236 

+ 

0.001 

0.73  +  0.02 

0.56  +  0.07 

0,45  +  0.06 

0.19 

50 

800 

32 

1.236 

+ 

0.001 

0.73  +  0.02 

0.55  +  0.06 

0.46  +  0.05 

0.18 

75 

560 

22 

1.298 

+ 

0.001 

0.67  +  0.02 

0.77  +  0.11 

0.15 

75 

600 

23 

1.298 

+ 

0.001 

0.69  +  0.02 

0.69  +  0.11 

0.38  +  0.11 

0.15 

75 

640 

24 

1.298 

+ 

0.001 

0.67  +  0.02 

0.77  +  0.09 

0.27  +  0.10 

0.16 

75 

680 

25 

1.298 

+ 

0.001 

0.68  +  0.02 

0.72  +  0.08 

0.33  +  0.08 

0.16 

75 

720 

26 

1.297 

+ 

0.001 

0.70  +  0.02 

0.64  +  0.08 

0.44  +  0.08 

0.18 

75 

760 

27 

1.297 

+ 

0.001 

0.71  +  0.02 

0.57  +  0.07 

0.52  +  0.07 

0.19 

75 

800 

28 

1.297 

+ 

0.001 

0.71  +  0.02 

0.57  +  0.06 

0.52  +  0.06 

0.19 
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Table  XVIII 


Coefficients 

of  X=X  +  X.p  + 
o  1 

2  2 
X£p  Jlnp  +  X^P 

for  neon 

(X 

in  cal/cm.s.K; 

p  in  amagat) 

t 

Pmax 

n 

a  +  a.  )xl04 
o  —  X 

0 

(X1  +  a.  )xl0^ 
1  *1 

(X£  +  jx,)x1010 

(X2  +  )xl010 

ax 

c 

am 

% 

25 

240 

16 

1.166  +  0.002 

0.9  +  0.1 

[1.1  +  0.8] 

[-6  +  5] 

0.24 

25 

320 

18 

1.167  +  0.002 

0.83  +  0.09 

[0.7  +  0.4] 

[-4  +  3] 

0.23 

25 

400 

21 

1.167  +  0.001 

0.79  +  0.07 

[0.4  +  0.2] 

[-2  +  2] 

0.22 

25 

480 

23 

1.168  +  0.001 

0.75  +  0.05 

[0.3  +0.2] 

[-1  +  1] 

0.22 

25 

520 

24 

1.167  +  0.001 

0.76  +  0.04 

0.31  +  0.12 

[-1.3  +  0.9] 

0.21 

25 

560 

25 

1.167  +  0.001 

0.78  +  0.04 

0.39  +  0.11 

-1.8  +  0.7 

0.21 

25 

600 

26 

1.167  +  0.001 

0.78  +  0.04 

0.38  +  0.09 

-1.7  +  0.6 

0.21 

25 

640 

27 

1.167  +  0.001 

0.79  +  0.03 

0.42  +  0.08 

-2.0  +  0.6 

0.21 

25 

680 

28 

1.167  +  0.001 

0.79  +  0.03 

0.44  +  0.07 

-2.1  +  0.5 

0.20 

25 

720 

29 

1.167  +  0.001 

0.80  +  0.03 

0.45  +  0.06 

-2.2  +  0.4 

0.20 

50 

240 

17 

1.236  +  0.002 

0.7  +  0.1 

[-0.1  +  0.6] 

[1  +  5] 

0.20 

50 

320 

20 

1.236  +  0.001 

0.70  +  0.07 

[-0.0  +  0.3] 

[1  +  2] 

0.18 

50. 

400 

22 

1.236  +  0.001 

0.71  +  0.05 

[  0.0  +  0.2] 

[1  ±  1] 

0.17 

50 

480 

24 

1.235  +  0.001 

0.78  +  0.05 

0.33  +  0.14 

[-1.4  +  0.9] 

0.19 

50 

520 

25 

1.235  +  0.001 

0.79  +  0.04 

0.37  +  0.11 

[-1.6  +  0.8] 

0.19 

50 

560 

26 

1.235  +  0.001 

0.78  +  0.04 

0.32  +  0.10 

[-1.3  +0.7] 

0.18 

50 

600 

27 

1.235  +  0.001 

0.79  +  0.04 

0.34  +  0.09 

-1.4  +  0.6 

0.18 

50 

640 

28 

1.236  +  0.001 

0.77  +  0.03 

0.28  +  0.08 

-1.0  +  0.5 

0.19 

50 

680 

29 

1.235  +  0.001 

0.78  +  0.03 

0.31  +  0.07 

-1.3  +  0.5 

0.19 

50 

720 

30 

1.235  +  0.001 

0.80  +  0.03 

0.36  +  0.06 

-1.6  +  0.5 

0.19 

75 

240 

14 

1.296  +  0.002 

0.8  +  0.1 

[1.1  +  0.6] 

[-6  +  4] 

0.16 

75 

320 

16 

1.297  +  0.001 

0.71  +  0.08 

[0.2  +  0.3] 

[-0  +  2] 

0.17 

75 

400 

18 

1.297  +  0.001 

0.73  +  0.05 

[0.3  +  0.2] 

[-1  +  1] 

0.16 

75 

480 

20 

1.298  +  0.001 

0.70  +  0.04 

[0.2  +  0.1] 

[-0.1  +  0.9] 

0.16 

75 

520 

21 

1.297  +  0.001 

0.70  +  0.04 

[0.2  +  0.1] 

[-0.2  +  0.7] 

0.15 

75 

560 

22 

1.297  +  0.001 

0.71  +  0.04 

[0.2  +  0.1] 

[-0.3  +  0.6] 

0.15 

75 

600 

23 

1.297  +  0.001 

0.73  +  0.03 

0.26  +  0.08 

[-0.8  +  0.6] 

0.16 

75 

.640 

24 

1.297  +  0.001 

0.71  +  0.03 

0.21  +  0.07 

[-0.5  +  0.5] 

0.16 

75 

680 

25 

1.297  +  0.001 

0.73  +  0.03 

0.25  ±  0.06 

[-0.8  +  0.5] 

0.16 

75 

720 

.26 

1.297  ±  0-001 

0.76  +  0.03 

0.34  +  0.07 

-1.4  +  0.5 

0.19 

Note:  When  X^  and/or  are  not  significant  individually  we  place  the  value  between 


brackets 
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4.4  THE  FIRST  DENSITY  COEFFICIENT  OF  THERMAL  CONDUCTIVITY  AND 
VISCOSITY  FOR  A  NUMBER  OF  GASES 


Hanley,  McCarty  and  the  author  applied  the  techniques  described 

in  Sections  4.2  and  4.3  to  the  thermal  conductivity  and  viscosity  of  a 

number  of  gases  and  tried  to  determine  experimental  first  density  coef- 

33 

ficients  X^  and  Unfortunately,  data  for  a  single  gas  that  are 

accurate  and  extensive  enough  to  yield  reliable  values  for  this  density 

coefficient  over  a  wide  temperature  range  do  not  exist.  Therefore,  to 

obtain  a  more  general  picture  of  the  temperature  dependence  of  the  first 

density  coefficient,  one  needs  to  coordinate  the  results  obtained  for 

different  gases  using  the  principle  of  corresponding  states.  For  this 

purpose  we  considered  the  noble  gases:  helium,  neon,  argon,  krypton  and 

xenon.  Hydrogen  and  nitrogen  were  also  included  to  obtain  additional 

information  at  more  extreme  reduced  temperatures. 

Whenever  possible,  we  prepared  a  set  of  tables  similar  to  those 

shown  for  the  thermal  conductivity  of  neon,  but,  of  course,  such  a  complete 

analysis  was  only  possible  in  a  limited  number  of  cases.  In  most  cases  we 

again  noted  a  tendency  for  the  logarithmic  equation  (4.2-5)  to  be  more 

consistent  with  the  experimental  data  than  a  pure  polynomial.  In  particular, 

39 

the  data  of  LeNeindre  et.al.  exhibited  parallel  behavior  to  that  shown 
for  neon.  In  many  cases  the  scatter  of  the  data  was  large  enough  to  mask 
the  difference  between  (4.2-5)  and  a  polynomial,  but,  even  then,  we 
preferred  the  equations  (4.2-4)  and  (4.2-5),  since  they  contain  the  first 
terms  of  the  theoretically  predicted  density  expansion. 

The  final  estimates  for  and  are  presented  in  Tables  XIX  and 
XX.  The  references  to  the  sources  of  the  experimental  data  used  in  the 
analysis  are  given  in  the  first  column.  For  experimental  data  that  were 
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given  as  a  function  of  pressure,  the  pressures  were  converted  to  densities. 

A  reference  to  the  equation  of  state  data,  when  such  a  conversion  was  needed, 
is  also  included.  Column  6  indicates  whether  the  final  estimate  of  A^  and 
H.  is  based  on  the  linear  equation  (4.2-4)  or  the  logarithmic  equation 
(4.2-5).  This  choice  was  usually  dictated  by  the  number  and  the  range  of 
the  data  at  a  given  isotherm. 

The  estimated  A^  and  together  with  their  standard  deviation  are 
listed  in  column  4  of  Tables  XIX  and  XX.  The  values  for  the  corresponding 
95%  confidence  interval  are  included  in  column  5.  When  the  data  were 
extensive  enough  to  permit  a  test  of  equation  (4.2-5),  empirical  values 
were  obtained  for  the  higher  order  coefficients.  However,  they  have  not 
been  included  here,  since,  as  mentioned  in  the  previous  Section,  their 
physical  meaning  is  uncertain. 

If  the  molecules  interact  with  a  spherically  symmetric  two-parameter 
potential,  if  the  internal  degrees  of  freedom  of  the  molecules  do  not  contri¬ 
bute  to  the  transport  of  momentum  and  energy,  and  if  the  translational 
motion  of  the  molecules  can  be  described  by  classical  mechanics,  the 

transport  properties  in  dimensionless  form  should  be  universal  functions  of 

2 

the  reduced  state  variables.  Thus,  A^  and  of  the  noble  gases  should 

obey  a  law  of  corresponding  states,  certainly  within  the  precision  with 

6  58 

which  these  coefficients  can  be  determined  at  the  present.  *.  _  One  may 

also  expect  that  the  corresponding  states  relation  for  the  transport 

properties  of  the  noble  gases  can  be  extended  to  include  the  viscosity  of 

33  59 

nitrogen  without  introducing  too  large  an  error.  ’  However,  it  is 

doubtful  whether  hydrogen  will  satisfy  a  corresponding  states  relation 

59 

with  the  noble  gases.  So  any  results  for  hydrogen  should  be  interpreted 
with  caution. 
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Dimensionless  first  density  coefficients  are  obtained  by  defining 


A*  =  A,  — 
1  1  rp 


STP 


(4.4-1) 


where  m  is  the  molecular  mass ,  and  r  and  e  are  length  and  energy  reducing 

3  60 

parameters  deduced  for  the  Lennard-Jones  potential.  ’  The  values  of 
these  reduced  first  density  coefficients  are  presented  in  the  last  column 
of  Tables  XIX  and  XX,  as  a  function  of  the  reduced  temperature  T*  =  kT/e. 
The  reduction  parameters  used  are  given  in  Table  XXI. 
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Table  XIX 

The  first  density  coefficient  X^  of  thermal  conductivity 


ref 

gas 

T 

(1^  + 

a.  )xl07 

A1 

(*x  ± 

A.  )xl07 
X1 

eq. 

T*  X*  +  l 

K 

cal/cm, 

s.K.am 

39 

He 

303.0 

2.33 

+  0.08 

2.33 

+  0.15 

(4.2-5) 

30.3 

6.94 

+  0.45 

40,41 

He 

195.2 

1.22 

+  0.07 

1.22 

+  0.14 

(4.2-4) 

19.5 

3.63 

+  0.42 

126.6 

1.40 

+  0.08 

1.40 

+  0.15 

(4.2-4) 

12.7 

4.17 

+  0.45 

77.9 

0.99 

+  0.06 

0.99 

±  0.26 

(4.2-4) 

7.79 

2.95 

+  0.77 

43.15 

1.01 

+  0.11 

1.01 

+  0.24 

(4.2-5) 

4.32 

3.01 

+  0.71 

42 

He 

260.0 

2.86 

+  0.07 

2.86 

+  0.17 

(4.2-4) 

26.0 

8.52 

+  0.51 

240.0 

2.51 

+  0.22 

2.51 

+  0.61 

(4.2-4) 

24.0 

7.45 

+  1.82 

220.0 

2.35 

+  0.20 

2.35 

+  0.51 

(4.2-4) 

22.0 

7.00 

+  1.52 

200.7 

2.05 

+  0.15 

2.05 

+  0.48 

(4.2-4) 

20.1 

6.11 

+  1.43 

180.0 

1.37 

+  0.06 

1.37 

+  0.15 

(4.2-4) 

18.0 

4.08 

+  0.45 

100.5 

0.88 

+  0.02 

0.88 

+  0.05 

(4.2-4) 

10.1 

2.62 

+  0.15 

87.5 

0.91 

+  0.04 

0.91 

+  0.11 

(4.2-4) 

8.75 

2.71 

+  0.33 

21.15 

0.442  +  0.003 

0.44 

+  0.01 

(4.2-4) 

2.12 

1.31 

+  0.03 

37 

Ne 

348.15 

0.76 

+  0.03 

0.76 

+  0.06 

(4.2-5) 

7.41 

2.26 

+  0.18 

323.15 

0.79 

±  0.04 

0.79 

+  0.08 

(4.2-5) 

6.87 

2.35 

+  0.24 

298.15 

0.79 

+  0.03 

0.79 

+  0.07 

(4.2-5) 

6.34 

2.35 

+  0.21 

•43 

Ar 

831.15 

0.99 

+  0.07 

0.99 

+  0.22 

(4.2-4) 

6.65 

2.02 

+  0.45 

775.15 

0.84 

+  0.05 

0.84 

+  0.16 

(4.2-4) 

6.20 

1.72 

+  0.33 

683.15 

1.18 

+  0.02 

1.18 

+  0.06 

(4.2-4) 

5.47 

2.41 

+  0.12 

674.15 

1.02 

+  0.09 

1.02 

+  0.19 

(4.2-5) 

5.39 

2.09 

+  0.39 

475.15 

1.24 

+  0.07 

1.24 

+  0.15 

(4.2-5) 

3.80 

2.54 

+  0.31 

413.15 

1.18 

+  0.08 

1.18 

+  0.16 

(4.2-5) 

3.31 

2.41 

+  0.33 

370.15 

1.20 

+  0.04 

1.20 

+  0.09 

(4.2-4) 

2.96 

2.45 

+  0.18 
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Table  XIX  (continued) 

The  £irst  density  coefficient  X^  of  thermal  conductivity 


ref 

gas 

T 

(x1  + 

a,  )xl0^ 

*1 

(Xi  + 

A,  )xl0^ 
A1 

eq. 

T* 

X*  +  l 

*X* 

A1 

K 

cal/cm. 

S.K.  dID 

A3 

Ar 

298.15 

1.06 

+  0.08 

1.06 

+  0.17 

(4.2-5) 

2.39 

2.17 

+  0.35 

44 

Ar 

348.15 

1.00 

+  0.06 

1.00 

+  0.13 

(4.2-5) 

2.79 

2.05 

+  0.27 

323.15 

1.02 

+  0.08 

1.02 

+  0.34 

(4.2-4) 

2.59 

2.09 

+  0.70 

' 

298.15 

1.02 

+  0.08 

1.02 

+  0.30 

(4.2-4) 

2.39 

2.09 

+  0.61 

45,46 

Ar 

193.15 

1.03 

+  0.04 

1.03 

+  0.50 

(4.2-4) 

1.55 

2.11 

+  1.02 

183.15 

1.15 

+  0.04 

1.15 

+  0.15 

(4.2-4) 

1.47 

2.35 

+  0.31 

173.15 

1.71 

+  0.12 

1.71 

+  0.31 

(4.2-4) 

1.40 

3.50 

+  0.63 

153.15 

1.95 

+  0.28 

1.95 

+  3.56 

(4.2-4) 

1.23 

3.99 

+  7.28 

39 

N2 

298.15 

1.21 

+  0.06 

1.21 

±  0.12 

(4.2-5) 

3.28 

2.25 

+  0.22 

47,48 

N2 

973.15 

1.36 

+  0.29 

1.36 

+  0.68 

(4.2-5) 

10.7 

2.53 

+  1.27 

873.15 

1.41 

+  0.20 

1.41 

+  0.47 

(4.2-5) 

9.83 

2.63 

+  0.88 

673.15 

1.35 

+  0.01 

1.35 

+  0.13 

(4.2-4) 

7.41 

2.51 

+  0.24 

573.15 

1.22 

+  0.07 

1.22 

+  0.22 

(4.2-5) 

6.31 

2.27 

+  0.41 

473.15 

1.25 

+  0.01 

1.25 

+  0.17 

(4.2-4) 

5.21 

2.33 

+  0.32 

398.15 

1.34 

+  0.01 

1.34 

+  0.13 

(4.2-4) 

4.38 

2.50 

+  6.24 

348.15 

1.31 

+  0.02 

1.31 

+  0.20 

(4.2-4) 

3.83 

2.44 

+  0.37 

45,48 

N2 

200.0 

0.98 

+  0.14 

0.98 

+  0.44 

(4.2-5) 

2.20 

1.82 

+  0.82 

180.0 

0.89 

+  0.22 

0.89 

+  0.70 

(4.2-5) 

1.98 

1.66 

+  1.30 

170.0 

0.96 

+  0.04 

0.96 

+  0.51 

(4.2-4) 

1.87 

1.79 

+  0.95 

160.0 

1.02 

+  0.05 

1.02 

+  0.62 

(4.2-4) 

1.76 

1.90 

+  1.15 

130.0 

2.61 

+  0.25 

2.61 

+  1.07 

(4.2-4) 

1.43 

4.86 

+  1.99 
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Table  XX 

The  first  density  coefficient  of  viscosity 


Ref 

gas 

T 

(n,  +  a  )xl07 
1  ~  nl 

'  (n.  +  A  )xl07 
1~  nl 

J. 

eq. 

T* 

(n*  +  A  *)xl0 
1  —  ri* 

K 

g/cm. 

iS.am 

49 

He 

293.15 

-0.024  +  0.007 

-0.024  +0.015 

(4.2-4) 

29.3 

-0.36  +  0.22 

50 

He 

373.15 

-0.21  +  0.03 

-0.21  +  0.13 

(4.2-4) 

37.3 

-3.10  +  1.92 

298.15 

-0.05  +  0.01 

-0.05  +  0.13 

(4.2-4) 

29.8 

-0.74  +  1.92 

223.15 

0.021  +0.005 

0.021  +  0.020 

(4.2-4) 

22.3 

0.31  +  0.30 

51 

Ne 

348.15 

0.48  +  0.07 

0.48  +  0.14 

(4.2-5) 

7.41 

1.36  +  0.40 

323.15 

0.43  +  0.04 

0.43  +  0.13 

(4.2-4) 

6.88 

1.22  +  0.37 

298.15 

0.55  +  0.05 

0.55  +  0.14 

(4.2-4) 

6.34 

1.56  +  0.40 

50 

Ne 

373.15 

0.52  +  0.05 

0.52  +  0.64 

(4.2-4) 

7.93 

1.47  +  1.81 

298.15 

0.66  +  0.04 

0.66  +  0.50 

(4.2-4) 

6.34 

1.87  +  1.41 

248.15 

0.68  +  0.05 

0.68  +  0.20 

(4.2-4) 

5.28 

1.92  +  0.56 

223.15 

0.76  +  0.07 

0.76  +  0.89 

(4.2-4) 

4.75 

2.15  +  2.52 

49 

Ne 

293.15 

0.48  +  0.07 

0.48  +  0.19 

(4.2-4) 

6.24 

1.36  +  0.54 

49 

Ar 

298.17 

2.15  +  0.07 

2.15  +  0.30 

(4.2-4) 

2.39 

2.19  +  0.31 

293.17 

2.03  +  0.09 

2.03  +  0.22 
__  f 

(4.2-5) 

2.35 

2.07  +  0.22 

52 

Ar 

298.15 

2.20  +  0.05 

2.20  +  0.12 

(4.2-5) 

2.39 

2.24  +  0.12 

50 

Ar 

298.15 

2.14  +  0.30 

2.14  +  1.29 

(4.2-5) 

2.39 

2.18  +  1.31 

248.0 

2.13  +  0.50 

2.13  +  2.15 

(4.2-5) 

1.98 

2.17  +  2.19 

223.0 

2.48  +  0.50 

2.48  +  2.15 

(4.2-5) 

1.78 

2.53  +  2.19 

53 

Ar 

348.15 

2.80  +  0.08 

2.80  +  0.34 

(4.2-4) 

2.79 

2.85  +  0.35 

323.15 

2.74  +  0.40 

2.74  +  0.92 

(4.2-5) 

2.59 

2.79  +  0.94 

298.15 

2.31  +  0.20 

2.31  +  0.46 

(4.2-5) 

2.39 

2.35  +  0.47 
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Table  XX  (continued) 

The  first  density  coefficient  of  viscosity 


Ref 

gas 

T 

(n.  +  a  )xl07 
!-  n1 

(n.  +  A  )  x  107 
l-  n1 

eq. 

T* 

(n?  +  A  .  )xl0 
l  —  n  j 

K 

g/cm. 

,  s.am 

54 

Kr 

348.15 

3.69 

+  0.30 

3.69  +  0.67 

(4.2-5) 

1.90 

2.02  +  0.37 

323.15 

4.81 

+  0.24 

4.81  +  0.76 

(4.2-4) 

1.77 

2.63  +  0.42 

298.15 

3.36 

+  0.19 

3.36  +  0.41 

(4.2-5) 

1.63 

1.84  +  0.22 

49 

Kr 

293.15 

3.47 

+  0.15 

3.47  +  0.64 

(4.2-5) 

1.60 

1.90  +  0.35 

49 

Xe 

298.15 

3.36 

+  0.30 

3.36  +  3.81 

(4.2-4) 

1.19 

1.14  +  1.30 

50 

N2 

373.15 

2.08 

+  0.15 

2.08  +  1.90 

(4.2-4) 

4.11 

2.75  +  2.51 

298.15 

1.18 

+  0.14 

1.18  +  0.36 

(4.2-5) 

3.28 

1.56  +  0.48 

248.15 

2.10 

+  0.30 

2.10  +  1.29 

(4.2-5) 

2.73 

2.78  +  1.71 

223.15 

1.95 

+  0.03 

1.95  +  0.13 

(4.2-5) 

2.45 

2.58  +  0.17 

49 

N2 

298.15 

1.48 

+  0.10 

1.48  +  0.32 

(4.2-5) 

3.28 

1.96  +  0.42 

293.15 

1.36 

+  0.25 

1.36  +  0.59 

(4.2-5) 

3.22 

1.80  +  0.78 

55 

N2 

323.15 

1.61 

+  0.30 

1.61  +  1.29 

(4.2-5) 

3.56 

2.13  +  1.71 

298.15 

1.70 

+  0.07 

1.70  +  0.14 

(4.2-5) 

3.28 

2.25  +  0.19 

56 

E2 

423.15 

0 

+  ? 

0  +  ? 

- 

12.8 

0  +  ? 

348.15 

0.16 

+  0.013 

0.16  +  0.05 

(4.2-4) 

10.6 

1.68  +  0.52 

57 

H2 

348.15 

0.173  +  0.004 

0.173  +  0.016 

(4.2-4) 

10.6 

1.82  +  0.17 
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Table  XXI 


Reduction  parameters  used  to  express  1*  and  nj  in  dimensionless  units. 


gas 


PSTP 


r 

o 

A 


K 


cal/cm.s .K.am 


g/cm.s.am 


He 

0.00017845 

2.63 

10.0 

2.979 

14.78 

Ne 

0.00089979 

2.72 

47.0 

2.980 

2.828 

Ar 

0.0017834 

3.41 

125.0 

2.045 

1.019 

Kr 

0.0037481 

3.62 

183.0 

2.306 

0.547 

Xe 

0.0058974 

3.96 

250.0 

,2.248 

0.340 

N2 

0.0012506 

3.68 

90.9 

1.862 

1.323 

H2 

0.00008994. 

2.87 

33.0 

1.065 

10.511 
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4.5  DISCUSSION  OF  THE  RESULTS 


We  have  proposed  some  criteria  to  check  whether  a  given  density 
function  is  a  possible  correct  representation  of  experimental  data  for  the 
transport  properties  of  gases.  According  to  these  criteria,  the  density 
function  suggested  by  the  theory,  which  Involves  a  term  logarithmic  in  the 
density,  is  shown  to  be  consistent  with  experimental  data. 

In  particular,  by  investigating  whether  the  coefficient  of  the 
linear  term  does  not  change  upon  adding  additional  terms  to  the  linear 
equation,  we  verified  whether  meaningful  values  for  the  first  density  co¬ 
efficient  are  deduced  from  the  experimental  transport  data.  Based  on  this 
analysis  we  report  experimental  values  for  the  first  density  coefficients 
X^  and  for  a  number  of  gases.  The  reduced  first  density  coefficients. 

X*  and  are  plotted  in  Figs.  13  and  14  as  a  function  of  the  reduced 

temperature  T*,  up  to  T*  =  30. 

Both  X*  and  nj  appear  insensitive  to  the  temperature  in  the  approximate 
range  2  <  T*  <  10.  The  theoretical  predictions  for  X^  and  developed  in 
Chapter  II  were  limited  to  a  gas  of  hard  spheres  and,  therefore,  cannot  be 
used  legitimately  to  describe  X^  and  for  a  real  gas.  Nevertheless,  we 
could  make  such  an  attempt  by  attributing  an  effective  hard  sphere  diameter 
a  to  the  real  molecules.  The  most  natural  choice  for  such  a  diameter  a  is 
the  requirement  that  the  dilute  gas  values  Xq  and  nQ  should  reduce  to  the 
hard  sphere  expressions  (2.2-13).  For  such  a  hypothetical  hard  sphere  gas 
the  first  density  coefficients  are  given  by  (2.8-1)  (together  with  the 
appropriate  collisional  transfer  contribution) .  If  we  apply  this  procedure 
at  temperatures  corresponding  to  T*  =  4  we  obtain  X*  'v  1.9  for  helium, 
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! 


Figure  13.  The  reduced  first  density  coefficient  X*  of  thermal  conductivity 
as  a  function  of  the  reduced  temperature  T*. 
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Figure  14.  The  reduced  first  density  coefficient  nj  of  viscosity  as  a 

function  of  the  reduced  temperature  T*. 
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A*  'v  1.8  for  neon,  A*  •v  1.7  for  argon  and  A*  ^  1.8  for  nitrogen. 
Corresponding  values  obtained  for  rj  are  ^  0.17  for  helium, 
nj  y  0.16  for  neon,  n*  ^  0.15  for  argon,  n*  ^  0.17  for  krypton, 
nj  ^  0.17  for  xenon,  nj  ^  0.15  for  nitrogen  and  n*  %  0.13  for 
hydrogen . 

From  Figs.  13  and  14  we  see  that  these  estimates  have  a 

reasonable  order  of  magnitude  in  the  insensitive  temperature  range. 

The  first  density  coefficients  appear  sensitive  to  temperature 

when  T*  <  2  and  when  T*  >  10.  In  fact,  it  is  clear  that  nj  will  pass 
7  58 

through  zero.  *  Since  this  observation  is  based  largely  on  experimental 

i  ■■ 

data  for  helium  (recently  Kao  and  Kobayashi  have  presented  further  evi¬ 
dence  that  n1  becomes  negative  for  helium*^),  some  estimated  values  for 
rij  for  hydrogen"^*"^  are  included  in  Fig.  14.  However,  it  does  not  seem 
possible  to  determine  accurately  the  reduced  temperature  at  which  nj  =  0. 

By  contrast,  Fig.  13  shows  that  in  the  same  temperature  range  A*  does  not 
change  sign. 

The  existence  of  the  insensitive  temperature  range  means  that  it 

is  Imperative  to  have  reliable  experimental  data  for  temperatures  outside 

the  range,  if  we  are  to  make  any  sensitive  comparisons  of  theory  with 

experiment.  At  this  time  we  rely  heavily  on  the  data  of  Ziebland  and 
45 

Burton  for  A*  at  low  reduced  temperatures,  and  we  are  very  dependent 
on  the  preliminary  work  at  NBS  for  A*  at  high  temperatures.  It  Is 
obvious  from  the  spread  of  the  confidence  Interval  bars  in  Figs.  13  and 
14,  that  the  coefficients  A^  and  should  be  known  with  a  better  precision. 
Thus,  more  reliable  data  are  required  for  the  whole  temperature  range.  The 
precision  of  A^  and  can  only  be  improved  if  the  individual  data  points 
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are  determined  with  a  standard  deviation  of  about  0.5Z  for  X  and  0.2% 

or  less  for  n>  Furthermore,  the  data  points  on  an  individual  Isotherm 

34 

should  be  numerous  and  well-spaced,  in  particular,  in  the  linear 
region.  It  goes  without  saying  that  systematic  errors  should  be 
minimized  as  much  as  possible. 
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APPENDIX  A.  SURFACE  INTEGRAL  FORM  OF  THE  TRIPLE  COLLISION  INTEGRAL 


The  fact  that  the  six  dimensional  configurational  triple  collision 


integral  can  be  converted  into  a  five  dimensional  surface  integral  was 
pointed  out  by  Green. ^ 

The  streaming  operators  «f(l.  .  .s)  are  mathematically  defined  as^ 

„  ,  -TO?"  tr/(0) 

0  (1.  .s)  =  lim  eT  (1.  .s)  =  lim  e  S.e  8 


where  ft (1.  .s)  is  a  "Hamilton" 


operator 


(A-l) 


s  .  a  s 

(1.  .s)  =  L  v.  .  -  T.  0  .  . 

S  i=l  1  9ri  i<j  ±j  ' 


iy<°> 

s 

of  s  particles 


and^T  '''/  that  part  of  which  is  associated  with  the  free  streaming 

5  S 


oC  (1.  .s)  =  E  v  .  • 

s  ,  .  x  or. 

i=J  i 

We  assume  that  for  large  values  of  r  the  streaming  operators 
become  independent  of 


(A- 2) 


(A- 3) 


lim  4—  <T„(1.  .s)  -Jjf(l.  .s)«/\l.  ,s)  -cT(l.  .s)  Q.  .s)  =  0  ,  (A-4) 

X-Xo  11  S 


s 


so  that 


62 


(o) 


(1.  .s)/(l.  .s)  -ef(l.  .s)#f'(o)(l.  .s)  =  E  0..(1,  . 

S  A  ✓  A  1 J 


(A-5) 


i<j 


+)  In  this  assumption  we  neglect  the  collisions  which  are  taking  place 
exactly  at  the  time  -t.  Therefore,  we  make  an  error,  referred  to  as  the 
"initial"  error  term  by  Dorfman  and  Cohen.  This  initial  error  term  is  less 
divergent  than  the  asymptotic  error  term  discussed  in  Section  3.1  and  does 
not  contribute  to  the  first  density  coefficients,  nor  to  the  coefficient  of 
the  logarithm. 
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In  particular 


(123) 


£ (123)  -  £ (123)  3*0)  123)  -  (ei2  +  813  +  623)<f{123)  ,  (A-6) 


*  (123)  tf* (a)  -/(a)^<o)(123)  =  0„J*(a)  , 


(A-7) 


where  a  is  either  12,  13  or  23.  Using  these  relations  the  operator  0(123) 
given  by  (2.4-15)  can  be  rewritten  into  the  form^ 


0(123  =  /dr21dr‘31^3(0)U(123)  -  U(123),^(o) 

-  Z  £{»tf(o)U(a)  -  U(a)^£(o)}U(0)]  . 
a*3  J  J  J 


(A-8) 


For  convenience  we  have  introduced  the  Ursell  operators 

U(a)  =«f(a)  -  1  , 


(A-9) 


U(123)  =/(123)  -£(  12)  -«f(13)  -  *T (23)  +  2  . 

The  summations  over  a  and  3  are  to  be  taken  over  the  pairs  12,  13  and  23. 

The  integrand  depends  only  on  the  relative  positions  and  the  integration 
can  be  carried  out  over  any  two  relative  position  vectors,  for  which  we  have 
chosen  ?2^  and  Since  in  (2.4-14)  the  operator  0(123)  operates  on  a 


function  of  the  velocities  alone,  the  term  U(123) 


vanishes  in  (A-8) . 


In  order  to  define  an  appropriate  transformation  for  the  integration 

«*  ^  _k 

variables,  we  consider  for  fixed  V^,  V2>  V3  the  positions  of  particles  1,  2 
and  3  which  are  related  to  each  other  by  free  particle  motion.  This  is 
schematically  illustrated  in  Fig.  Al,  where  the  time  t  associated  with  the 
free  motion  of  the  particles  is  taken  to  increase  when  tne  diagram  is  read 
from  bottom  to  top. 
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Figure  A  1.  Diagram  representing  the  free  trajectories  of  particles  1,  2  and  3. 


However,  since  the  £  operators  transform  the  present  phases  of  the  particles 
to  phases  in  the  previous  history  of  the  motion,  we  shall  retrace  the 
trajectories  backwards  in  time  and  always  read  the  diagrams  from  top  to 
bottom  as  in  the  main  text. 

We  remark  that  the  three  free  particle  trajectories  have  at  least 
one  intersection",  by  which  we  mean^ that  along  the  free  trajectories  the 
relative  distance  of  at  least  one  pair  of  particles  becomes  equal  to  o  for 
some  value  of  t.  Otherwise,  the  particles  are  not  dynamically  correlated 
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and  the  Ursell  operators  in'  (A-8)  would  vanish  for  all  positions'  along  these 
trajectories.  Suppose  that  when  the  trajectories  in  Fig.  A1  are  retraced 
from  top  to  bottom,  the  distance  is  the  first  distance  which  becomes 
equal  to  c.  Then,  the  positions  of  the  particles  along  the  free  tra¬ 
jectories  can  be  characterized  by 


*21  =  **21^  =  W  =  r^i  (0)  +  tV^1  f 

*31  =  *31 (t)  =  exp(t^(o))  ^(0)  =  r31(0)  +  tV31  , 


(A-10) 


where  the  notation  r^CO)  indicates  that  the  distance  r91 (t)  is  equal  to 


21 


c  for  t  =  0.  We  take  r^CO),  r^(0)  and  t  as  the  new  integration  variables. 
In  this  new  coordinate  system  d^^CO)  is  a  two  dimensional  surface  element. 
The  Jacobian  of  this  transformation  is 


d*21d*31  =  V21'*21^d*21^d*31^dt  * 


(A-ll) 


From  the  expression  (2.4-15)  for  0(123)  we  see  that  the  integrand 
will  give  contributions  in  the  three  regions  for  which,  respectively,  r2i<0> 
r3^<o  and  r32<o.  Therefore,  the  upper  limit  for  the  integration  over  t  can 
be  taken  at  any  value  t  >  0  corresponding  to  phases  of  receding  particles. 

The  lower  limit  has  to  be  taken  at  a  sufficiently  negative  value  of  t,  so 
that  the  free  trajectories  will  no  longer  intersect  when  retraced  further 
backwards  in  time.  As  a  consequence,  all  Ursell  operators  vanish  at  the 
lower  limit  of  the  t-integration. 

In  this  new  coordinate  system 

Jf^(o)  _  3_  (A-12) 

^3  3t 
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since  for  any  function  g  (r21’^31’Vl’^2,V3>: 


-  7?  7?  7? 


<o)  T'  3 


gCr21’r31,Vl*V2’V3)  =  '  g(r21(0),r31(0),V1,V2,V3) 

t /(0)  * 

8t  e  g(r21(0),r31(0),V1,V2,V3)  =  —  8(r2i»r3i»vi»v2*v3) 


Therefore,  as  pointed  out  by  Green, ^  equation  (A-8)  transforms  to 

*  ^  *  (A"13) 
0(123)  =  /drfc21(0)d?31(0)V21.r21(0)  /dt  F|^UC123)  -  l  Z{f^U(a)  -  U(a)|^-}U(e)J  . 

In  the  first  term  of  (A-13)  the  t-integration  can  be  carried  out  immediately. 
The  second  term  can  be  simplified  by  partial  integration. 

0(123)  =  /df21(0)df^1(0)V21.?^1(0)  [u(123)  -  Z  Z/dt{|^U(a)Ju(B)]  .  (A-14) 

To  elucidate  the  result  for  the  second  time  integral  we  shall 
consider  as  an  example 


O'  =./dt[~U(13)]U(23)  +  /dt[|jU(23)]U(13)  .  (A-15) 

We  consider  the  three  possible  successions  of  two  binary  collisions  in 

Figs.  A2,  A3  and  A4  separately.  We  first  evaluate  (A-15)  for  phase  points 

for  which  one  encounters  along  the  free  trajectories  the  two  collisions 
'  '  +) 

indicated  in  Fig.  A2  .  In  the  figure  a  circle  indicates  a  collision  between 

+)  In  diagram  A2  (when  read  from  top  to  bottom)  we  disregard  the  possibility 
that  the  collision  between  2  and  3  may  cause  a  deflected  path  of  3  to  collide 
with  the  original  path  of  1.  The  contribution  of  such  an  event  is  taken  into 
account  by  diagram  A3  (upon  interchanging  1  and  2) . 
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Figure  A  2,  Diagram  of  successive  collisions  leading  to  a  contribution 

U(13)U(23)  to  the  surface  integral. 

the  particles  whose  trajectories  are  enclosed.  In  the  case  under  considera¬ 
tion  both  particles  1  and  3  and  particles  2  and  3  are  aimed  to  collide.  We 
have  indicated  the  free  trajectories  by  dotted  lines  when  they  have  to  be 
distinguished  from  the  actual  trajectories,  if  we  would  take  into  account 
the  interaction  between  1  and  3.  Following  the  particles  along  their  free 
trajectories,  particles  1  and  3  are  inside  each  other  (r^(t)<a)  during  the 
time  interval  t2<t<t^and  particles  2  and  3  are  inside  each  other  during  the 
time  interval  t.<t<t_. 
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For  the  evaluation  of  the  first  integral  in  (A-15)  we  have  to  study 
the  variation  of  U(13)  as  a  function  of  t  along  the  free  trajectories.  For 
t>t^U(13)  has  a  certain  value  which  is  constant,  since  the  operators 
ultimately  operate  on  functions  of  the  velocities  alone.  For  t2<t<t^  the 
operator  fif  (13)  creates  infinite  velocities;  this  contribution  vanishes  as 
a  result  of  the  normalization  of  the  distribution  function  and  U(13)  can 
be  taken  as  -1.  Finally,  for  t<t2>  particles  1  and  3  are  no  longer  aimed 
to  collide  and  U(13)  =  0.  Thus 


Similarly 


t>t. 


U(13)  =  U (13) 


V^I 

:  U (13) 

=  -1 

(A-16) 

t<C2 

:  F (13) 

=  0 

fc>t3 

:  U (23) 

=  U(23) 

t4<t<t3 

:  U (23) 

=  -1 

(A-17) 

t<fc4 

:  U  (23) 

=  0  . 

It  follows  from  (A-16)  that  the  first  integrand  in  (A-15)  only 
contributes  for  t  =  t^  and  t  =  t^  : 

/dt[|^U(13)]U(23)  = 

tj+6  t2+6 

lim  f  dt[f-U(13)]U(23)  +  lim  /  dt [|-U(13) ]U(23)  . 
tj-6  6-H-O  t2-6  ot 

For  hard  spheres  we  can  take  6  arbitrarily  small,  so  that  U(23)  does  not 
vary  during  the  interval  6 .  Thus 

tj+6  t2+6 

/dt[^-U(13)]U(23)  =  U(13)U(23) |  +  U(13)U(23) | 

Ot  r  . 
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The  second  term  in  (A-15)  can  be  integrated  in  a  similar  way  so  that  we 
obtain 


O'  =  U(13)U(23)t  +6  -  U(13)U(23)t  +  U(13)U(23)t  +g  -  U(13)U(23)t 

+  U(23)U(13)t  +6  -  U(23)U(13)t  _fi  +  U(23)U(13)t  +6  -  U(23)U(13)t 
3  3  u4  4 


(A-18 ) 


where  the  indices  denote  the  values  of  t  for  which  the  operators  have  to 
be  evaluated.  From  (A-16)  and  (A-17) 


U(13)U(23)t  +6  =  U (13)U (23)  , 
D(13)l](23)ti_6  .  u(23)  > 

U(13)U(23)t  _6  -  U(23)U(13)t  _{  -  0  , 


so  that  (A-18)  reduces  to 


O’  =  U(13)U (23)  +  U(23)  +  U(13)U(23)  .  .  +  U(23)U(13)..  ^ 

t2+6  t3+5 

-  U(23)U(13)  .  +  U(23)U(13)  _  . 

t„-o  t,+o 


(A-19) 


In  order  to  evaluate  the  remaining  terms  in  (A-19)  we  distinguish 

three  different  cases. 

Case  I:  t,<t-<t„<t 
4  3  2  l 

In  this  case  the  two  collisions  are  spatially  separated .  From  (A-16) 
and  (A-17)  we  conclude 


U(13)U(23)  fi  =  -  U (23)  , 

U (23)U (13) t  +6  =  U(23)U(13)t  _6  =  U(23)U(13)t  _{=  0  , 
3  3  4 
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so  that  (A-19)  becomes 

0*  =  U(13)U(23)  .  (A 

Case  II:  t^t^t-j^ 

In  this  case  the  collision  between  2  and  3  starts  while  1  and  3  are 


still  overlapping.  In  this  case 

1  h  ■ 

U(13)U(23)t  +g  =  +  1 

;  U(23)U(13)t  +g  =  -U(23) ; 

U (23)U(13)  .  =  +  1 

t3-6 

;  U(23)U(13)t  +g  =  0  , 

4 

so  that  (A-19)  reduces  to  (A-20). 

Case  III:  t^t^t^^ 

In  this  case  the  collision  between  2  and  3  not  only  begins  while  1  and 
3  are  overlapping,  but  2  and  3  have  also  separated  before  1  and  3  leave  their 
overlapping  configurations.  We  have 

11(13)11(23)^  =  0  ;  U(23)U(13)t  +g  =  -  U(23)  ; 

U(23)U(13)t  _g  -  +  1  ;  U(23)U(13)t  +g  =  +  1, 

■3.  4 

so  that  (A-19)  reduces  again  to  (A-20) . 

Thus  we  conclude  that  for  phase  points  represented  by  Fig.  A2  the 
integral  (A-15)  leads  to  the  result  (A-20)  regardless  whether  or  not  the  collisions 
are  spatially  separated,  provided  that  the  collision  between  2  and  3  begins  after 
the  collision  between  1  and  3  started,  when  the  diagram  is  retraced  from  top  to 
bottom. 

Using  the  same  method  one  can. prove  that  also 


O'  =  U(13)U(23)  t 


(A-21) 


145,-, 


AEDC-TR-69-68 


for  phase  points  represented  by  the  diagram  of  Fig.  A3.  When  the  integration 
is  carried  out  for  phase  points  represented  by  Fig.  A4,  one  finds^ 

O'  =  0  .  (A-22) 

Thus  events  of  this  kind  do  not  contribute  to  the  surface  integral,  although  ' 
the  operator  U(13)U(23)  itself  would  not  vanish  for  t>tQ. 

The  contribution  of  any  diagram  containing  two  other  successive 
binary  collisions  is  obtained  by  permutation  of  the  particle  numbers.  .  Thus 
for  hard  spheres  (A-14)  reduces  to  the  surface  integral  as  given  by  Green^ 

A  .  J  A 

0(123)  =  /dr21(0)dr31(0)  V^.r^CO)  T(123)  ,  (A-23) 


Figure  A  3.  Diagram  of  successive  collisions  leading  to  a  contribution 

U(13)U(23)  to  the  surface  integral. 
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Figure  A  4.  Diagram  of  successive  collisions  not  leading  to  any  contribution 
to  the  surface  integral,  although  U(13)U(23)  ^  0  for  t>tQ. 

with 

T(123)  =  U(123)  -H  U(<x)U(3)  .  (A- 24) 

a>3 

The  notation  a>3  indicates  that  the  product  U(a)U(B)  should  only  be  evaluated 
for  successive  binary  collisions  involving  pairs  a  and  3  for  which  the 
a-collision  occurs  prior  to  the  B-collision  when  the  diagrams  are  retraced 
backwards  in  time.  For  example,  if  a=13  and  3=23,  U(13)U(23)  should  only  be 
evaluated  for  the  collision  sequences  shown  in  Fig.  A2  and  A3  (regardless,  of 
course,  of  the  occurrence  of  other  collisions). 
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,To  evaluate  the  triple  collision  operator  T(123)  in  more  detail 
we  introduce  the  following  notation.  Each  binary  collision  in  a  collision 
sequence  representing  the  dynamical  history  of  a  phase  point  will  be  in¬ 
dicated  by  listing  the  pairs  whose  interaction  has  to  be  taken  into 
account  successively  for  that  particular  collision  to  occur.  For  example, 
in  Fig.  A2  the  first  collision  is  denoted  by  (13)  and  the  second  collision 
by  (23).  In  Fig.  A3  the  first  collision  is  again  denoted  by  (13),  but  the 
second  collision  is  now  denoted  by  (13) (23),  since  the  collision  between 
2  and  3  occurs  after  the  interaction  between  1  and  3  has  been  taken  into 
account.  With  this  notation  (A-24)  can  be  written  more  specifically  as 

T (123)  -  U(123)  -IE  S (a)U(B)  +  E  E  U(6)  (a*0)  (A-25) 

(cO>(ct)(S)  (a)>(6) 

where  the  summation  in  the  second  term  is  to  be  taken  over  all  pairs  of 
collisions  (a)  and  (a)(6)  for  which  (a)  occurs  prior  to  (a)(6),  and  the 
summation  in  the  third  term  over  all  pairs  of  collision  (a)  and  (0)  such 
that  (a)  occurs  prior  to  (6) . 

The  triple  collision  operator  T(123)  can  be  decomposed  into  a  sum 
of  operators,  each  of  which  is  associated  with  a  particular  triple  collision 
event : 

T (123)  =  E  E  T  (a1;a2)  *  (A-26) 

U  perm  y 

Here  we  have  indicated  the  possible  pairs  of  particles  by  a^  (i  =  1,2,3). 

The  operators  T^(a^;a2)  are  defined  in  Section  2.5  for  the  specific  example 
a^  =  12  and  a2  =  13;  a2  refer  to  the  pairs  of  particles  involved  in  the 
first,  resp.  second  collision.  The  first  summation  in  (A-26)  is  to  be  taken 
over  the  six  diagrams  p  =  Rl,  R2,  Cl,  C2,  HI,  H2  and  the  second  summation 
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over  the  six  permutations  of  the  particles  1,  2  and  3. 

The  validity  of  the  addition  theorem  (A-26)  for  an  arbitrary  phase 
point  can  be  demonstrated  as  follows.  Let  (a^)  be  the  first  real  collision 
in  the  diagram  representing  the  dynamical  history.  If  there  is  a  second 
real  collision,  we  assume  that  to  be  (a^) (a2)  .  As  a  consequence,  (a^)> (a^) (a^)  . 
We  also  have  to  consider  hypothetical  collisions  of  the  type  (a)  and  (a)(6). 

If  there  is  only  one  hypothetical  collision  of  the  type  (a) ,  we  suppose  it  is 
^  there  are  two  hypothetical  collisions  of  this  type,  we  assume 
(°t2^>^a3^*  Then  from  (A-25)  T(123)  can  be  written  as 


T (123)  =  [  § (123)  -/(c^)  - /(a1)U(a2)]  -/(a^U^) 


(a1)> (a^) (a3) 


-  E  E  a)U(6)  +  U(cO,  (a^0) , 
(a)> (a)  (6)  J 


(A-27) 


in  which  the  conditions  concerning  the  time  order  of  the  events  are  indicated 
where  necessary.  The  sum  of  the  first  three  terms  between  the  brackets  is 
only  different  from  zero  if  there  are  three  or  more  real  collisions,  thus 
if  the  conditions  for  either  Rl(a^;a2)  or  Cl(a^;a2)  are  satisfied.  In  that 
case 


S  (123)  -/(a,)  -,T(a1)U(a2)  =/(123)  -JV^eTO^), 

which  is  just  the  contribution  TR1(“1;«2^  or  TCl^al;a2^'  The  next  term  c/V-^U^) 
is  different  from  zero,  if  there  is  a  collision  (c^Ka^)  such  that  («1)>  (a^  (a3) . 
This  implies  also  that  (a^) > (a^) (a2)> (a^) (a3) ,  since  otherwise  (a^) (a3)  would 
be  a  real  collision  contrary  to  our  assumption.  Thus  the  conditions  for  diagram 
are  satisfied  with  contribution  =  “  *^“1^  ^“3)*  Similarly, 

one  easily  verifies  that  the  terms  -  J'(a)U(3),  with  (a)>(a)(6)  and  are 


.149 


AEDC-TR-69-68 


only  different  from  zero  If  the  conditions  for  diagram  R2(a^;a),  when  3=a^, 
or  for  diagram  C2(a^;a),  when  are  satisfied,  leading  to  the  contribution 

T^2 («!»“)  =  -«T(a)U(a1)  or  T^a^a)  =  -  J*(a)U(8).  Finally,  the  last  term 
U(a^)  does  not  vanish  if  there  is  a  collision  (a3).  As  a  consequence  of  our 
assumption  this  implies  (a2)> (a^) ,  so  that  the  conditions  for  diagram  H2(ot^;ci2) 
are  satisfied  with  the  contribution  =  U(a^).  This  concludes  the 

proof  that  T(123)  can  be  represented  as  a  sum  of  operators  associated  with 
the  basic  triple  collision  diagrams,  as  expressed  by  (A-26) . 

The  result  (A-26)  for  T(123)  has  to  be  substituted  into  expression 
(A-23)  for  0(123).  Since  the  integrals  in  (2.4-16)  for  X^  and  are 
symmetric  in  all  three  particles,  it  is  sufficient  to  evaluate  (A-26)  only 
for  one  permutation  of  the  particles  and  to  multiply  the  result  by  6.  Thus, 
we  replace  0(123)  in  (2.4-16)  by 


0(123)  =  6  0(12;13)  =  6  E  /dr^Ojdr^O)  V^.r^O)  T  (12;13)  .  (A-28) 

y  v 

A 

In  (A-28),  r2]/0)  ant*  r^(0)  represent  the  relative  position 
vectors  at  the  time  of  the  first  collision  in  the  diagrams  of  Fig.  4.  Let 
be  the  perihelion  vector  of  the  first  collision,  k2  the  perihelion  vector 
of  the  second  collision,  and  t  the  time  between  the  first  and  the  second 
collision.  To  specify  the  collision  cylinder  associated  with  the  second 
collision,  we  have  to  make  a  distinction  between  the  diagrams  of  the  first 
kind  in  which  this  collision  is  real  and  the  diagrams  of  the  second  kind  in 
which  this  collision  is  hypothetical.  Therefore,  we  define 


V31U)  =  ^-(V^.^)^,  (y  =  Rl,  Cl,  HI), 
V^d)  =  V31,  (y  =  R2,  C2,  H2)  . 


(A-29) 
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Then  V^(l)  is  the  relevant  relative  velocity  to  be  considered  for  the 
second  collision.  In  terms  of  these  new  variables,  (A-28)  transforms  to 

0.(12;  13)  =  S  o4/dk1/dk2  /dTV21.‘k1V21(l).k2Tij(12;13).  (A-30) 

It  should  be  kept  in  mind  that  the  integral  in  (A-23) ,  and  consequently 
in  (A-30) ,  is  to  be  carried  out  over  phase  points  corresponding  to  receding 
particles .  This  implies  that  in  the  evaluation  of  (A-30)  one  has  to  verify 

•A 

for  any  given  k^,  k2 ,  and  t ,  that  the  particles  will  not  suffer  any  collision 
in  the  future. 
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APPENDIX  B.  DYNAMICS  OF  A  COLLISION  BETWEEN  TWO  HARD  SPHERES- 


Consider  two  hard  spheres  1  and  2  with  unit  diameter.  The  relative 
position  of  2  with  respect  to  1  is  r!^  ■=  r^  -  ^  and  the  relative  velocity 

^  ^  m ^  A 

V21  =  v2  “  vi  at  t^ie  initial  time  t  =  0.  The  "impact  vector"  b2^  (see 
Pig.  Bl)  of  the  collision  between  1  and  2  is  related  to  the  initial  position 


and  initial  velocity  v91  by 


21 


T*  (v21  11  r21>  *  V21 

b21  *  - 2 

V21 

and  the  perihelion  vector  k  of  the  collision 
k  =Jb21  +  *21  1,1  ‘  b212  ' 


(B-l) 


(B-2) 


The  two  particles  will  experience  a  collision,  if  they  approach  each  other 


v21'r21  <  0  ■ 


(B-3) 


and  if  the  impact  parameter  b2^  is  smaller  than  the  diameter  (cr  =  1) 


-•»  v2  2 

Cv2i  x  rn)  -  v21  <  0  . 


(B-4) 


The  time  t  at  which  the  collision  will  take  place  is 


t  =  - 


V21*r21  +  ^  ~  b21 
V21 


(B-5) 


The  velocities  v^,  v^  after  the  collision  are  related  to  the  initial  velocities 
vr  v2  by 


=  v1  +  (v21-k)k  ( 


^2  =  ^2  ”  ^21'"^k  * 


(B— 6  ) 
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